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PREPAGE 


This  Memorandum  contains  descriptive  material  and  reference  Infor¬ 
mation  necessary  to  understand  and  use  TACTICS,  a  computer  program 
which  mathematically  simulates  the  flight  trajectories  of  as  many  as 
three  different  vehicles  simultaneously.  The  program  has  been  in  use 
at  Rand  and  elsewhere  for  about  a  year  and  a  half,  principally  in 
studies  of  aircraft  and  missile  performance  in  air-to-air  combat  and  in 
air-to-surface  missile  (ASM)  and  surface-to-air  missile  (SAM)  simula¬ 
tions.  Because  the  simulation  model  is  primarily  useful  as  a  research 
tool  for  studying  interceptor-target  guidance  and  interceptor  trajec¬ 
tories  in  general,  emphasis  has  been  placed  on  providing  versatility 
and  flexibility  for  solving  a  wide  variety  of  problems. 

Magnetic  tape  copies  of  the  program  have  been  sent  upon  request 
to  various  facilities  of  the  USAF  and  USN  and  contractors  engaged  in 
research  in  related  fields,  as  well  as  to  the  armed  forces  (or  affil¬ 
iated  institutions)  of  Japan,  Canada,  and  Germany.  The  descriptive 
material  should  be  of  inter'' to  those  concerned  with  related  fields 
of  effort,  while  the  reference  information  is  necessary  for  a  thorough 
understanding  and  effective  use  of  the  model. 


SUMMARY 


TACTICS  is  a  computer  program  written  in  FORTRAN  IV  which  mathe¬ 
matically  simulates  the  dynamics  of  flight  in  three-dimensional  space 
of  as  many  as  three  vehicles  simultaneously.  The  purpose  of  this  Mem¬ 
orandum  is  to  acquaint  prospective  users  with  the  capabilities  and 
basic  theory  of  the  program  and  to  serve  as  a  reference  manual  for 
those  who  wish  to  use  the  program. 

The  first  part  of  the  Memorandum  contains  the  description  and 
theory  of  operation  and  is  oriented  toward  those  with  a  mathematical 
or  technical  background.  The  second  part  is  concerned  with  how  to  use 
the  program,  i.e.,  how  to  provide  input  data,  select  options,  and  de¬ 
velop  a  flight  program.  Wherever  possible,  FORTRAN  symbols  are  rele¬ 
gated  to  this  second  part.  A  number  of  illustrative  examples  of  a 
wide  variety  of  problems  are  given  in  detail,  including  data  and 
FORTRAN  listings,  in  order  to  facilitate  the  use  of  TACTICS  and  the 
obtaining  of  results  without  detailed  knowledge  of  the  inner  workings 
of  the  program.  In  many  cases  it  should  be  possible  to  set  up  a  spe¬ 
cialized  problem  by  modifying  or  combining  various  features  of  the 
examples. 
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SYMBOLS 


FORTRAN 

Notation 

Symbol 

Definition 

A(I,J) 

^A 

Unit  vector  normal  to  velocity  vector 

V  and  also  in  the  horizontal  plane 

AB00ST(I) 

Boost  acceleration^  assumed  to  he  an 
average  value  (ft/sec^) 

AC0M(I) 

Absolute  magnitude  of  vehicle  lateral 
acceleration  (ft/sec^,  printout  in  g*s) 

AC0MA(I) 

^Ch 

Horizontal  component  of  lateral  accelera¬ 
tion  (g's) 

AC0MD(I) 

^CV 

Vertical  component  of  lateral  acceleration 
of  vehicle (1)  (g's) 

ALPHA(I) 

a 

Angle  of  attack  (deg) 

ALPHAO(I) 

a 

o 

Zero-lift  angle  of  attack  (deg). 

ALT(l) 

h 

Altitude  of  the  vehicle  (ft) 

A0UT<I) 

a 

0 

Absolute  magnitude  of  output  lateral 
acceleration  of  vehicle  (i)  (g's) 

A0UTA(I) 

®ch 

Horizontal  component  of  output  lateral 
acceleration  of  vehicle  (1)  (g's) 

A0UTD(I) 

^oV 

Vertical  component  of  output  lateral 
acceleration  of  vehicle  (i)  (g's) 

AREA(I) 

A 

2 

Reference  area  (ft  ) 

ASMAX(I) 

Smax 

Structural  lateral  acceleration  limit 

(g'8) 

AZMAX(I) 

n 

max 

Maximum  azimuth  angle  (gimbal  limit) 

(deg) 

AZMUTH(l) 

^12 

Azimuth  angle  in  aircraft  coordinates, 
vehicle  2  with  respect  to  vehicle  1 
(deg,  +  for  right,  -  for  left) 

Ald.J) 

\l 

Unit  vector  normal  to  the  axis  of  the 
aircraft  and  also  in  the  horizontal  plane 

BANK(l) 

^B 

Aircraft  bank  angle  defined  in  rela¬ 
tion  to  wind  axis  system  (see  Figs.  7 
and  8)  (deg) 

BC0N(I) 

dC  /d(cf) 

Coefficient  used  with  parabolic  approxi- 

matlon  for  drag  coefficient  at?  a  function 
of  lift  coefficient 


Subscripts  "i"  which  indicate  vehicle  numbers  have  been  omitted 
throughout  the  list  for  notational  convenience. 
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FOmiUC 


Notation 

S/aibol 

Definition 

EE^Uti^GCl) 

»12 

Bearing  angle  in  aizcT-afc  coordinates, 
vehicle  2  with  respect  to  vehicle  1  (deg'^ 

3EASNG(2) 

®13 

Bearing  angle  in  aircraft  coordinates, 
vehicle  3  vith  respect  to  vehicle  1  (deg) 

BEAUiGO) 

®23 

Bearing  angle  in  aircraft  coordinates, 
vehicle  3  vith  respect  to  vehicle  2  (deg) 

SEAX!iG(») 

®21 

Bearing  angle  in  aircraft  coordinates, 
vehicle  1  vith  respect  to  vehicle  2  (deg) 

B£A5NG(5) 

®31 

Bearing  angle  in  aircraft  coordinates, 
vehicle  1  vith  respect  to  vehicle  3  (deg) 

BEASNG(6) 

®32 

Bearing  angle  in  aircraft  coordinates, 
vehicle  2  vith  respect  to  vehicle  3  (deg) 

B£IA(I) 

S 

Ballistic  coefficient  (ib/ft^) 

CDOC^Nd) 

r* 

''D 

Zero-lift  drag  coefficient 

CLMAX(I) 

o 

^Laax 

Maxlaam  aerodynaaic  lift  coefficient 

C0DSAG(I) 

Sjf 

Aerodynaaic  drag  coefficient 

C0LIFT(I) 

C 

L 

Aerodynaaic  lift  coefficient 

CWD#r(I) 

« 

W 

Tiae  rate  of  change  of  weight  (Ib/sec) 

DCI.J) 

h 

Unit  vector  norssal  cojthe  velocity  vector 
and  noraal  to  vector  1^,  forsing  a  right- 
hand  systea 

DATA(200) 

Initial-conditiOT  input  deta  anglng  froo 
I  to  200  (so  far,  only  143  ai  used) 

DELV 

AV 

Boost  velocity  of  aisslle  (ft/sec) 

DEH3(I) 

P 

Air  density  (slug/ft^) 

DKAG(I) 

D 

Drag  force  (lb) 

DTMIN 

Kininiina  value  for  integration  step  size, 
used  only  for  deteraination  of  alss  dis¬ 
tance  (sec) 

DT0 

Starting  value  fur  integration  step  size 
(sec) 

DTP0 

Time  interval  for  printing  output  (sec) 

DVPHI(I) 

Ay 

Assumed  error  in  y  aiming  (deg) 

DVTU(I) 

AO 

Assumed  error  in  for  aiming  (deg) 

Dld.J) 

^1 

Unit  vector  normal  to  the  axis  of  the 
aircraft  and  also  in  a  plane  normal  to 

the  horizontal  plane 
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FORTRAN 

Notation 

Synbol 

Definition 

£LEV(1) 

Elevation  angle  in  aircraft  coordinates, 
vehicle  2  with  respect  to  vehicle  1  (deg, 
+  for  above,  -  for  below) 

ELEV(2) 

Elevation  angle  in  aircraft  coordinates, 
vehicle  3  with  respect  to  vehicle  1  (deg, 
+  for  above,  -  for  below) 

EEEVO) 

^3 

Elevation  angle  in  aircraft  coordinates, 
vehicle  3  with  respect  to  vehicle  2  (deg, 
+  for  above,  -  for  below) 

ELEy(4) 

Elevation  angle  In  aircraft  coordinates, 
vehicle  1  with  respect  to  vehicle  2  (deg, 
+  for  above,  -  for  below) 

£LEV(5) 

Elevation  angle  in  aircraft  coordinates, 
vehicle  1  with  respect  to  vehicle  3  (deg, 
+  for  above,  -  for  below) 

ELEV(6) 

^32 

Elevctlon  angle  In  aircraft  coordinates, 
vehicle  2  with  respect  to  vehicle  3  (deg, 
+  for  above,  -  for  below). 

ELEVMXd) 

^■ax 

Maxinua  elevation  angle  i%iabal  liid.ts) 
(deg) 

EXTR(2G> 

Extra  qucntltles  (real)  In  C^9#K  for 
optional  use  (the  first  six  are  part  of 
standard  printout) 

G 

8 

2 

Hass  conversion  (32.174  ft/sec  );  also 
used  as  constant  gravitational  attraction 
for  flat-earth  option 

GF0RC 

Total  lateral  acceleration,  including 
gravitational  effects  specified  for  a 
maneuver  (g'a) 

GF^RCE(I) 

F 

n 

Total  lateral  acceleration,  including 
gravitational  effects  (g*s) 

HM 

Actual  Integration  step  size  used  by  the 
program  (sec) 

HKIM 

Minimum  specified  value  for  Integration 
step  size  used  In  variable  Adams-Moulton 
mode  Integration  (sec) 

HMX 

Maxiisum  specified  value  for  Integration 
step  size  used  In  variable  Adama-Moulton 
mode  Integration  (sec) 

IAERi$ 

Flag  Indicating  aerodynamic  option 

lEXTRA(lO) 

Ten  extra  integer  quantities  for  optional 

use  (C0MH0N  package) 
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FORTRAN 

Notation 

Syo^ol 

Definition 

ILAUN 

Flag  indicating  sequence  of  steps  Involved 
in  Missile  launch  (see  Appendix  H) 

IHISS 

Flag  indicating  that  the  prograa  should 
continue  after  finding  point  of  closest 
approach  (aiss  distance) 

LMPLSE(I) 

T 

Specific  iapulse  of  rocket  aotor  (sec) 

INERF,  INERT 

Flag  indicating  that  a  vehicle's  velocity 
is  in  relation  to  a  rotating  or  non¬ 
rotating  (inertial)  coordinate  systea 

EPSL^ 

Threshold  value  used  in  on-off  control 
lavs  for  stability 

ERTEST 

Maxiaua  a^ovabie  relative  truncation 
error  for  Adaas-Houlton  variable  step 
size  integration  node 

IPRINT(20) 

Flag  indicating  printout  option 

ZR#T8 

Flag  indicating  option  for  rotating  or 
nonrotatlp^  earth 

IST0RE 

Flag  indicating  that  position  and  ve¬ 
locity  values  are  to  be  stored  at  the 
tine  of  launch 

ITAU(I) 

Flag  indicating  nuaber  of  vehicle  first- 
order  tiase  lags 

ITMR 

Flag  indicating  thrust  option 

JAIM0S 

Flag  specifying  whether  initial-con¬ 
dition  value  of  velocity  ie  expressed 
in  Mach  nuaber  or  ft/sec 

JINTEG 

Flag  specifying  integration  aode 

JP«}L,  KP0L,  LP0L, 
NP0L,  NP0L 

Flags  used  in  F0L1CY  subroutine 

JVEH(I) 

Flag  Indicating  aerodynaaic  option  for 
tables 

KINTEG 

Flag  indicating  whether  a  round-earth 
or  flat-earth  option  is  to  be  used 

KLAUN 

Deciaal  fraction  of  aissile's  aaxiaum 
range  at  which  it  is  to  be  launched 

LAMDAO(I) 

Navigation  constant  for  closed-loop 
guidance  routines 

LAT(I) 

Geocentric  latitude  of  the  vehicle  (deg) 

LATO 

% 

Latitude  of  the  local  coordinate  system 
origin  (deg; 

FORTRAN 

Notation 

SyrAol 

Definition 

LEVEL(I) 

Flag  set  in  subroutine  STRLVL  tised  to 
coaDunicate  to  P0LICT  that  the  vehicle 
velocity  vector  is  in  a  'horizontal  plane 
within  a  tolerance  of  0.002  rad 

LIFT(I) 

L 

Aerodynaalc  lift  force  (lb) 

A 

Longitude  of  the  vehicle  (deg) 

L^NGO 

h 

O 

longitude  of  the  origin  of  the  local 
coordinate  syatea  (deg) 

JtXCHd) 

H 

Mach  Duaber 

HACHMX(I) 

M 

max 

Placard  Halt  (aaxiaue  Mach  nuid>er  for 
vehicle) 

M1S3<R 

Missile  range  to  target  within  which 
progran  will  autocaticall}*  initiate 
process  for  alas  distance  coi^utatian 
(ft) 

H^E 

Flag  indicating  captive  flight  optlcn 

NPRINT 

Flag  Indicating  auaber  of  sections  of 
output  t($  be  printed 

#HEGA{I,J) 

u 

Angular  rate  vector  of  vehicle  (i)  with 
respect  to  x,  y«  z  coordinate  fraae 
(rad/sec) 

^HEGAB(1,J) 

“b12 

Angular  rate  bias  tem  used  for  predic¬ 
tive  guidance,  vehicle  2  with  respect 
to  vehicle  1  (rad/sec) 

AflEGAB(2,J) 

“b13 

Angular  rate  bias  tem  used  for  predic¬ 
tive  guidance,  vehicle  3  with  respect 
to  vehicle  1  (rad/sec) 

^HEGAB(3,J) 

“B23 

Angular  rate  bias  tem  used  for  predic¬ 
tive  guidance,  vehicle  3  with  respect 
to  vehicle  2  (rad/sec) 

0i4EGAB(4,J) 

“B21 

Angular  rate  bias  tem  used  for  predic¬ 
tive  guidance,  vehicle  1  id.th  respect 
to  vehicle  2  (rad/sec) 

0M£GAB(5,J) 

“b31 

Angular  rate  bias  tem  used  for  predic¬ 
tive  guidance,  vehicle  1  with  respect 
to  vehicle  3  (rad/scc) 

^1EGAB(6,J) 

“b32 

Angular  rate  bias  teno  used  for  predic¬ 
tive  guidance,  vehicle  2  with  respect 
to  vehicle  3  (rad/sec) 

0MEGAE 

(J 

A 

Earth's  angular  rate  of  rotation 

(7.29211585  x  10-5  rad/aec) 


-xviii- 


FORTRAN 

Notation 

Syod>ol 

Definition 

^1EGAR(1,J) 

®R12 

Relative  angular  rate  vector,  vehicle  2 
with  respect  to  vehicle  1  (rad/sec) 

tfKEGAR(2,J) 

“’Ria 

Relative  angular  rate  vector,  vehicle  3 
with  respect  to  vehicle  1  (rad/sec) 

^MEGARO.J) 

“r23 

Relative  angular  rate  vector,  vehicle  3 
with  respect  to  vehicle  2  (rad/sec) 

^!EGAR{4,J) 

"r21 

Relative  angular  rate  vector,  vehicle  1 
with  respect  to  vehicle  2  (rad/sec) 

^1EGAR(S,J) 

“r31 

Relative  angular  rate  vector,  vehicle  1 
with  respect  to  vehicle  3  (rad/sec) 

^ffiGAK(6,J) 

“r32 

Relative  angular  rate  vector,  vehicle  2 
with  respect  to  vehicle  3  (rad/sec) 

PHID^(I) 

V 

Tiiae  rate  of  change  of  <p  (rad/sec) 

PRES (I) 

P 

Air  pressure  (Ib/ft^) 

Q(I) 

q 

Dynanlc  pressure  (Ib/ft  ) 

R(I,J) 

r 

Range  vector  from  x,  y,  z  origin  to 
vehicle  (1)  (ft) 

RAD 

ir/180'“ 

Factor  for  converting  degrees  to 
radians  (rad/deg) 

RD^Td) 

^12 

Range  rate,  vehicle  2  with  respect  to 
vehicle  1  (ft/sec) 

RD0T(2) 

'13 

Range  rate,  vehicle  3  with  respect  to 
vehicle  1  (ft/sec) 

RD0T(3) 

• 

^23 

Range  rate,  vehicle  3  with  respect  to 
vehicle  2  (ft/sec) 

RD0T(4) 

^21 

Range  rate,  vehicle  1  with  respect  to 
vehicle  2  (ft/sec) 

RD0T(5) 

^31 

Range  rate,  vehicle  1  with  respect  to 
vehicle  3  (ft/sec) 

RD0T(6) 

^32 

Range  rate,  vehicle  2  with  respect  to 
vehicle  3  (ft/sec) 

PJAUN 

Range  at  which  missile  is  to  be  launched 
used  in  LEADCL  subroutine  (ft) 

BMIHAX 

Maximum  range  of  missile  (ft) 

R^LL(I) 

Aircraft  roll  angle  (deg)  defined  in 
relation  to  aircraft  axes  coordinate 
system 

R0LLL 

Total  angle  through  which  vehicle  Is  to 
roll  (deg) 

FORTRAN 

Notation 

Symbol 

Definition 

R0LLR8(Z) 

Time  rate  of  change  of  aircraft  roll 
angle  (deg/sec) 

RR(I,J) 

R 

Geocentric  range  vector  (ft)  directed 
radially  from  the  earth's  center 

RRD0T(1,J) 

R 

Velocity  of  the  topocentric  range  vector 
with  respect  to  Inertial  space  (ft/sec) 

RRELd.J) 

^12 

Relative  range  vector,  vehicle  2  rela¬ 
tive  to  vehicle  1  (ft) 

RREL(2,J) 

Relative  range  vector,  vehicle  3  rela¬ 
tive  to  vehicle  1  (ft) 

RREL(3,J) 

^^23 

Relative  range  vector,  vehicle  3  rela¬ 
tive  to  vehicle  2  (ft) 

RREL(4,J) 

*^21 

Relative  range  vector,  vehicle  1  rela¬ 
tive  to  vehicle  2  (ft) 

RREL(5,J) 

*^31 

Relative  range  vector,  vehicle  1  rela¬ 
tive  to  vehicle  3  (ft) 

RREL(6,J) 

*^32 

Relative  range  vector,  vehicle  2  rela¬ 
tive  to  vehicle  3  (ft) 

RO 

Average  radius  of  the  earth  (20,902,287  ft) 

RO(I) 

^0 

Position  of  the  local  x,  y,  z  coordinate 
system  origin  with  respect  to  inertial 

« 

space  (ft) 

ROD0T(I) 

"o 

Velocity  of  the  local  x,  y,  z  coordinate 
system  with  respect  to  inertial  space 
(ft/sec) 

SGAMA(I) 

0 

Absolute  magnitude  of  the  angle  between 
the  velocity  vector  of  the  Interceptor 
and  the  line  of  sight  from  interceptor 
to  target  (rad) 

1,L0PE(I) 

<1C^  /da 

Slope  of  versus  a  curve 

S0UND(I) 

3 

Speed  of  sound  (ft/sec) 

TABC0N(I) 

T 

ab 

Constants  to  be  used  for  afterburner 
thrust  (lb) 

TAU(I,J) 

T 

First-order  approximation  for  overall 
missile  response  time 

TBURNl 

First-stage  burning  time  (sec) 

TBURN2 

Second-stage  burning  time  (sec) 

TEMP(I) 

T. 

Temperature  (degrees  Fahrenheit) 

FORTRAN 

Notation 

Symbol 

Definition 

TGUIDE(I) 

Time  interval  missile  is  to  fly  un¬ 
guided  after  launch  (sec) 

THC0N(I) 

T 

M 

Constants  to  be  used  for  military  thrust 
(lb) 

THD0T(I) 

0 

Time  rate  of  change  of  theta  (rad/sec) 

THR0TL(I) 

K 

Constant  used  to  multiply  thrust  to 
represent  throttle  setting 

THRUST (I) 

T 

Propulsive  thrust  force  (lb) 

TIME 

t 

Running  time  (sec) 

TLAUN(I) 

Launch  time  (sec) 

T0TAL 

Time  value  at  which  program,  is  to  stop 

UNITE (I, J) 

Unit  vector  directed  along  the  lift 
vector 

UNITLL(I,J) 

^LL 

Unit  vector  directed  eastward  along  the 
local  parallel  of  latitude 

UNITPP(I,J) 

^P 

Unit  vector  directed  northward  along 
the  local  meridian  of  longitude 

UNITPVd.J) 

h 

£nit  vector  normal  to  velocity  vector 

V  and  along  the  net  lateral  acceleration 
vector 

UNITR(K,J) 

\  ■ 

Unit  vector  directed  along  the  relative 
range  vectors  and  using  the  same  sub- 
ccrlpt  notation 

UNITRRd.J) 

Unit  vector  directed  radially  fr£m  the 
earth's  center  along  the  vector  R 
(RR(I,J)) 

UNITTd.J) 

^T 

Unit  vector  directed  along  the  longi¬ 
tudinal  axis  of  the  vehicle  (1)  and  also 
a8sumed_to  be  coincident  with  the  thrust 
vector  T 

UNITV  d,J) 

Unit  vector  corrected  along  velocity 
vector  V 

Vd,J) 

V 

• 

Velocity  vector  (ft/sec) 

VU0Td) 

V 

Rate  of  acceleration  or  deceleration 
along  the  trajectory,  l.e.,  changes 

In  speed  (ft/sec^) 

VREL(1,J) 

^12 

Relative  velocity,  vehicle  2  with  re¬ 
spect  to  vehicle  1  (ft/sec) 

VREL(2,J) 

^]3 

Relative  velocity,  vehicle  3  with  re¬ 
spect  to  vehicle  1  (ft/sec) 

Notation 


Syabol 


Definition 


VREL(3,J) 

^23 

Relative  velocity,  vehicle  3  with 
spect  to  vehicle  2  (ft/aec) 

re- 

VREL(4,J) 

^21 

Relative  velocity,  vehicle  1  with 
spect  to  vehicle  2  (ft/sec) 

re- 

VREL(5,J) 

^31 

Relative  velocity,  vehicle  1  with 
spect  to  vehicle  3  (ft/sec) 

re- 

VREL(6,J) 

'''32 

Relative  velocity;  vehicle  2  with 
spect  to  vehicle  3  C ft/aec) 

re- 

WBURN(I) 

Weight  of  missile  at  burnout  (lb) 

WD0T(I) 

W 

Time  rate  of  change  of  weight  (Ib/sec) 

WEIGHT(I) 

W 

Weight  (lb) 

WO(I) 

w 

0 

Initial  weight  of  vehicle  (1)  (lb) 

r 

C' 

& 


'f 


Part  1 


DESCRIPTION  AND  THEORY  OF  OPERATION 


1 


I.  INTRODUCTION 


TACTICS  is  a  computer  program  i#rltten  In  FORTRAN  for  use  in  simu¬ 
lating  the  kinematics  and  dynamics  of  motion  of  three  vehicles  in  three- 
dimensional  space.  It  was  developed  primarily  as  a  research  tool  for 
use  in  detailed  explorations  of  the  mechanics,  geometry,  and  vehicle 
performance  characteristics  of  an  interceptor-target  engagement.  The 
output  is  a  step-by-step  time  history  of  variables  relating  to  each  of 
three  vehicles'  position,  velocity,  acceleration,  applied  forces,  atti¬ 
tude  or  orientation,  and  aerodynamics.  While  the  program  is  highly 
versatile,  its  most  important  capabilities  relate  to  interceptor-target 
guidance  and  intercept  trajectories  in  general.  Since  the  flights  of 
three  vehicles  may  be  represented  simultaneously,  the  program  can  be 
used  to  simulate  aerial  combat  between  aircraft  (e.g.,  a  two-on-one 
engagement,  or  a  one-on-one  with  missile  launch).  There  are  also  a 
number  of  other  possibilities,  such  as  (1)  using  one  or  more  of  the 
vehicles  to  represent  an  alr-to-surface  missile  (ASM) ,  surface-to-air 
missile  (SAM) ,  or  surface-to-surface  missile  (SSM) ,  (2)  using  the  tar¬ 
get  to  represent  an  ICBM  reentry  vehicle  (RV) ,  (3)  using  vehicles  1  and 
2  to  represent  first-  and  second-stage  boosters,  or  (4)  using  one  ve¬ 
hicle  to  represent  an  orbiting  satellite. 

So  far,  the  model  has  been  used  primarily  in  simulating  fighter- 
versus-flghter  and  raisslle-versus-fighter  duels.  Because  of  the  number 
of  maneuver  routines  created  for  this  specialized  purpose,  the  model's 
present  development  in  this  area  is  farther  advanced  than  in  other 
areas.  However,  ASM  and  SAM  simulations  are  currently  being  performed 
in  connection  with  other  Rand  projects,  and  the  program  has  been  used 
successfully  in  satellite-intercept  problems.  The  model's  usefulness 
and  potential  capabilities  are  expected  to  grow  with  each  new  applica¬ 
tion. 

In  the  organization  of  TACTICS,  primary  emphasis  was  placed  on 
versatility  in  order  to  accommodate  a  broad  spectrum  of  problems.  At 
the  same  time,  it  was  thought  that  the  program  should  be  easy  to  use 
and  easily  adaptable  to  refinements  and  new  features.  Accordingly,  it 
was  built  in  modular  building-block  form,  with  many  subroutines  that 
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cc  be  replaced,  modified,  or  simply  disregarded  at  the  option  of 
the  us nr.  In  fact,  new  building  blocks  are  welcomed,  since  they  in¬ 
crease  the  program's  potential  problem-solving  capability.  The  con¬ 
struction  of  TACTICS  is  in  many  respects  comparable  to  that  of  the 
(1) 

ROCKET  program,  although  the  two  have  different  purposes.  Many 
ideas  were  borrowed  from  ROCKET  in  the  formulation  of  TACTICS,  espe¬ 
cially  its  most  importanc  distinguishing  feature,  that  of  allowing  the 
researcher  to  obtain  results  without  a  detailed  knowledge  of  the  inner 
workings  of  the  program. 

The  program  is  able  to  simulate  the  flight  or  trajectory  of  almost 
any  type  of  vehicle  in  almost  any  mixture  of  vehicles.  Various  guidance- 
law  subroutines  are  available  for  simulating  terminal-homing  or  command- 
guided  trajectories  (e.g.,  biased  proportional  and  proportional  naviga¬ 
tion,  lead  collision,  and  pursuit  and  lead  pursuit).  A  variety  of  open- 
ended  control-law  subroutines  are  also  available  for  simulating  aero¬ 
dynamic  maneuvers  such  as  turning,  diving,  climbing,  and  combinations 
thereof.  The  library  of  guidance-  and  control-law  subroutines  has  been 
growing  and  will  probably  continue  to  grow  as  new  problems  are  encountered. 
In  fact,  experience  has  shown  that  it  is  convenient  to  incorporate  into 
one  of  these  control-law  subroutines  the  unique  features  of  a  particular 
vehicle  or  problem  (e.g..  Sidewinder  or  a  hypothetical  aircraft  or  mis¬ 
sile  design) . 

Options  are  available  for  considering  the  earth  either  flat  or 
round  and  either  fixed  or  rotating,  but  the  gravitational  field  is  a 
simple  inverse  square  law  field.  Each  of  the  three  vehicles  has  three 
wind  axes  associated  with  its  attitude  or  orientation  reference  system. 

The  vehicles  may  be  either  fixed  or  in  motion  with  respect  to  the  earth 
but  not  below  its  suvface.  The  miss  distance  or  point  of  closest  ap¬ 
proach  between  two  of  three  vehicles  may  be  calculated  and  the  problem 
run  automatically  terminated  if  so  desired. 

As  mentioned  earlier,  the  step-by-step  time  history  of  the  engage¬ 
ment  is  limited  to  the  consideration  of  the  performance  of  three  vehicles 
at  one  time.  However,  two  devices  have  been  incorporated  for  extending 
this  capability  to  more  than  three  vehicles  by  means  of  sequential  com¬ 
putations.  These  devices,  designated  "recall"  and  "restore,"  are 


-5- 


illustraced  In  Fig.  1.  Consider  a  vehicle  1  which  launches  a  missile 
2  at  a  target  3.  At  some  subsequent  time  or  event  (e.g.y  a  miss)  it 
may  be  desirable  to  recall  vehicle  2  and  place  it  in  captive  flight 
on  either  1  or  3  for  a  subsequent  quasi-vehicle  4.  (It  is  not  neces¬ 
sary  that  4  have  the  same  characteristics  as,  or  even  resemble,  2.) 

To  illustrate  the  restore  feature,  consider  the  same  example  except 
that  at  some  time  or  event  subsequent  to  the  launching  of  vehicle-  2 
(e.g.,  a  hit,  miss,  or  ground  Impact)  we  wish  to  restore  the  launch- 
time  situation.  After  restoration  has  occurred,  events  may  proceed, 
and  a  new  launch  (4,  5,  etc.)  may  take  place  at  a  subsequent  time  or 
event.  Or  perhaps  branching  is  desired,  i.e.,  some  characteristic  or 
parameter  is  altered  and  4,  5,  etc.  are  to  be  launched  under  the  same 
restored  conditions  of  time  and  geometry. 

Consider  the  factors  vital  to  any  program  for  simulating  flight 
trajectories.  As  Illustrated  in  Fig.  2,  it  is  necessary  to  (1)  read 
in  data  as  initial  conditions,  (2)  calculate  the  geometry,  (3)  specify 
the  applied  forces  involved,  (4)  integrate  the  equations  of  motion,  and 
(5)  output  the  answers.  This  is  the  basic  framework  of  TACTICS,  to 
which  embellishments  are  added  (e.g.,  model  atmosphere,  coordinate  trans¬ 
formations,  and  aerodynamic  computations).  Although  several  options  are 
provided  for  modes  of  integration  and  output  form,  the  framework  may 
be  considered  Inflexible  except  for  specifying  the  applied  forces.  In 
this  case,  maximum  flexibility  is  provided  by  an  arrangement  comparable 
to  plug-in  modules;  that  is,  the  applied  forces  and  the  tlme(s),  event(s), 
or  situatlonCs)  dictating  their  application  are  contained  within  two 
subroutines  which  may  be  specialized  to  deal  with  a  particular  problem. 

To  illustrate,  each  particular  problem  is  defined  by  initial-condition 
input  data  and  by  a  number  of  P0LICY  statements,  which  are  logical  ex¬ 
pressions  dictating  the  control  laws  governing  the  flight  of  each  vehi¬ 
cle.  They  are  usually  conditionally  based  on  time  or  on  geometric  and 
kinematic  relationships.  (The  POlLICY  subroutine  is  discussed  in  more 
detail  in  Section  II.) 

A  rudimentary  description  of  the  operational  principles  necessary 
to  simulate  an  intercept  problem  will  now  be  useful.  As  a  starting 
point,  imagine  a  vehicle  in  three-dimensional  space  having  a  position 


aenS  a  vieldci'2}'  iSefisaj  fifficial-cscBflfCioa  fc^sc  tSaca.  If  ao  acod- 
cra!£fe(&*  or  a^lSie::^  forces  are  flaeolwo,  e&e  slaue  bSsccsy  or  cnjee- 
tces7  of  S&e  vsSdcle  vrili  o&vaoasljr  &e  a  seraSg^s:  ULae.  Eowrner.  la 
ebe  gesMsal  ccajeca^sv  sase.  ebere  vill  be  a  cec  ancelecarfsa  Cz,  j,  z) 
(&3e  tro  grsTieaciosal.  aerograaeac,  ad  ebnsr  forces,  ubere  sbe  cosvtae^ 
foadeatcm  bcSaMcs  zSue  retalsaac  force  as2  z,  j,  z  is  f  «'  aa. 

If  f C  is  asnoed  e&*£  all  foroes  caa  S:«  defitsed  ad  specified,  she 
Krajeccezy  slsalacfsa  {problem  is  re&sced  zo  ill  isz^^sasios  c&e  set 
acceleraCiaa  as  a  ferae rfea  of  sfae  co  ebsaia  velsciB^  ad  (2^  inEeg^'^'' 
isf  velociCj'  CO  flbeaia  ptssicioa.  If  ^ae  spiled  aec  forc.%  were  coa- 
scaac  or  a  simple  ferae rigo  of  cise,  coe  czajeecosv  siamalacioo  wscsld  Ese 
scraisseforward  aad  relaclvelj’  simple.  la  e!se  seaccal  ease,  parcicsalarlj 
zbea  clooed-loop  (aidasce  aad  aecodjsamics  are  iisvolved,  the  forces  mav 
be  otiaplieated  foogticas  of  posicioo,  velociey,  geamegry,  ^od  cixe  (aad 
aajp  iawoire  cbe  b^tarlo;;  or  predicted  Idisrior  of  sone  ocber  Tehicle(s»; 
bence  saaerical  ioccgzacioa  cechaiqses  owsc  Emk  used  zo  calctalacc  die 
crajecCczj  stepvise  tsslog  cine  iacreaeac  £c.  It  should  he  clear  freas 
Che  preceding  discossica  chat  Co  start  a  problem  zua,  iapat  data  far 
isicial  posiCioQ  and  velocity  nusC  be  supplied.  The  hulk  of  all  ocher 
ii^mC  data  vill  pertain  to  psraneters  associated  vich  the  definicica 
or  calculation  of  the  forces  to  be  applied  (e.g.,  lift,  drag,  thrust, 
and  gravity).  Representing  the  notion  of  three  different  vehicles 
nultipUes  the  ii^iat  requiresents,  of  course,  but  net  necessarily  by 
a  factor  of  three.  The  definition  of  the  intercept  problem,  in  terss 
of  logic,  ti&ing,  and  geonetric  events,  or  cf  any  one  of  these,  dictates 
vhat  types  of  ccctrcl  lavs  are  to  be  applied  under  which  circumstances. 
This  definition  is  contained  in  a  POLICY  subroutine  that  lit  a  sense  is 
also  an  input  to  the  program.  The  forces  or  force  functions  which  are 
to  be  called  upon  by  POLICY  arc  contained  in  the  library  of  control-law 
subroutines.  It  is  not  expected  that  this  library  vill  ever  include 
every  conceivable  force  function  or  guidance  law,  especially  since 
TACTICS  is  a  research  tool  for  experimentation.  However,  experience 
has  shown  that  new  subroutines  may  usually  be  conveniently  generated 
by  modifying  those  already  on  hand. 
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The  P0iS£:¥  porcioo  cf  ebe  profraaD,  is  sspplied  cbe  user, 

coaslscs  of  logical  searwreass  ani  expressless  diccaElss  the  guldance- 
or  Lcaczol-lav  sohreutiaes  citac  govexs  the  flight  of  each  vehicle. 

These  sc&Eoaclacs  defise  the  varioos  forces  that  are  to  be  applied  to 
a  vehicle  to  perfora  a  certain  aaaecver  or  to  guide  In  accordance  vith 
seme  prescribed  gaidance  Imr.  Ihc  result  of  all  forces  vill  be  a  net 
accelerafcise.  9ce  need  not  be  fan! liar  with  the  ■athcnatics  Involved 
Ic  guidance,  aercdvaanics.  propulsion,  etc.  to  write  a  F0LICT  routine. 
Ssvevec,  a  large  somber  of  options  are  available,  and  the  user  should 
be  fanillar  vith  toem  and  vltb  the  library  of  available  guidance-  or 
control-law  roctlses.  To  illustrate  the  use  of  a  F0L1CT  siAroutlne, 
consider  the  problex  of  sinolatinf.  the  flight  of  an  aircraft  fron  take¬ 
off  to  lasdlsg.  Clearlj,  a  nusber  of  control  lavs  defining  the  forces 
and  hence  the  net  acceleratlcn  would  be  needed  for  even  the  most  ele¬ 
mentary  flight  plan.  This  subject  is  discussed  in  further  detail  in 
Sections  IV  and  Ix.  Consider  the  net  accelerations  associated  with 
each  of  the  followliig: 

o  Clinb  (vertical  plane). 

c  Straight  flight,  i.e.,  no  turning  coiq>onent  of  acceleration. 

o  Turn  (horizontal  plane). 

o  Dive  (vertical  plane). 

Moreover,  assume  that  subroutines  are  available  fpr  describing  the 
above  trajectories  (CLIMB,  STRFLT,  RTURN,  LTURK,  and  DIVE,  respectively). 
A  P0LICY  routine  is  required  to  define  the  conditions  dictating  the 
transition  from  cne  maneuver  to  another.  For  example,  CLIMB  may  be 
called  for  a  given  time  interval  or  until  some  specified  altitude  or 
other  criterion  is  reached;  then  STRFLT  is  called,  and  so  on. 

With  respect  to  the  preceding  illustration,  the  most  obvious 
example  of  choice  of  option  is  the  selection  of  the  CLIMB,  STRFLT,  etc. 
routines.  In  the  argument  listings  for  these  routines  there  are  also 
options  pertaining  to  the  following  conditions: 
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o  The  vehicle  the  lau  is  to  govern  (1,  2,  or  3). 

o  T1:e  angnicude  (and  perhaps  direction)  of  the  propulsive  or 
thrust  force. 

o  The  nuaber  of  g's  coaeanded  for  the  clinb,  dive,  or  turn 
aancuvers. 

a  The  use  of  tabulated  values  or  analytic  functions  in  calcu¬ 
lating  aerodynanic  forces. 

Appendix  D  contains  a  list  of  optional  subroutines  currently  available 
and  instructions  for  their  use.  Section  IX  gives  examples  of  problem 
runs  vith  corresponding  POLICY  subroutines. 

So  far  this  subject  has  been  discussed  in  the  context  of  choices 
to  be  made  within  the  POLICY  routine.  Numerous  other  options  may  be 
seJected  by  reading  in  flags  or  constants  as  part  of  initial-condition 
data.  Ecwever,  it  is  sometimes  desirable  to  override  these  Initial 
instructions  in  POLICY  if  during  a  problem  run  a  situation  arises  that 
requires,  perhaps,  a  change  in  printout  frequency  or  integration  step 
sice. 
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III.  INPUT-OUTPUT  GEOMETRY 


In  order  Co  construct  a  problem  run,  it  is  necessary  to  define  the 
problem  by  reading  In  Initial-condition  data  and  establishing  a  flight- 
control  program  (P0LICY  subroutine) .  Before  a  detailed  description  of 
this  process  Is  possible.  It  Is  essential  to  have  a  basic  acquaintance 
with  certain  fundamentals  of  the  coordinate  system  and  intercept  geometry 

ABSOLUTE  AND  RELATIVE  POSITION  GEOMETRY 

Figure  3  shows  two  position  vectors,  r^^  and  r^,  in  a  space  referred 
to  a  three-dimensional  x,  y,  2  coordinate  system  of  arbitrary  origin. 

The  position  of  point  1  may  be  defined  either  in  Cartesian  coordinates, 
^1’  ^1*  ^1*  spherical  coordinates,  r^,  0^^,  Because  there  are 

certain  advantages  to  each  form  of  coordinate  system,  TACTICS  can  con¬ 
vert  one  to  the  other,  so  that  all  vectors  are  expressed  in  six-element 
form,  as  in  the  following: 

r^  -  (x,  y,  2,  r,  0,  9) 

Initial-condition  position  data  may  be  read  in  either  Cartesian 
or  spherical  coordinates  with  respect  to  the  reference  frame.  Generally, 
the  former  is  the  most  convenient  for  this  purpose.  For  example,  it  is 
easy  to  visualize  that  one  vehicle  is  initially  at  an  altitude  of  ft 
and  arbitrarily  placed  at  Xj^  »  0,  y^  =  0,  while  perhaps  another  vehicle 
is  initially  at  an  altitude  of  ft  and  displaced  horizontally  from 
the  first  by  X2,  y2  ft. 

Once  initial  conditions  have  been  established  and  the  vehicles 
located  in  a  reference  coordinate  frame,  it  is  generally  true  that 
the  absolute  displacements  of  the  vehicles  with  respect  to  this  frame 
are  no  longer  of  primary  interest.  Intercept  problems  are  mainly 
functions  of  relative  geometry,  i.e.,  relative  position  and  velocity 
or  components  thereof. 

In  Fig.  3,  the  vector  rj^2  J^epresents  the  range  and  direction  of 
point  2  relative  to  point  1.  For  three  vehicles,  there  is  the  vector 
r^,  not  shown.  More  generally,  then,  the  relative  position  vectors  are 
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These  quantities,  as  veil  as  all  other  quantities  pertaining  to  rela¬ 
tive  geoaetry  (usually  the  entire  intercept  probles),  are  calculated 
and  used  extensively  within  the  progras  and  fom  an  Isportant  portion 
of  the  printout.  As  an  exaaple,  guided-sassile  control  laws  and  cri¬ 
teria  for  missile  launch  are  usually  dependent  on  relative-position 
and  angular-rate  data.  Because  relative-position  information  is 
visualized  and  interpreted  more  conveniently  in  spherical  coordinates, 
it  is  printed  out  in  this  fom. 

In  siunmary: 

o  Initial-condition  position  data  may  be  specified  in  either 
Cartesian  (y.,  y,  z)  or  spherical  (r,  0,  9)  coordinates  with 
respect  to  the  reference  frame.  Units  of  measurement  are 
feet  and  degrees. 

o  All  relative  geometry  (e.g.,  computed  within  the 

program  and  is  printed  in  spherical  coordinates. 


ABSOLUTE  AND  RELATIVE  VELOCITY  GEOMETRY 

Input  velocity  data  and  position  data  are  treated  almost  identi¬ 
cally  to  one  another.  Figure  4(a)  shows  a  velocity  vector  associated 
with  the  position  vector  r^^.  For  initial  conditions,  may  be  entered 

in  either  Cartesian  (V  ,  V  ,  V  )  or  spherical  (V,  6  y)  form  (see 

X  y  z  V 

Fig.  4(b)).  Conversion  from  one  form  to  the  other  cakes  place  within 
the  program.  The  unit  of  measurement  is  feet  per  second. 

As  in  the  preceding  discussion  of  position  data,  the  relative 
velocity  vectors  are  calculated  within  the  program  and  occupy  signifi¬ 
cant  portions  of  the  printout.  Similarly,  the  relative  velocities  are 
given  by 
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1  -  1,  2,  3 
j  “  1,  2,  3 
i  *  j 


and  so  on  for  etc.  As  with  position  vectors,  coordinate 

transformations  are  performed  so  that  these  vectors  are  also  expressed 

slx-dlmenslonally  In  terms  of  x,  y,  z,  V,  6^,  and  y.  It  is  frequently 

necessary  in  the  computations  to  call  upon  both  absolute  (r.  and  r.)  and 
•  XX 

relative  (r^^  and  r^^)  quantities.  The  angular  rates  of  rotation  of 
the  range  vectors,  both  absolute  and  relative,  are  of  basic  importance 
to  most  guidance  or  trajectory  problems.  Expressed  in  vector  notation. 


1  ri  1  i 


where  r^  =  scalar  time  rate  of  change  of  the  vector  r^ 


1^^  »  a  unit  vector  along 


u)^  ®  angular-rate  vector  orthogonal  to  both  r^  and 


Similarly,  for  the  relative  velocity  vectors 


hi  ■  'ijlry  +  “tj  'ij 


It  is  important  to  note  that  for  approaching  vehicles  the  relative 

f 

range  rate  r^^j  is  negative;  for  vehicles  which  are  separating,  it  Is 
positive.  Also,  the  point  of  closest  approach  or  minimum  miss  distance 
between  two  vehicles  occurs  when  the  absolute  value  of  this  relative 
range  rate  la  zero. 

Changing  to  spherical  notation  for  convenience,  let  us  consider 

s  • 

the  vector  velocities  r^  and  r^  or  and  ,  respectively.  The 
following  relationship  for  the  angular-rate  vector  representing  the 
rotation  of  the  vector  r^^  Is  applicable: 


“iJ  "  ’'ij 
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Note  that  with  these  basic  relationships  it  is  possible  to  re¬ 
solve  vector  velocities  into  two  useful  components,  one  along  r^^  and 
the  other  transverse  to  Also,  the  angular  rotation  of  the  line 

of  sight  (LOS)  between  r^  and  r^  is  known  both  in  direction  and  mag¬ 
nitude.  This  latter  quantity  is  of  prime  importance  in  most  terminal 
(or  command)  missile  or  aircraft  guidance  applications.  In  actual 
practice,  it  is  usually  a  quantity  determined  by  rate  gyro  measurements 
(in  missiles)  or  by  processing  9,  (p  angle  data  measured  by  a  ground 
radar. 

So  far,  relationships  have  been  given  for  resolving  the  velocity 
vectors  into  components  parallel  and  perpendicular  to  the  r .  and  r 

i  Ij 

values  and  for  calculating  the  angular-rate  vectors.  In  subsequent 
discussions  of  guidance  laws  and  acceleration,  two  useful  vectors  as¬ 
sociated  with  velocity  will  become  Important.  Imagine  a  plane  normal 
to  a  velocity  vector  "r  having  spherical  coordinates  V,  0^,  y.  It  is 
convenient  and  useful  to  define  a  unit  vector  1^  common  to  this  plane 
and  the  horizontal  x-y  plane,  and  a  unit  vector  T_  common  to  this  plane 

•  L* 

and  the  vertical  plane.  Thus,  i*  or  V  and  the  unit  vector  1^  will  form 
a  right-hand  orthogonal  system  with  1^  and  Tj^,  as  shown  in  Fig.  5.  The 
following  relationships  are  applicable  for  determining  the  components: 


l^jj  =  -  sin  0„  1^^  =  -  sin  Y  cos  8„  1,,^  =  cos  y  cos  0 


V  Dx 


V  ^Vx 


^Dz  =  ^ 


^Vz  =  ^ 


ELEVATION,  AZIMUTH.  AND  BEARING-ANGLE  GEOMETRY 

An  Important  consideration  in  most  Intercept  problems  is  the 
direction  or  orientation  of  the  LOS  between  one  vehicle  and  another. 

A  simple  example  of  the  calculation  is  a  description  of  relative  tar¬ 
get  position  by  a  pilot  in  a  cockpit  in  terms  of  numbers  on  a  clock 
face  (12  o’clock,  straight  ahead;  3  o’clock,  directly  to  the  right; 
etc.)  with  elevation  designated  as  high  or  low.  Obviously,  this 


Fig,  5  Unit  vectors  normal  to  velocity  vector 
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deseriptisn  wculd  zcouire  alteratios  if  the  aircraft  were  suddenly  to 
anlce  a  sobstaacial  <£iscge  £a  attitude  or  orientation  by  pitching,  roll- 
yjwing.  Siaiilarly,  aACncS  calculates  and  prints  out  aziauth 
Cclock— face  poaitioa}  and  elevation ( (hi gh~low)  orientation  of  the  lOS 
ftou  ea^  vehicle  relative  to  the  other  tvo.  It  also  prints  out  a  bear¬ 
ing  angle  defined  here  as  the  angle  Between  the  ICa  and  the  longitudi- 
axis  of  the  vriiicle.  Azixuth  angle  is  oeasured  in  the  pl(uxe  formed 
ihe  ri|^t  wing  and  the  longitudinal  axis  as  a  pilot  would  view  the 
relative  geoaetzy.  Slwilarly,  elevation  3ng}.e  Is  oeasured  in  the  plane 
foraed  by  the  loogitudinai  axis  and  the  lir.e  chrou^  the  top  of  the 
cockpit.  Azinoth  angle  ranges  frea  0  to  1180  deg,  with  positive  to  the 
ri^t.  Slevatloc  angle  ranges  troa  0  to  ±90  deg,  with  positive  upward, 
as  shofra  ic  Fig.  6.  Hote  chat  the  bearing  angle  is  defined  as  the  to¬ 
tal  angle  between  the  lOS  and  longitudinal  axis  and  hence  is  independent 
of  roll  angle  (see  Appendix  A  for  definition);  this  is  a  very  useful 
feature^  since  it  is  confusing  to  visualize  azlouch  and  elevation  angles 
when  a  vehicle  is  saneuverlng. 

For  those  prealeas  which  involve  sensors — e.g.,  infrared,  radar, 
or  optical — If.  is  likely  that  constraints  will  be  Imposed  on  naxiamm 
values  of  elevation  and  on  aziauth  look  angles  due  to  mechanical,  elec¬ 
trical,  cr  optical  Units.  Provisions  have  accordingly  been  sade  tc 
li^ut  any  such  constreints  as  Initial-condition  data. 

Mathsaatical  details  and  derivations  pertaining  to  aircraft  atti¬ 
tude  or  orientation  angles  and  elevation,  azinutb,  and  bearing-angle 
geonetry  are  contained  in  Appendix  A. 

In  auasary: 

o  Associated  with  each  vehicle  are  two  LOS  oriented  toward  Che 
other  two  vehicles. 

o  The  vehicle's  coordinate  system  may  be  regarded  as  a  longitudi¬ 
nal  axis  through  the  airframe,  a  line  through  the  right  wing, 
and  a  line  through  the  top  of  the  cockpit. 

o  The  orientation  of  the  LOS  with  respect  to  this  coordinate 
system  is  resolved  in  terms  of  azimuth,  elevation,  and  bear¬ 
ing  angles. 
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IV.  DEFIN1!«;  THE  FORCES  AND  ACCELERATIONS 

COKIROL  LAWS 

TACTICS  integrates  the  equations  of.  motion  defined  by  the  three 
components  of  net  acceleration  associated  with  each  of  the  vehicles. 

A  resultant  vector  force  F  will  define  an  acceleration  a.  The  forces 
of  primary  interest  which  will  sum  to  this  resultant  force  F  may  be 
categorized  as  (1)  gravitational,  (2)  aerodynamic  lift  and  drag,  and 
(3)  propulsive.  If  we  assume  that  these  forces  are  defined  in  magni¬ 
tude  and  direction,  it  is  a  straightforward  procedure  to  add  them, 
resolve  them  into  components,  and  determine  the  net  acceleration.  Fact 
control-law  subroutine  may  be  considered  as  a  modular  unit  where  this 
process  or  its  equivalent  is  performed,  the  output  being  the  three  com¬ 
ponents  of  net  acceleration  applicable  to  a  particular  vehicle,  as 
shown  schematically  in  the  block  diagram  below. 


Aerodynamic  force 

Control 

Law 

_2i - ► 

Gravitational  force 

Propulsive  force  _ 

Vehicle  (i) 

Mass 


The  important  point  to  note  in  this  diagram  is  the  correspondence 
of  the  net  sum  of  forces  to  the  net  acceleration.  Indeed,  there  are 
many  control  laws  in  which  the  net  acceleration  is  by  definition  the 
starting  point,  and  all  control- law  computations  are  primarily  con¬ 
cerned  with  finding  the  correspondence  between  the  forces  that  would 
be  necessary  to  create  such  an  acceleration. 

To  Illustrate,  consider  a  hypothetical  control  law  in  which  all 
components  of  net  acceleration  are  by  definition  sero.  Assume  that 
propulsive  and  gravitational  forces  are  determined.  The  function  of 


Che  control  lav  in  this  sl^>le  case  Is  to  deCeralne  Che  aerodynamic 
lift  and  drag  forces  necessary  to  guarantee  the  postulated  output  con¬ 
dition,  i.e.,  zero  acceleration. 


LATERAL  ACCELERATION  AND  ITS  COMPONEOTS 

Many  reference  texts  differ  on  the  definition  and  usage  of  the 
terra  “lateral  acceleration."  Throughout  this  Memorandum,  it  is  defined 
as  a  vector  quantity  in  a  plane  normal  to  the  velocity  vector  or  flight 
path  of  the  vehicle.  Since  the  velocity  vector  V  may  be  oriented  in 
any  direction  in  the  general  case,  the  above  definition  does  not  con¬ 
strain  the  lateral  acceleration  vector  to  any  particular  direction  or 
to  any  plane  other  than  Che  one  normal  Co  V. 

Referring  to  Fig.  5,  1  and  1  are  also  by  definition  in  tit:  same 

A  \J 

plane  normal  to  V.  Accordingly,  it  is  convenient  to  resolve  the  lat¬ 
eral  acceleration  a  into  two  components,  a^^  and  a^. 


a 


V 


=  a  ♦ 


(7) 


By  definition,  the  component  a^^  now  represents  a  turning  acceleration 
in  the  horizontal  plane  and  a^  represents  a  climbing  or  diving  accel¬ 
eration  in  the  vertical  plane  (caution;  a^  is  not  necessarily  in  the 
vertical  direction).  At  this  point,  it  is  necessary  to  distinguish 
between  the  specified  or  commanded  values  of  a  and  its  components  (a^ 
and  Sy)  and  the  modified  or  output  values.  In  simulating  guidance  of 
aircraft  or  missiles,  real-world  considerations  often  necessitate  the 
modification  of  commanded  values,  because  of  constraints  such  as  struc¬ 
tural  or  aerodynamic  limitations,  time  lag,  and  noise.  This  subject 
will  be  discussed  further  in  Section  V. 

For  the  moment,  let  us  assume  that  commanded  values  are  constrained 
or  modified  so  that  output  values  result.  With  subscripts  used  to 
denote  the  difference,  a  commanded  value  of  lateral  acceleration  is 
designated  as  a^  with  components  and  a^^,  and  a  modified  or  output 


value  as  a^  with  components  a^j^  and  a^^.  The  third  component,  V,  is 
assumed  to  he  a  function  only  of  propulsion  and  aerodynamic  drag  forces 
with  no  constraints  or  modification.  The  total  resultant  acceleration 
V,  describing  the  motion  or  trajectory  of  a  vehicle  represented  as  a 
point  mass,  is  then 


V  =  a  ,  1.  +  a  „  1,  ’I-  V  1„ 
oh  A  oV  D  V 


where  1^  is  the  unit  vector  along  V  and  V  is  the  time  rate  of  change 
in  the  magnitude  of  V.  Expressed  in  inertial  Cartesian  (x,  y,  z)  form, 
which  is  more  convenient  for  numerical  integration  of  the  equations  of 
motion,  the  components  are 


X  =  V  •  i  =  a  .  1.  +  a  ,,  1^  +  V  1„ 

oh  Ax  oV  Dx  Vx 


y  =  V  *  j  =  a  ^  1.  +  a  „  U  +  V  1„ 

oh  Ay  oV  Dy  Vy 


z  =  V  •  k  =  a  ,,  1^  +  V  1„ 
oV  Dz  Vz 


where  i,  j,  and  k  are  the  unit  vectors  along  the  x,  y,  z-coordinate 
frame  axes,  respectively,  and  1^^,  1^^,  etc,  are  the  components  of 
the  unit  vectors  given  by  Eq.  (6). 

To  represent  the  flight  of  three  vehicles  in  motion  simultaneously 
18  differential  equations  must  be  integrated,  9  involving  the  accelera- 
tions  V(l)  and  9  the  velocities  V(l).  TACTICS  expresses  these  differ¬ 
ential  equations  in  two  different  forms:  One  uses  the  flat-earth 
representation  (less  complex  and  faster)  and  the  other  uses  a  round 
rotating  or  nonrotating  earth  (essential  for  space  applications,  long 
ranges,  or  high  speeds).  The  expressions  in  Eq.  (9)  are  the  basic 
acceleration  equations  which  are  integrated  for  the  flat-earth  repre¬ 
sentation.  Further  details  on  numerical  integration  methods,  round- 
earth  form  of  the  equations,  coordinate  transformations,  and  deriva¬ 
tions  are  included  in  Appendix  H.  Next,  the  correspondence  between 
net  lateral  acceleration  and  the  pro;vitatlonal ,  propulsive,  and  aero¬ 
dynamic  forces  which  cause  this  acceleration  will  be  considered.  At 
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thls  point  it  is  important  to  stress  that  a  vehicle  maneuver  may  be 
defined  in  either  of  two  ways,  depending  on  the  particular  problem. 

A  complete  force-acceleration  correspondence  must  be  established  by 
specifying  either  (1)  the  lateral  acceleration  and  sufficient  other 
information  about  the  forces,  or  (2)  the  forces  and  sufficient  other 
Information  about  the  lateral  acceleration. 

FORCES  AND  FORCE  RESOLUTION 

The  various  forces  assumed  to  be  acting  on  the  center  of  gravity 
(c.g.)  of  each  vehicle  ate  shown  in  Fig.  7.  In  vector  form,  the  summa¬ 
tion  of  these  forces  is 


mV-mg  +  L  +  T  +  D  (10) 

where  V  *  vehicle  acceleration 

m  »  mass  of  the  vehicle  (m  •  W/g) 

g  =  acceleration  due  to  gravity  (unit  vector  in  the  -z  direction) 

L  =•  lift  force  (defined  as  normal  to  the  velocity  V) 

T  =  resultant  thrust  force 

D  =  dreg  force  (unit  vector  in  the  -V  direction) 

This  is  a  general  vector  expression  where  the  positive  or  negative  signs 
are  accounted  for  by  the  orientation  of  the  unit  vectors  associated  with 
each  term.  For  notatlonal  convenience,  the  1  subscript  which  refers  to 
a  particular  vehicle  has  been  omitted.  The  above  equation  may  be  re¬ 
solved  and  separated  into  two  equations,  one  representing  the  lateral 
acceleration  normal  to  V  and  the  other  representing  the  V  acceleration 
in  the  direction  of  V.  Before  doing  this,  certain  important  definitions 
or  assumptions  are  necessary: 

o  The  thrust  force  T  acting  through  the  c.g.  is  also  coincident 
with  the  longitudinal  body  axis  of  the  vehicle, 
o  The  angle  of  attack  a  is  taken  to  be  the  angle  between  T  and 
the  tangent  to  the  flight  path  V. 
o  All  maneuvers  consist  of  "coordinated  turns,"  defined  by  the 
condition  that  a  plane  passed  through  the  longitudinal  body 
axis  and  including  the  yaw  axis  must  also  include  the  vector  V. 
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Ihe  sicuation  for  a  coordicated  tuzs  in  level  fll^c  is  sboim  ia  Fig.  7. 
Hie  vectors  L,  T,  and  V  are  co-plaaar  and  T  is  coincidesc  vitb  the  lon¬ 
gitudinal  axis. 

For  the  particular  case  of  a  level-flight  turn,  the  bank  angle 
is  shown  to  be  the  angle  between  L  and  the  vertical,  but  in  the  general 
case  ir  is  the  angle  betveen  L  and  the  reference  vector  1^  shown  in 
Fig.  3.  C!!ote  that  tiiis  is  a  vied  axis  or  velocity  reference,  since 
1^  is  normal  to  V.)  Kith  these  definitions  and  assunptions  in  ssind, 

£q.  (10)  say  now  be  re^irolved  into  two  scalar  equations  by  suaalsg  forces 
and  accelerations  a?,ong  V  (direction  1,^  and  in  a  plane  nomal  to  V, 
using  the  sysbol  to  denote  t!ie  direction  of  this  resulting  noiaal  or 
lateral  acceleration  (see  page  22  for  definition  of  lateral  acceleration). 


V  *  ^  (T  cos  a  -  D  -  W  sin  y) 


a 

o 


[ 


J  (L  +  T  sin  c) 


(11) 

(12) 


where  a^  is  the  absolute  magnitude  of  net  lateral  acceleration  a^  and 
a  is  the  angle  of  attack,  as  defined  above. 

Two  exas^les  will  illustrate  the  basic  relationships  in  order  to 
avoid  possible  aisunderstanding  of  definitions  or  terminology. 


Case  1;  Determination  of  Forces  Reeuired 
For  a  Specified  Lateral  Acceleration 

Assume  a  simulation  of  a  climbing  turn  arbitrarily  defined  by  a 
lateral  acceleration  component  a^^^  in  the  horizontal  plane  and  a  posi¬ 
tive  a^y  component  in  the  vertical  plane.  If  the  maneuver  is  initiated 
from  level  flight  (i.e.,  V  in  the  horizontal  plane),  the  a^^  cooqionent 
and  Ip  will  initially  be  oriented  toward  the  vertical  -  direction  but  not 
thereafter.  Figure  8  shows  the  vectors  and  their  scalar  magnitudes. 

The  total  net  lateral  acceleration  will  be 

n~7~2 

a  -/a  ,  4-  a  „ 
o  •,/  oh  oV 

V 


(13) 
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The  dorms.!  force  F^,  defined  as  the  sum  of  aerodynamic  lift  and  propul¬ 
sive  forces  normal  to  V,  is 

F  »  L  +  T  sin  a  *  ~y(.g  cos  y  +  a  (14) 

n  g  V  '  oV  oh 

The  direction  of  the  vector  L  is  given  by 

j  _  »[%h  ^  ^  ^  \] 

'■  * 

Note  that  the  gravitational  force  component  normal  to  V  is  asso¬ 
ciated  with  the  g  cos  y  term. 

Assuming  a  small,  the  initial  (y  »•  0)  lift  force  L  required  for 
the  maneuver  is 


The  direction  of  a  is 
o 

(a  ,  1  +  a  ,,  1_) 

1^  =  — 

“o 

and  the  bank  angle  is 

h  -  (r+THST  ■  i) 

The  minus  sign  is  arbitrarily  assigned  to  make  negative  for  a  left 

D 

turn,  i.e.,  a^j^  is  positive  for  6^  increasing  (shown  in  Fig.  4). 
Accordingly,  in  terms  of  bank  angle  the  components  of  a  become 


®oh  *  ■  W  '^'b 


The 
the  lift 


force  F^  is  defined  as  the  force  acting  in  the  direction  of 
vector  L.  Expressed  in  g's,  it  is 


(21) 


^  ^  L  -f  T  sin  g 
W  W 


From  Fig.  5, 


cos  i;»g  = 


Sov  +  g  cos  Y 


Hence, 


?  a  „  +  g  cos  Y 
n  _  oV  _ _ 

W  g  cos  ij»_ 


For  a  horizontal  turn  in  level  flight,  as  shown  in  Fig.  7,  this 
reduces  to  the  simple  expression 

F  . 

W  cos  \l>„ 

Case  2;  Determination  of  Acceleration  Components 
Resulting  From  Specified  Forces 

Assume  that  aerodynamic  and  propulsive  forces  may  be  determined 
utilizing  some  prescribed  control  law  and  that  it  is  necessary  to  deter¬ 
mine  the  acceleration  vector  components  a^j^,  and  V  from  Eq.  (10), 


which  is  repeated  below; 

V  =  g  +  J  (L  +  T  +  D)  (2< 

Ignoring  for  the  moment  the  effects  of  time  lag,  aerodynamics,  or 
structural  constraints  and  referring  to  Eqs.  (6),  the  acceleration 
components  are 

.  j_  _ 

V  -  V  ♦  1„ 


“oh  ■  ’  •  1a 


“oV  ■  •  Id 


The  orientation  of  the  lateral  acce.l  eratlon  vector  a  is  given 


by  Eq.  (16) . 


« 
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AERODYNAMIC  FORCES 

Aerodynaialc  computations  may  involve  a  large  number  cf  input  param¬ 
eters  corresponding  to  the  flight  characteristics  of  a  particular  air¬ 
frame.  On  the  other  hand,  there  are  trajectory  simulation  problems 
where  aerodynamic  effects  are  not  pertinent  and  no  aerodynamic  compu- 
tatloi .  are  necessary  for  a  particular  vehicle.  An  example  is  a  point- 
mass  target  moving  along  a  path  where  the  motion  is  arbitrarily  postu¬ 
lated,  say,  a  straight-line  path  at  constant  speed. 

In  the  general  case,  expressions  for  force  magnitudes  (lift  L  and 
drag  D)  are 


L  =  C  Aq 

I)  =  CpAq  (26) 

q  =  §pV^ 


where  =  lift  coefficient 
Cp  =  drag  coefficient 

2 

A  »  aerodynamic  reference  area  of  the  vehicle  (ft  ) 

2 

q  =■  dynamic  pressure  (Ib/ft  ) 

p  =  air  density  (slug/ft^)  (see  Appendix  I  for  a  definition  of 
model  atraosphere)_^^^ 

In  simulating  the  flights  of  s^ctf-lQ.^ircraft  designs,  provisions 
have  been  made  for  Incorporating  tables  to  describe-.t^  interactions 
among  lift  coefficient  drag  coefficient  C^,  angle  of  art'atk--a.,._ 
maximum  lift  coefficient  and  Mach  number.  Interpolation  rou¬ 

tines  are  employed  to  provide  the  equivalents  of  functional  relation¬ 
ships.  Alternatively,  provisions  have  been  made  for  using  the  follow¬ 
ing  analytic  expressions: 

C,  =  -3 —  (a  -  a  ) 

L  da  o 


dC 


D 


d(cj) 


s 

\ 


(27) 
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taSaere  «5C  /ds  *  slope  o£  t&e  C  »  F(c:}  ctsrsre 
£» 

a  «  zero-life  as^c  oi  accadk.  (dez) 
o 

»  zes©-ilfc  dsas  c5«££lcl«ct 
a 

cC^dff^  «  a  coef€l£iaai  osed  v£c&  a  parabolic  fesecioa  fer  drag 

cecfficiesr,  i.e.,  ^  treaced  as  a  cos- 

acasr  or  a  foacdsa  of  Xacs  raanber} 

Refra^  Is  as.  hew  eo  prosids  ao^  cse  rabulared  calces  zeed  er'  hsw  co 
s^ecr  Ehe  varices  cp^ioss  available  ars  provided  is  Secrioss  IS 
chrzcgb  XIII  aad  Appeadijc  D  of  ehir  Xeanargradiifla. 

PgO?EISiy£  HaCFi; 

As  prevloosl^  siesricseds  propoisive  forces  zre  sssamed  co  ace  cs 
che  c.g.  of  Che  v^icle  asd  co  be  coiccideac  wich  the  vehicle's  losgi- 
rodinai  axis.  As  ia  asrCcdToaicic  cco^cacioas,  cables  aav  be  izu»>rpo- 
raced  co  describe  che  i&cerrelacicKsships  beewaea  siilicarv  afeer- 
btiraer  chrasc,  foel  flow,  as^  jfach  so3d>ec.  Ihere  are  also  previsioas 
for  reading  in  cc-nscaaC  val^zes  for  chrosc  as  ispuc  daca.  A  varying 
chrese  condicios  sny  be  siscilaced  by  che  siaple  ezpedienC  of  xnlcip'^ 
ing  each  cl!n:at  value  by  a  "chroccle''  paraaecer  (norsally  set  to  1.0 
unless  ochexvise  specified  m  P0LIC1).  Addlclonal  flexibilicy  is  also 
available  by  incerporacisg  chrusc  values  or  fusctiocal  relaciccshlpa 
ic  maneuver  scbroucines  specialized  to  cbe  parcicular  problem  or  to 
vehicle  characreristirs. 

GRAVITATIONAL  FOSCES 

TACTICS  is  auCosacically  zcC  for  Che  aia^lesc  assuzsytion  of  a 

2 

gravitational  force — a  fiat  earth  causing  1  g  or  52.174  ft/^ec  acting 
doKnwj»rJ  in  the  negative  z-directicn.  However,  options  are  provided 
for  considering  the  core  conplicated  cases  of  a  round  earth,  rotating 
or  nonrocating,  so  that  centrifugal,  Coriolis,  and  Inverse  square  law 
ef'ccts  nay  be  included  If  pertinent  to  the  prebien,  e.g. ,  in  space 
applications  and  hypersonic  flight.  The  basic  expressions  for  the 
tlat-ea’'l.h  representation  are  given  by  Eq.  (9),  and  the  gravitational 


t 


force  Cieras  «ppe«rs  in  a  racber  acralghtforward  way  (which  is  £0  be 
expected,  since  %  is  treated  as  a  ccostant  vector,  i.e.,  in  nagnlCude. 
and  direction)  in  £^s.  (11)  through  (16),  All  expressions  and  deriva- 
ticns  relating  to  the  round-earth  representation  are  discussed  in 
Appendix  3  because  cf  the  large  mmber  of  details  involving  coordinate 
transfcriBatioss  asd  reference  vectors.  Ibe  gravitational  acceleration 
in  this  case  Is  taken  to  be  rotating  and  varying  in  aagnltude,  as 


where  »  gravitational  force 

V  -  1-407645  10^°  ft^/sec^ 

R  s  radius  vector  from  the  geocenter  to  the  vehicle 

R  Ir j 
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V.  AERODYNAMIC,  STRUCIURAL.  AND  TIME-IAG  CCKiSTRAiyiS 

AESCDYMAMIC  ASP  STRUCIlJaAL  COXSTRAISIS 

As  laeationed  la  Section  1V»  in  sianlating  vehicle  flight  perfor- 

aance  it  is  usually  necessary  to  ispose  constraints  or  lixits  cn  the 

forces  applied  to  the  vehicle.  Two  piinary  liaitations  are  (1)  aerch- 

dtp^txsia,  l.e.,  a  llxitation  on  lift  coefficient  to  soae  specified 

value  representing  a  boundary  cn  flight  stability,  ?cd  (2)  stiucturalt 

i.e.,  a  load  lloit  lapsed  by  possible  daacage  to  the  airfrase  or  cost- 

ponents  (or  possibly  to  a  huaan  being).  The  aerodynasic  constraint 

iopcses  the  following  condition  on  the  aiagnitude  of  the  noraal  force 

F  (see  Bq.  (14)): 
n 


F 

n 


2  C,  Aq  + 
iMax  ^ 


T  sin  o 


nax 


(29) 


where  C,  is  a  naxisma  value  for  G,  and  a  is  a.  saxisun  value  for 
Laax  I,  3iax 

a.  Since  C,  and  c  are  functionally  related,  either  one  my  be  specified 
and  the  other  calculated.  Similarly,  tee  structural  constraint  is 

F  (30) 

n  Snax 

where  a^  is  the  structural  for  hunan)  acceleration  lizsit.  It  sh-'^uld 
Smx 

be  noted  that  by  definition  the  force  F  is  noraal  to  the  velocity  vec- 

o 

tor  V.  To  be  precise,  the  structural  linit  a^^^  should  be  considered 
as  noraal  to  the  vehicle’s  iougitudlnai  axis,  which  is  separated  frea 
V  by  the  angle  u.  In  the  TACTICS  prograa  and  In  the  derivations  which 
follow,  this  difference  is  ignored  and  is  assuzaed  to  be  a  lateral 

acceleration  linit,  by  previous  definition  also  noraal  tc  V.  Designat¬ 
ing  the  specified  or  cosputed  value  of  the  nonaal  force  as  the 

aercdynaaic  limit  as  F  ,  and  the  structural  limit  as  F  ,  TACTICS 
•'  na  ns 

takes  he  applied  normal  force  to  be 


F  -  min(F  F  ,  F  ) 
n  nC  na  ns 


(31) 


f 
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wiere  ain  is  the  ■laiaua  aagnltade  of  Che  three  values.  Khen  a  con¬ 
straint  value  or  is  taken,  the  initial  assunptlon  is  that 
the  specified  or  coaputcd  direction  of  the  commanded  acceleration,  i.e., 
Che  unit  vector  Ij^,  reaains  constant  but  that  the  Magnitudes  of  the 
forces  and  corresponding  accelerations  should  be  cot^tible  with  the 
constraining  value.  1£  a  is  designated  as 


■ax 


[§Jain  (F  ,  ?  ) 

iWI  '  na*  ns' 


(32) 


the  problem  is  to  deceradne  the  corresponding  aagnitude  of  the  lateral 
acceleration  eagnicude  a^.  The  applicable  expression  is 


nax 


•  o  1 


g  cos  Y 


(33) 


from  wiiich  the  following  equation  is  derived: 

3^  +  2  g  cos  Y  (1^  •  Ijj)  a^  +  (g  cos  Y>^  -  *  0  (34) 


If  the  dot  product  is  denoted  as 


cos  6  «  (1^  •  Ijj) 


(35) 


the  sol^jtion  for  a  is 
o 


2  2  2  2 

a  ■  -g  cos  Y  CO*  5  t  la  -g  ccs  y  *in  6 

o  *  -1/  nax  * 


(36) 


To  solve  the  quadratic  equation,  the  following  conditions  sust  be 
satisfied: 


o  The  value  of  sust  by  definition  be  positive,  since  It 

represents  the  aagnitude  of  a  vector.  (For  multiple  positive 
solutions  the  largest  aagnitude  is  taken.) 
o  The  discriiclnant.  should  not  be  negative. 

Difficulties  in  maintaining  the  ''Irectlon  of  a  (i.e.,  1.)  arise  when 

o  1 

a  becomes  small  (e.g.,  a  wingless  vehicle  with  small  C,  ).  These 
teax  V  e>  »  e-  Lmax 

special  conditions  are  handled  as  follows: 


The  block  labeled  “Llalt"  represents  the  possible  laq>osition  of  aero¬ 
dynamic  or  structural  constraints  as  discussed  previously.  Since  cer¬ 
tain  arbitrary  assuaptlons  are  involved  in  each  of  these  processes,  the 
purpose  of  the  diagram  Is  to  emphasize  the  aodclar  building-block  form; 
alternate  subroutines  adapted  to  a  particular  problem  nay  be  substituted. 

The  time-lag  transfer  function  describes  the  performance  of  the 
vehicle  hardware  mechanization  in  terms  of  control-system  response, 
gyro  prediction,  tracking  system,  filter  circuits,  etc.  Here  the  main 
generalization  of  the  simulation  occurs,  for  the  question  arises  as  to 
the  pertinence  of  this  aspect  of  vehicle  performance  to  the  problem 
solution.  Unless  ^  specific  hardware  design  is  being  investigated. 


VI.  CaiDAXCE  AND  COItIBOL  LAMS 


Generally,  the  call  for  a  specified  control  latf  is  conditional 

ispoD  the  fclfillaent  of  one  or  nore  ge(»etric,  klneeatic,  or  tine  con- 

ditlocs.  Ihe  tern  "control  lau"  as  used  here  asms  defining  a  cocnanded 

■ri  lateral  acceleration  vector  aL  associated  vith  the  vector  velocity  V. 

f  c 

For  those  lavs  involving  aircraft  or  aissile  guidance,  such  as  lead  col¬ 
lision  or  proportional  navigation,  the  calculation  of  the  vector  quac- 
’  tity  is  dependent  upon  the  kioeaatic  state  of  the  systea  (position, 

\  1  velocity,  and  perhaqts  acceleration) ;  these  nay  be  classified  as  closed- 

loop  guidance  levs.  On  the  other  hand,  there  are  open-loop  control  lavs, 
c.g. ,  turn,  dive,  cliii>,  etc.,  in  which  the  klnexatic  state  of  the  aya- 
'  tea  is  not  inliclt  within  the  law  itself.  This  difference  should  be- 

coae  aore  apparent  as  these  lews  are  suaaarized  later. 

TWenty-four  different  control  laws  are  described  in  this  Meaorandua, 
aany  of  which  require  lengthy  explanations  and  derivations.  For  con¬ 
venience,  only  a  fctw  illtistrative  exaaples  are  given  below;  a  detailed 
'  listing  of  all  control  lavs  and  their  derivations  is  given  in  Appendix  C. 

A  brief  reviev  of  zeros  and  of  the  basic  applicable  acceleration 
equations  follows: 

1.  The  coananded  lateral  acceleration  vector  a^  is  by  definition 
in  a  plane  noraal  to  V  (sec  page  21)  and  is  subject  to  possible  aero¬ 
dynamic,  structural,  and/or  time-lag  constraints  resulting  in  the  for¬ 
mulation  of  the  vector  a  (see  Section  V). 

_  o 

2.  The  vector  a  is  resolved  inlo  two  components  by  the  unit 

_  _  o 

vectors  1^  and  1^  so  that 


a  V.  “  ® 
oh  o  A 

a  ■  a  •  1_ 


5 


3 4  If  sU  forces — grsvitatioxial,  propulsive,  and  aercdynaaic — 

» 

are  specified,  the  total  nvc  acceleration  V,  consisting  of  the  cob- 
ponents  a^^,  a^^,  and  V,  say  be  deteraiced.  On  the  other  hand,  if  a^y 
and  a  .  (gravitational  and  propulsive  forces)  arc  specified,  the  aero- 

Oil 

dynaalc  forces  aay  be  determined  to  colculate  the  corresponding  coa- 
penent  along  V  (see  Eq.  (25)). 

4.  The  basic  acceleration  equations  (flat-earth)  are  given  by 
Eq.  (9)  in  Cartesian  (x,  y,  z,)  coordinates,  frea  which  the  velocities 
and  positions  are  detereined  for  each  vehicle  after  nuaerical  integra¬ 
tion. 

The  guidance  lavs  given  in  the  illustrations  below  are  alacst  all 
of  the  type  in  which  the  vector  a^  and  propulsive  forces  are  specified 
and  the  corresponding  aerodycasic  forces  are  to  be  detemined. 

OPEK-LOOP  CONTROL  LAWS 

Straight  Flight 

The  coBsanded  lateral  acceleration  a^  Is  zero.  The  vehicle  will 

o 

fly  a  straight-line  path  (but  not  necessarily  "straight  and  level**). 
However,  an  acceleration  or  deceleration  along  this  path  may  occur  due 
to  the  thrust -drag  relationship.  The  guidance  law  is 


Ivl  -  0 


1  »  1 
1  D 


Captive  Flight 

This  routine  Is  used  as  a  device  to  zero  out  coisputatlons  and 
printed  values  for  vehicles  which  are  in  captive  flight.  There  are 
three  taudes: 

o  Vehicle  2  locked  to  vehicle  1. 

o  Any  vehicle  (1,  2,  or  3)  may  be  locked  to  the  zero  origin 
doe.,  zero  position,  velocity,  and  acceleration). 

0  Vehicle  2  locked  to  vehicle  3. 


\ 
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Launch 

This  control  law,  which  simulates  the  launch-boost  phase  of  a 
missile  flight,  may  be  applied  to  any  of  the  three  vehicles.  The  call 
for  "launch"  is  usually  based  on  seme  criteria  stated  in  F0LICY  (e.g., 
range,  range  rate,  geometry,  accelerations,  time,  and — most  Importantly — 
combinations  thereof).  When  this  routine  is  called,  the  boast  velocity 
AV  must  be  specified  as  a  constant.  The  coimnanded  lateral  acceleration 
is  gravitational  only: 


=  -g  cos  Y  Ijj 


(41) 


Left  or  Right  Turn 

The  commanded  lateral  acceleration  vector  a^  is  in  the  horizontal 

plane  and  has  a  constant  value  as  specified  when  calling  the  routineCs) . 

The  two  routines  (left  and  right)  are  identical  except  for  an  algebraic 

sign  corresponding  to  the  direction  of  the  tuini  (±  1^).  It  was  decided 

to  specify  the  magnitude  of  the  turning  acceleration  in  terms  of  the 

resultant  normal  acceleration,  which  is  expressed  in  g's  or  F  /W  (see 

n 

Eq.  (21)).  Accordingly, 

Sq  -  g ^(F^/W)^  -  cos^  Y  (±  1^)  (42) 

CLOSED-LOOP  CONTROL  LAWS 


Proportional  Navigation 

The  commanded  lateral  acceleration  a  is  proportional  to  the  space 

i./ 

rate  of  rotation  of  the  LOS  between  missile  and  target.  Expressed  in 
vector  notation. 


ac  -  X  1^ 


(43) 


t 


-39- 


where  A  -  the  "navigation  constant"  (nay  be  treated  as  either  a  con¬ 
stant  or  a  variable) 

V  *  vehicle  speed 

=  relative  angular-rate  vector  as  defined  in  Section  III 
(o)^  =  In  Eq.  (5)) 

ly  *  unit  vector  along  the  missile  velocity  vector  V,  l.e., 

\-v/v 

The  direction  of  the  acceleration  is  defined  by 


1 ,  *  a>  y-  — 
1  r 


The  commanded  acceleration  a.  may  be  resolved  into  horizontal  and  ver- 
tlcal  components  by 


‘Ch 


CV 


(44) 


Missile  (X) 

This  routine  is  provided  as  an  example  of  how  to  incorporate  all 
significant  aerodynamic,  propulsive,  and  guidance  characteristics  of 
a  hypothetical  guided  missile  design  into  a  single  package.  Numerical 
values  which  are  unique  to  the  configuration  and  which  prestunably  will 
not  be  varied  may  be  Hated  in  the  routine  rather  than  supplied  as 
input  data  for  each  simulation  run  (e.g.,  reference  area  A,  initial 
weight  W^,  burning  rate  W,  etc.).  Moreover,  specialized  analytic  func¬ 
tions,  e.g.,  linear  or  polynomial  curve  fits,  nay  be  used  for  aero- 
dj'namic  C^,  a,  and  relationships  as  well  as  for  propulsion  charac¬ 
teristics.  The  guidance  law  for  Missile  (X)  is  a  modified  form  of  pro¬ 
portional  navigation;  again,  for  convenience,  all  necessary  relation¬ 
ships  are  incorporated  within  the  routine,  which  makes  it  unnecessary 
to  call  upon  the  proportional  navigation  routine.  For  further  details 
see  Appendix  C. 


VII.  CONCLUSIONS 


The  basic  framework,  organization,  input-output  integration,  flow, 
etc.  of  TACTICS  are  considered  complete.  However,  in  accordance  with 
its  purpose  as  a  research  tool,  it  is  open-ended  and  subject  to  adapta¬ 
tion  for  each  new  problem;  in  this  sense,  it  will  never  be  complete. 

This  adaptive  process  is  simple  and  flexible  because  of  the  available 
options  and  because  the  main  variables  defining  a  problem  can  be  treated 
externally  by  modular  units  (e.g.,  control- law  and  P0LICY  subroutines). 
Hundreds  of  simulation  problem  runs.  Involving  air-to-air  combat,  SAM, 
and  AS!!  applications,  have  been  performed  with  a  wide  variety  of  input- 
data,  control-law,  and  P0LICY  options.  Experiments  have  also  been 
performed  for  space  applications  using  the  two-body  equation  of  motion. 
There  is,  however,  no  guarantee  that  the  program  will  work  perfectly 
for  all  cases  in  spite  of  all  the  check-out  and  operational  experience. 
The  number  of  possible  configurations  and  combinations  is  extremely 
large,  and  it  is  unlikely  that  all  will  ever  be  tried.  As  with  all 
computer  problems,  skepticism.  Intuitive  reasoning,  and  cross-checking 
are  necessary. 


Part  2 


OPERATING  THE  PROGRAM 


^nz.  mrsmrcnes 


lErz  st/cdoas  stfafeb  sollou  saeold  be  osssliSered  aa  cpezzcics  rnaa&.l 
Saz  ebe  pzogssai.  'Sb*  reag^'s  ffgiriHlarltty  ultb  FsIsOUS  I?  is  jssaaed, 
fees:  ksawle£g/t  as  grogsMBclsg  ez?«siease  is-  acsessagy- 

SecS^OQ  3X  deals  trlth  lasnaiadsg  a  sdbroatiae  empased  o£ 

Isgicsl  scacewass  vaich  diceisirc  efea  ccaczol— lav  st£trcctizaes  if^esrsios 
sba  fligbc  of  each  vehicle,  as  eaglainaa  ia  Seccica  II.  Seceica  Z 
eolsiss  bew  £o  sec  cp  a  pcohlem  rca  hjr  ceadiag  ia  inicial-coodlicica 
£aca  sad  seLacds^  cptioas  for  iscegtacioo  merhoigs,  g:-Einccac,  cable 
valcses,  ecc. 

After  the  printed  cccpuc  sections  axe  explained  and  Ulc&crated 
is  Sectloa  XZ,  several  saiiple  problexs  are  described  ia  detail  ia  Sec- 
Cioa  XIZ.  lae  appeadlaes  coacaia  a  list  of  FG3I3A3!  lasCrcctlcss  for 
calling  optional  subroutines  (lacludiag  an  explanation  of  flags  and 
argument  variables),  listings  of  subroatlses,  and  an  explanation  of 
aerodynamic  tables  and  fomats.  Careful  study  of  the  illustrative 
exas^les  in  Section  XII  is  recosnended,  since  they  aay  serve  as  con¬ 
venient  guides  for  several  types  of  probless. 


PRECEDING  PAGE  BLANK 
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■a.  gcaagnASiaac  a  gfeSCT  SEnaHaggEss 


la  Sestica  H  crae  pssrpose  aaS  ftaasgioeaS  oseszeios  of  c&c  P^UOT 
sa&gpcslBe-  v^e  dg^cri&esS.  £si  sUcftrsicloQ  vxs  {ires  of  za  eleacatsrj 
flisSac  pSaa  tral'S-jr^  for  sascssrers  so^  ss  clsatg?,  scraijac  aad 

dive.  O&Tiosslj',  pollo^  cscj-tioss  based  era  sooe  crlceria  are  taecessary 
1EO  carry  ooc  ssch  a  flig&r  p.  »,  l.e.,  Co  dfccace  cbe  erassicfera  froos 
ersc  aaacsTer  co  asocber.  ‘Sbe  PCStlSAS  socadoc  for  cbe  elcoescary 
mazaesnrer  st£»ragfcines  s&eadcacd  above  is 

o  cmsi  (I,  Gs^,  USSt,  TI330 
o  Si3FLT  (I.  I£ES0,  Il£2) 
o  DIVEl  (I,  GF^SC,  IAEM,  USS) 
o  SZB51  (1,  GFiNCC,  lAE^,  ITEK) 

vbere  I  =  vehicle  Co  vhich  the  lav  applies  (1,  2,  or  3) 

IAE30  =  aa  integer  to  indicate  how  aerodynaaic  c<mmtations  are  to 
be  carried  out 

ITHR  *  an  integer  to  indicate  how  propulsion  cosputations  are  to 
be  carried  out 

GF0RC  =  nuanber  of  g's  as  shown  in  Fig.  8  and  given  by  Eq.  (21)) 

required  in  the  naneuver 

The  arguments  IAER0  and  ITHR  specify  whether  table  values,  analytic 
expressions,  or  other  alternatives  are  to  be  used  (a  complete  descrip¬ 
tion  is  given  in  Appendix  D). 

Assume  that  the  flight  plan  meat.oned  above  is  required  to  siruu- 
late  an  actual  takeoff  of  vehicle  1,  an  F-104  aircraft,  and  that  the 
criteria  that  might  be  used  are  the  following: 

o  Time;  TIME  (straight  flight  at  TIME  -  0) 
o  Speed:  V(l,4)  (sufficient  speed  for  climb) 
o  Flight-path  angle:  V(l,6)  (to  level  off) 
o  Altitude:  R(l,3)  (prior  to  turning) 
o  Heading:  V(l,5)  (proper  course  angle) 

The  FORTRAN  notation  for  position  and  velocity,  expressed  in  both 
Cartesian  and  spherical  coordinate  systems,  is  as  follows: 


BC1,1)  -  X, 

B(1.4)  -  jr,) 

2(1,2)  -  y. 

BO.5)  - 

2(1,3)  » 

2(1,6)  »  (5^ 

^a»l)  * 

V(l,2)  *  7^ 

Vfi.s)  - 

Vi 

V(l,3)  * 

^(1.6)  -  Y, 

All  ?Cil3AS  nat&tioa  for  posicioa,  relative  poslcica,  velocity, 
and  relative  velocity  is  parallel  to  the  ccove. 

A  typical  initial  P^LICif  statement  for  velu.cla  1  at  time  zero 
aight  thee  be 


CALL  STRFLT  (1,2,2) 

where  I  «  1  prefers  to  vehicle  1) 

IA£R^  =  2  \£efets  Co  aerodynanic  tables  ) 

ITHR  =  2  ('/cfers  to  military  thrusc  tables^) 

The  successive  statements  might  then  be  the  following: 

Climb  at  0.5  g  to  20-deg  flight-path  angle 
o  IF  (V(l,4)  .GT.  120.0)  CALL  CLIMBl  (1,0. 5, 2, 2) 
o  IF  (V(l,6)  .GT.  20.*RAD)  CALL  STRFLT  (1,2,2)* 

Begin  leveling  off  at  1000  ft 
o  IF  <R(1,3)  .GT.  1000.0)  CALL  DIVE1(1,0.5,2,2) 

When  within  0.1  deg  of  level  flight,  make  a  0.5-g  turn  to  a 
heading  of  150.0  deg  and  resume  straight  flight 
o  IF  0/(1, 6)  .LT.  0.1*RAD)  CALL  RTRN1(1,0.5,2,2) 
o  IF  (V(l,5)  .GE.  150.0*RAD)  CALL  STRFLT  (1,2,2) 

This  elementary  illustration  defines  a  policy  and  flight  path  for 
one  vehicle.  For  intercept  trajectwry  problems,  policies  for  two  or 


^F-lOA  tables  will  be  loaded  for  vehicle  1. 

^The  symbol  RAD  is  used  to  convert  degrees  to  radians. 
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misre  vcSiicles  tacsc  be  sfm&larly  defl-zei.  Saiistlsg  %  rcfereace  ta  a. 
v^iicie  Cl  »  1,  2.  or  3)  in  ?tfLICT  resoles  in  «  zero  acederadoa 
de£is£c£cn  for  e£:ac  vehicle,  vlsicl!  aoi  be  appzcrprlete  far  cosscant- 
or  zero-vclocirj  Ce.g..  sroaao-eargec)  cases.  Seccion  XII  presents 
more  tjpical  ana  empiex  eicazpies  oi  ?![fLsC7  sobroetiaes  tfich  sample 
probiess. 

5o  far  this  sf&Jecc  has  been  discussed  In  the  context  of  choices 
to  be  made  wichia  the  ?#£lCf  routine.  SCucerous  other  options  nay  be 
selected  by  reading  in  flags  or  constants  as  part  of  iaitial-condltisn 
data.  Bovever,  it  is  sometimes  desirable  to  overrld't  these  initial 
instructions  in  P#LIC¥  if  during  a  problea  run  a  situation  arises  that 
requires,  perhaps,  a  change  in  frequency  of  printout  or  Integration 
step  size.  The  foUouing  subsections  list  possible  options  in  addition 
to  these  aentioned  above. 

STOPPIiiG  THE  PgQGRAM 

The  ptograa  automatically  temlnatcus  after  finding  the  closest 
miss  distance  between  vehicles  2  and  1  or  2  and  3  unless  otherwise 
specified  by  the  flag  IHISS.  If  no  isissile  is  Isv^ched,  the  program 
stops  when  ruimlng  time  (TIME)  becomes  greater  than  maximum  specified 
tine  (T0TAL);  this  value  is  set  in  initial-condition  data  or  in  P0LICY. 

If  the  program  is  to  continue  after  finding  the  miss  distance, 
flag  IMISS  =  1  must  be  set  in  P0LICY.  This  is  used  when  a  second  mis¬ 
sile  is  to  be  launched  or  if  the  vehicles  are  to  continue  on  their  tra¬ 
jectories  (IMISS  must  be  reset  for  every  launch).  For  example,  if  the 
user  wants  the  program  to  stop  after  finding  the  second  missile's  miss 
distance,  IMISS  must  be  reset  to  zero. 

In  summary: 

o  IMISS  >  0.  The  program  finds  the  closest  miss  distance  of 
the  missile  and  then  stops. 

o  IMISS  ■  1.  The  program  continues  after  finding  the  miss  dis¬ 
tance.  IMISS  must  be  reset  for  each  launching. 

o  If  no  missile  Is  launched,  the  program  stops  when  time  is 
greater  than  T0TAL  (DATA  64) . 
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ISIECRAnOS 

Four  different  cypes  of  intregrstlon  can  be  used  (see  Section  XIZI). 
The  selection  is  aade  by  netting  the  JINIEG  flag  as  an  initial  condition 
(DAIA  122) . 

If  JINIEG  *  0,  varlable-ncep  Adans-Sfoulton  predictor-corrector 
integration  is  used.  Oth  iz  values  vhich  are  necessary  vhen  using  this 
type  of  integration  are  ESJEST  and  HKIK.  ERIEST  (DAX&  123)  is  the 
truncation  error  test  for  variabls  integration;  if  not  specified  by 
the  user,  it  is  autonatically  set  at  l«0E-O5.  HHIN  (DATA  135)  is  the 
ainisua  step  size  for  the  integration.  The  progras  sets  IfiflH  *  DI0 
(the  initial  integration  step  set  in  DATA-  136)  imless  otheivlse 
specified. 

If  JINTEG  1,  fixed-step  Kunge-Kutta  integration  is  used;  if 
JINIEG  «  2,  fixed-step  Adsms-^ulton  integration  is  used.  In  these 
cases  the  step  size  depends  on  the  value  read  in  for  DI0.  If  JINTEG  « 

3,  the  progras  integrates  on  a  variable  step  size,  controlled  to  allnw 
printout  exactly  at  specified  Inteanrals. 

In  sunmary: 

o  JINTEG  »  0;  Variable-step  Adams-Houlton  integration, 
o  JINTEG  «  1:  Fixed-step  Runge-Kutta  Integration, 
o  JINTEG  ■  2;  Fixed-step  Adams-Moulton  integration, 
o  JINTEG  “  3;  Variable-step  Adams-Moulton  integration  with 
exact  printout. 

NUMBER  OF  VEHICLES  USED 

If  the  problem  does  not  use  one  or  more  of  the  vehicles,  its 
printout  values  can  be  set  to  zero  by  calling  CAPFLT(I,M0DE) ,  where  I 
indicates  which  vehicle  is  not  to  be  used  and  MODE  is  set  at  2.  (See 
Example  3  in  Section  XII.) 

TIME  LAGS 

Subroutine  LAG  represents  a  transfer  function  simulating  the  time 
lag  between  input  command  AC0M  and  output  response  A0UT.  This  function 


Bay  vary  considerably  in  cosplexiCy.  froa  a  sickle  one-to-one  corre¬ 
spondence  to  a  highly  cosrplicated  presentation  of  a  nissile  guidance 
and  control  loop.  The  degree  of  rseiisn  (and  hence  coaplexity) 
depends,  of  course,  on  the  specific  problem  and  Its  significance  to 
the  final  results.  Accordingly,  the  LAG  subroutine  should  be  considered 
a  flexible  aodule  which  can  be  izodified  for  a  particular  problea.  A 
currently  available  LAG  routine  represents  the  tise  response  of  a  vehi¬ 
cle  as  the  product  of  as  many  as  three  first-order  (exponential)  time 
lags.  Unless  specified  by  the  user,  no  time  lag  is  used;  A0UT  AC0M 
unless  structurally  (A5MAX)  or  aerodynamicaliy  (CI^lAX)  constrained. 

In  order  to  introduce  time  delay,  flag  TAU(l)  must  be  set  in  initial- 
conditions  data  to  indicate  the  number  of  time  lags  desired,  and 
TAU(I,J)  must  be  set  equal  to  the  value  of  the  time  lags.  See  DATA 
65-78  in  Appendix  E  for  details  on  reading  the  time  constants.  Example 
3  in  Section  XII  is  an  illustration  of  the  way  a  time  lag  is  used  for 
a  SAH. 

RECALL 

This  feature  enables  the  recall  of  the  missile  once  it  has  been 
launched.  For  details  on  the  use  of  the  recall  option,  see  Fig.  24 
and  Example  4  In  Section  XII. 

RESTORE 

The  restore  option  enables  the  user  to  restore  all  numerical 
values  existing  at  launch  time  if  a  hit  or  miss  has  occurred.  See 
Fig.  27  and  Example  5  in  Section  Xll  for  details  on  its  use. 

ROUND  EARTH.  ROTATING  OR  NONROTATING 

The  FORTRAN  flags  IRT,  IRF,  IR0T8,  IPRINT(I),  INERF,  and  INERT 
are  used  to  select  the  round-earth  options.  In  defining  ^jatellite 
motion,  velocity  components  are  usually  given  with  respect  to  a  non- 
rotating,  inertial  coordinate  frame.  However,  in  defining  aircraft 
or  missile  motion  it  is  convenient  to  use  velocity  components  expressed 
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:<lth  respect  to  a  local  earth-fixed  rotating  fratse  of  reference.  Of 
course,  if  the  earth  is  considered  as  nonrotating,  there  Is  no  difference. 

These  options  are  controlled  as  follows: 

o  IR0T8  ^  0:  Nonrotating  earth. 

11^8  »  1:  Rotating  earth. 

o  INERF  *  0:  Vehicle  1  (fighter)  velocity  expressed  In  relation 
to  local,  earth-fixed,  rotating  frante. 

INERF  *  i:  Vehicle  1  velocity  expressed  in  relation  to  local, 
earth-fixed  nonrotating  Inertial  fraine. 

o  INERT  >  0  or  1:  Corresponding  statements  apply  to  vehicle  3 
(target)  velocity. 

The  location  of  the  origin  of  the  local  earth-fixed  frame  is 
specified  by  Initial-condition  data  in  terms  of  latitude  LATO  and 
longitude  L0NGO.  This  is  essential  to  properly  account  for  rotational 
effects.  The  x-,  y-,  and  z-axes  of  this  coordinate  frame  are  taken  to 
form  a  right-handed  system  oriented  in  the  following  way:  The  z-axls 
is  coincident  with  the  local  vertical  or  radius  of  the  earth  at  the 
point  LATO,  L0NGO;  the  y-axls  is  in  the  local  horizontal  plane  directed 
eastward;  and  the  x-axls  is  directed  southward  or  northward.  Initlal- 
coadltlon  position  data  for  each  vehicle  must  be  supplied  in  terms  of 
this  local  frame  of  reference,  and  thei'e  are  three  coordinate  systems 
which  might  be  used: 

o  Cartesian  (x,  y,  z)  coordinates  with  respect  to  the  origin 
defined  by  lATO,  L0NGO  (set  IRF  or  IRT  -  0  and  use  DATA  2-4 
or  41-43  respectively) . 

o  Spherical  (r,  0,  <p)  coordinates  with  respect  to  the  origin 
defined  by  LATO,  L0NGO  (set  IRF  or  IRT  -  1  and  use  DATA  7-9 
or  46-48  respectively) . 

o  Geocentric  latitude,  longitude,  and  altitude  above  a  reference 
p.pheroid  (set  IRF  or  IRT  -  2  and  use  DATA  111-113  or  116-118 
respectively) . 

Mathematical  details  and  derivations  pertaining  to  the  coordinate  trans¬ 
formations  and  equations  of  motion  are  contained  in  Appendix  B. 
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X.  IKPUT  FORM  AND  ORDER  0?  INPUT 

Input  data  include  various  problem  constants,  parameters ^  control 
flags  (for  options),  and  initial-condition  values  for  position  and 
velocity.  There  are  two  sections  of  data:  (1)  the  aerodynamic  tables 
for  specific  aircraft  or  missiles  and  (2)  the  main  sec,  which  gives 
initial  conditions,  flags,  and  constants.  An  input  form  is  shown  in 
Fig.  9  for  convenience  in  specifying  the  main  set  of  data.  A  brief 
explanation  is  given  here,  but  it  is  likely  that  the  examples  given  in 
Section  XII  will  more  clearly  illustrate  the  use  of  the  form.  Each 
line  refers  to  a  data  card  (the  examples  will  show  that  in  general 
most  of  the  lines  and  spaces  may  be  left  blank) . 

The  first  line,  l.e.,-  the  first  data  card,  specifies  whether 
aerodynamic  tables  are  to  be  used  for  a  vehicle.  It  gives  values  to 
the  flags  JVEH(I),  I  =  1,2,3.  If  JV£H(I)  »  1,  tables  for  vehicle  (I) 
will  be  used;  if  JVEH(I)  =  0,  tables  for  vehicle  (I)  will  not  be  used. 
(If  not,  the  appropriate  aerodynamic  and  propulsion  constants  may  be 
entered  in  DATA  80-94.)  The  spacing  for  the  JVEH  card  is  as  shown  in 
Fig.  10.  If  cables  are  to  be  used,  they  are  read  immediately  following 
the  JVEH  card,  as  shown  in  Fig.  11. 

A  title  card  describing  the  run  (F0RMAT  12A6)  is  placed  between 
the  aerodynamic  tables  and  the  main  set  (this  card  is  required) .  The 
locations  (l.e>,  data  numbers)  and  descriptions  of  all  input  data  are 
given  in  Appendix  E.  However,  the  input  form  has  spaces  only  for  those 
locations  most  likely  to  be  used^  with  spare  spaces  available  at  the 
bottom  of  the  page.  (There  are  provisions  for  200  data  locations,  but 
only  143  are  currently  being  used.) 

AERODYNAMIC  AMD  PROPULSTc-  TABLES 

The  user  is  •  co  approximate  the  aerodynamics  necessary  Co 
describe  the  flight  of  an  aircraft  by  means  of  tables.  If  tables  are 
used,  they  furnish  such  information  as  angle  of  attack  and  drag 

* 

Figure  9  shows  the  input  form  reduced  in  size.  A  full-sized  ver¬ 
sion  is  included  at  the  end  of  the  Memorandum  for  the  reader’s  use. 


Title  card 

A/tain  set  of  data  (DATA  1-143) 

Data  for  second  run 
JVEH  flag  card 

Tables  (if  different  from  first  run) 

Title  card 

Main  data  that  differs  from  that  of  first  run 


Data  for  third  run,  elc. 


Fig.  1 1  -  Order  of  input 
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coefficlent  as  a  function  of  lift  coefficient  and  Mach  nundber  and  thrust 
as  a  function  of  altitude  and  Mach  number. 

When  calling  a  control-law  subroutine,  the  flags  IAER0  and  ITHR 
are  used  to  select  the  aerodynamic  and  propulsion  operations  pertaining 
to  the  choice  of  cable  values,  analytic  functions,  or  constants.  In¬ 
structions  for  calling  these  options  are  given  In  Appendix  D. 

Appendix  F  describes  the  form  and  organization  of  aerodynamic  and 
propulsion  tables.  A  complete  set  of  tables  is  also  shown  as  an  example 
In  Fig.  35,  Appendix  F.  Hie  model  atmosphere  Is  calculated  by  analytic 
expressions  (from  Ref.  4)  as  given  in  Appendix  I  (there  are  no  provisions 
for  using  table  values  for  model  atmosphere  representation) . 

AERODYNAMICS  AND  PROPULSION  WITHOUT  TABLES 

There  are  several  possibilities  for  handling  aerodynamic  or  pro¬ 
pulsion  calculations  for  the  JVEH(I)  «»  0  option: 

o  Constants  may  be  read  in  the  appropriate  locations  (80-94)  and 
Eq.  (27)  used  for  C^,  C^,  and  a  computations  by  specifying  the  argument 
IAER{f  *=  1  when  calling  the  control-law  subroutine,  (See  Appendix  D 
for  instructions.) 

o  Constants  may  also  be  used  In  appropriate  locations  (95-109) 
to  specify  values  of  propulsion  and  fuel  flow  characteristics,  e.g., 
thrust,  specific  impulse,  burnout  weight,  or  boost  acceleration.  Thrust 
is  set  equal  to  the  data  value  by  using  the  ITHR  flag  as  explained  In 
Appendix  D  for  ITHR  ■  3,  4,  or  5.  Vehicle  weights  are  determined  from 
the  expression 

W  ««  Wo  -y*  W  dt 

where  W  may  be  a  value  determined  from  fuel  flow  tables  (ITHR  ■  1,  2) 
or  a  constant  input  value  as  specified  In  data  locations  144-146  when 
using  the  ITHR  >•  3,  4,  5  option.  However,  In  order  to  simulate  a  con¬ 
stant-acceleration  boost  phase  of  a  missile,  the  program  will  automat¬ 
ically  calculate  vehicle  weight  from  the  expression  W  ■  W  exp  a_(t-t,  )/gI 

O  D  L 
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where  W  »  the  current  weight  at  time  t 
*  the  Initial  weight 
t  >>  the  launch  time 

I4 

»  the  boost  acceleration,  considered  to  be  an  average  value 

D 

I  B  the  specific  impulse  of  the  rocket  motor 
This  automatic  alternative  computation  is  initiated  when  an  input  data 
value  for  ag  is  supplied  in  locations  99,  104,  and  109  (as  applicable 
to  vehicle  1,  2,  or  3) .  Provisions  are  also  made  for  supplying  missile 
burnout  weights  in  locations  98,  103,  and  108,  but  no  computations  are 
performed  on  these  quantities  within  the  main  body  of  the  program  itself; 
the  locations  are  provided  merely  as  a  convenience  for  formulating  sub¬ 
routines  requiring  this  form  of  input  data. 

o  For  those  cases  where  aerodynamic  computations  for  a  particular 
vehicle  are  not  significant  to  Che  problem,  it  is  not  necessary  to  read 
in  either  tables  or  constants  for  that  vehicle.  The  IAER0  ■  3  argument 
is  used  (see  Appendix  D)  and  the  vehicle  is  assumed  to  be  a  point  mass 
moving  at  constant  speed  (but  not  necessarily  at  zero  lateral  acceleration). 

o  For  specialized  aerodynamic  or  propulsion  characteristics  (e.g.. 
Missile  (X),  Appendix  C)  analytic  expressions  and  necessary  constants 
may  be  Incorporated  into  the  control-law  subroutine. 

MAIN  SET  OF  DATA 

The  main  set  of  data  specifies  such  initial  conditions  as  vehicle 
position  and  velocity,  aerodynamic  constaitt;s,  structural  and  attitude 
limits,  time  lags,  program  flags  (see  Appendix  H)  and  constants.  Data 
numbers  with  corresponding  program  variables  and  descriptions  are  given 
in  Appendix  E.  There  are  143  separate  entry  spaces  in  the  main  data. 

If  a  value  is  not  read  into  a  space,  that  entry  is  automatically  taken 
to  be  zero;  therefore,  only  nonzero  data  need  be  specified. 

The  format  for  reading  in  the  main  data  is  (Al,  13,  5E14.8).  The 
initial  value  on  the  card  is  the  number  of  the  first  data  entry  on  the 
card.  Entries  which  follow  this  value  on  the  card  must  be  in  sequence. 

A  minus  sign  is  placed  in  the  first  column  of  the  last  data  card  to  in¬ 
dicate  that  the  entries  for  that  case  are  finished.  Any  data  following 
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this  card  are  £or  a  new  case  (see  Fig.  10).  For  sample  data,  refer  to 
the  numerical  exeimples  in  Section  XII. 


Initial  Position  and  Velocity 


As  shown  on  the  input  form,  data  spaces  are  available  for  those 
flags  that  indicate  how  the  position  and  velocity  of  the  vehicles  are 
to  be  entered:  in  spherical  or  Cartesian  coordinates  or  (in  the  case 
of  position  only)  in  latitude  and  longitude.  Data  entries  which  follow 
are  for  specifying  the  positions  and  velocities  of  vehicles  1  and  3. 
Vehicle  2  is  initially  assumed  to  be  a  missile  attached  to  vehicle  1; 
however,  data  locations  30  through  35  (not  shown  on  the  input  form) 
are  available  for  setting  initial  position  and  velocity  conditions 
separately  for  vehicle  2.  That  is,  if  data  entries  are  made  in  any 
one  of  these  locations  (30-35) ,  TACTICS  will  start  simultaneously  com¬ 
puting  the  trajectories  of  all  three  vehicles;  otherwise,  vehicle  2  is 
initially  attach»'.d  to  1  ?!nd  computations  will  be  performed  for  only 
two  vehicles  (unless  or  unti]  there  is  a  P0LICY  call  for  launching  2) . 


Aerodynamic  Constants 


If  analytic  functions  are  to  be  used  for  aerodynamic  computations, 
the  equation  constants,  i.e.,  those  applicable  to  Eqs.  (26)  and  (27), 
must  be  entered.  Data  spaces  are  assigned  for  such  constants,  as  shown 
on  the  input  form. 


Other  Data 

Spaces  for  structural  and  aerodynamic  constraints  are  shown  on  the 
input  form.  Maximum  elevation  and  azimuth  angle  gimbal  limits  may  be 
specified  by  DATA  138-1A3.  Up  to  three  consecutive  vehicle  time  lags 
can  be  entered  (see  Section  V).  Program  flags  and  constants  should  be 
set;  typically,  these  would  specify  minimum  integration  step,  type  of 
integration,  and  total  time  of  run.  If  certain  necessary  values  for 
program  operation  are  not  entered,  TACTICS  will  automatically  assume 
"default"  values,  print  an  Informative  message,  and  continue. 


c 
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RUNNING  MORE  THAN  ONE  CASE 

Any  manber  of  consecutive  cases  can  be  submitted  for  a  single  run. 
As  indicated  previously,  ihe  last  card  of  the  data  for  a  case  has  a 
minus  sign  in  the  first  column  indicating  the  end  of  the  case.  To 
submit  a  second  case,  only  data  differing  from  values  in  the  first 
case  need  be  entered.  For  each  case  a  JVEH  card  and  comment  card  must 
be  entered  before  the  main  data.  With  reference  to  Fig.  10,  note  that 
although  the  JVEH  card  contains  all  zeroes,  the  same  tables  are  used 
in  the  second  case  as  in  the  first. 


XI.  OUTPUT 


The  output  from  a  typical  computer  run  consists  of  labeled  and 
unlabeled  Initial-condition  data  values  and  optional  sections  of  print¬ 
out.  Initial-condition  values  will  automatically  be  printed,  but  the 
printing  of  the  optional  sections  must  be  specified  In  F0LICY.  The 
printing  of  unlabeled  data  consists  of  a  listing  of  all  locations, 

DATA  1-200,  as  shown  In  Fig.  12.  Labeled  Initial  conditions  for  these 
values  are  as  shown  In  Fig.  13.  Table  1  lists  sets  of  available  op¬ 
tional  output  sections,  the  most  Important  of  which  Is  the  main  set 
or  "standard  output"  giving  position,  velocity,  acceleration,  weight 
information,  and  other  basic  quantities  printed  at  specified  time 
Intervals;  these  intervals  are  controlled  by  the  variable  DTP0  (DATA  63). 
The  other  sections  concern  information  about  aerodynamics,  attitude, 
etc.  as  Indicated  in  Table  1.  In  order  to  specify  these  optional  sec¬ 
tions  in  P0LICY,  the  variable  NPRINT  is  set  equal  to  the  required  num¬ 
ber  of  output  sections. 

A  typical  output  specification  In  P0LICY  would  be  the  following: 

NPRINT  ■  3:  Three  sections 

IPRINT(l)  •  1;  Standard  output 

IPRINT(2)  •  2;  Aerodynamics 

IPRINT(3)  ■  3:  Attitude  angles 

A  sample  page  of  output  for  the  NPRINT  *  3  option  is  shown  in 
Fig.  14.  Most  of  the  items  shown  are  self-explanatory  except  perhaps 
the  following: 

1.  The  integration  step  size  shown  at  the  top  of  Fig.  14  is  the 
last  integration  step  taken  before  printout.  If  the  JINTEG  ■  J  mode 
of  Integration  is  being  used,  two  values  will  be  printed:  the  step 
taken  to  reach  printout  time  (TP0)  and  the  step  the  routine  would  have 
used  If  there  had  not  been  an  Immediate  print  requirement.  (The  step 
taken  to  reach  printout  must  necessarily  be  the  smaller  of  the  two.) 

2.  The  units  used  are  feet  for  distance,  degrees  for  angles, 
pounds  for  force  and  weight,  g's  for  acceleration,  and  radians  per 
second  for  angular  rates. 
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Fig.  13  —  Sample  initial  data  values  (labeled) 
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3.  GAMD0T  and  THD0T*C0S  (GAM)  are  the  angular  rates  of  rotation 
of  the  vehicles'  velocity  vectors. 

4.  0MEGAR  refers  to  the  angular  rates  of  rotation  of  the  LOS 
RREL(I,J).  0MEGAB  output  is  provided  for  guidance  laws  that  use  a 
biased  angular-rate  term  (see  Appendix  C) . 

5.  An  EXTR  F0RMAT  (E16.8)  is  provided  for  six  extra  quantities 
(middle  right  of  standard  output  package).  This  Is  extremely  useful 
for  debugging  and  printing  additional  Information. 

6.  A  CLMAX  or  ASMAX  print  (adjacent  to  acceleration  quantities) 
will  occur  whenever  these  respective  limits  are  exceeded.  If  no  in¬ 
put  data  value  has  been  given  to  CLMAX  or  ASMAX,  these  values  are 
assumed  to  be  infinite  (e.g.,  10°). 

7.  PHRREL  and  THRREL  (bottom  right)  refer  to  the  orientation 
angles  cp  ax'.d  6,  respectively,  of  the  LOS,  i.e.,  RREL(I,J). 


Table  1 

OUTPUT  SECTIONS 


Section 

Title 

Description 

1 

Standard  output 

Position,  velocity,  acceleration,  and 
other  basic  quantities 

2 

Aerodynamics 

Lift,  maximum  lift  and  drag  coefficients, 
angle  of  attack,  lift,  and  drag 

3 

Attitude  angles 

Roll,  bank,  elevation,  azimuth,  and 
bearing  angle 

4 

Round-earth 

coordinates 

Latitude,  longitude,  and  altitude 

5 

Initial  conditions 

Input  data 

6 

Atmosphere 

Air  pressure,  temperature,  air  density, 
speed  of  sound,  and  Mach  number 

9 

Angular  rates 

Angular-rate  vectors  and  unit  vectors 

^A’  \lv 

NOTE:  Sections  7,  8,  and  10  through  18  are  blanks  available  to  the 
user  if  different  types  of  output  are  desired. 
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XII,  EXAMPLES 

Five  examples  indicating  Che  type  of  program  which  can  be  run  on 
TACTICS  and  the  way  Che  program  can  be  used  are  given  in  this  section. 
Each  problem  is  defined  by  Initial-condition  input  data  and  by  a  F0LICY 
subroutine  which  dictates  the  control  laws  governing  the  flight  of  each 
vehicle.  The  examples  show  how  a  F0LICY  subroutine  is  set  up  and  what 
input  data  are  necessary  for  a  particular  problem.  The  exan^les  se¬ 
lected  cover  the  different  options  available  in  TACTICS  by  employing  a 
number  of  the  special  features  included  in  the  program  (e.g.,  restore 
and  recall) .  These  examples  range  from  the  simple  case  of  vehicles 
flying  maneuvers  with  no  missile  to  the  more  compllcrited  one  of  launch¬ 
ing  two  missiles. 

EXAMPLE  1:  AIRCRAFT  MANEUVERS.  NO  MISSILE 

The  flights  of  two  vehicles  in  different  maneuvers  ace  simulated; 
no  missile  is  Involved.  In  setting  up  this  example^  as  with  all  the 
others,  two  steps  are  necessary;  (1)  developing  the  F0LICY  subroutine 
and  (2)  entering  initial-condition  data. 

P0L1CY  Subroutine 

The  following  POLICY  statements  define  the  problem: 

Vehicle  1  (Constant-Speed  Aircraft:  No  Aerodynamic 
or  Propulsion  Computations) 

0  Straight  flight. 

0  After  0.1  sec,  pull  a  4-g  climb. 

o  When  vehicle  has  climbed  30  deg,  level  off  to  horizontal  and 
continue  on  straight  flight. 

Vehicle  2 

o  Not  used. 
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Vehiele  3  (Constant-Speed  Aircraft:  Ho  Aerodynamic 
or  Propulsion  Computations} 

o  Fly  straight  flight. 

o  After  1  sec,  perform  a  barrel  roll  pulling  4  g's  and  rolling 
60  deg/sec. 

o  Uhen  roll  is  coapleted,  perfors  4-g  diving  turn  with 
130-deg  roll. 

These  stateiaents  are  translated  into  FORTRAN  expressions  to  forau- 
late  the  PtLlCt  subroutine.  Figure  15  is  a  listing  of  the  routine. 

The  type  of  output  desired  is  first  specified  in  POLICY  (see 
Section  XI).  In  this  example,  only  the  standard  output  and  attitude 
angles  are  to  be  printed,  since  no  aerodynamic  or  propulsion  computa¬ 
tions  are  being  performed. 

Commands  governing  the  first  vehicle's  flight  are  given.  There 
are  different  control  lavs  and  criteria  for  changing  maneuvers.  For 
each  maneuver  the  type  of  aerodynamic  and  thrust  computation  must  be 
specified  in  the  argument  following  the  maneuver  name  (see  Appendix  D) . 
For  this  example,  no  aerodynamic  or  propulsion  computation.^  ero  in¬ 
volved,  so  IAER0  »  3  and  ITHR  *5. 

CLIMBl  (simple  climb)  is  called,  specifying  4  g's  in  the  argument, 
until  vehicle  1  has  climbed  30  deg.  The  criterion  used  in  this  case 
is  V(l,6),  the  flight-path  angle  STRLVL  is  called  to  level  off 

the  vehicle.  Flag  LEVEL(l)  =  2  Indicates  that  the  aircraft  is  hori¬ 
zontal  again. 

Since  the  second  vehicle  is  not  being  used,  CAPFLT(2,2)  is  called, 
which  sets  all  values  pertaining  to  the  vehicle  equal  to  zero. 

The  motion  of  the  third  vehicle  is  defined  by  calling  STRFLT 
(IAER0  =  3,  ITHR  “  5)  and  then  changing  to  BRLRLl  (barrel  roll).  In 
the  argument  listing  for  BRLRLl,  the  number  of  rolls  (1),  the  g’s  to 
be  pulled  (4),  and  the  roll  rate  (60  deg/sec)  are  specified.  Flag 
IR0LL(3)  =  2  Indicates  that  the  number  of  rolls  required  Is  completed. 
This  is  used  as  the  criterion  in  P0LICy  for  switching  to  subroutine 
RTRN5,  a  4-g  diving  turn  with  a  13G-deg  roll.  See  Appendix  D  for 
further  instructions  on  the  calling  of  individual  maneuvers. 
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C  SAMPLE  POLICY  FOR  FLYING  AIRCRAFT  MANUEYERS*  NO  MISSILE 

r 

C****  SPECIFY  OUTPUT 
NPRINT=2 
IPRINT(1)=1 
IPRiNT(2)>:3 
C 

C  FIGHTER  COHMANOS  ********** 

GO  TO  f 110tl20«130tl40)tJPOL 
I 10  CONTINUE 

IF  ITINE  .G  O.i)  CO  TO  120 


CALL  STRFLTI1»3*51  PFIGHTE 

GO  TO  190  SEC,  I 

120  CONTINUE 

IF  fVfl,6t  «GT.  30.0«RA0I  60  TO  130 
CALL  CLIMBlf 1.4«0,3,5I  ♦FIGHIE 

JP0L»2  CLIMBE 

GO  TO  190  BEGINS 

130  CONTINUE 

IF  ILEVELII}  .EC«  2)  GO  TO  140 
CALL  STRLVLllt3,5,LEVELI  AWHEN  L 

JPCL»3  OFF  ro 

GO  TO  190  FLIES 

140  CONTINUE 

CALL  STRFLTflt3«5} 

JP0L»4 
190  CONTINUE 


♦FIGHTER  FLIfS  STRAIGHT  FLIGHT  FOR  .1 
SEC,  2ER0  THRUST 


♦FIGHTER  PULLS  4  G  CLIMB  UNTIL  IT  HAS 
CLIMBED  30  PEG,  AT  WHICH  TIME  IT 
BEGINS  LEVELING  OFF 


♦WHEN  LEVELin«2,  VEHICLE  HAS  LEVELED 
OFF  TO  H0RI20NTAL  iVIIeSI^O)  AND 
FLIES  STRAIGHT 


i!  1 


c  **********  missile  COMMANDS  ********** 

GO  TG  1210,220,2301, KPOL 
210  CONTINUE 

call  CAPFLTIZ,2}  AMISSILE  not  being  used  so  QUANTITIES 

GO  TO  290  ZEROED  OUT 

220  CONTINUE 
230  CONTINUE 
290  CONTINUE 
C 

c  ***^******  target  commands  ********** 

GO  TO  (310,320, 330), LPOL 
310  CONTINUE 

IF  (TIME  .GE.  1.0)  GO  TO  320  ATARGET  FLIES  STRAIGHT  FLIGHT  FOR 
CALL  STRFLTI3,3,5)  1  SEC 

GO  TO  390 

320  CONTINUE  ATARGET  THEN  MAKES  ONE 

CALL  BRLRL1(3,1.0,IR0LL,4.0,60.0,3,5)  BARREL  ROLL  (4G,  ROLL 

IF  UR0LL(3)  .EO.  2)  GO  TO  330  60  OEG/SEC).  iR0LL(n«2 

LP0L«2  INDICATES  ROLLS  ASKED 

GO  TO  390  FOR  ARE  COMPLETED 

330  CONTINUE 

CALL  RTRN5(3,4.0, 130.0,3,5)  ATARGET  PULLS  4  G  DIVING  TURN 
LP0L=3  V»ITH  130  DEG  ROLL 

390  CONTINUE 


RETURN 

END 


Fig.  15  —  Semple  P0.ICY  sobrouHne  for  aircraft  maneuvers, 

no  missiles 


1 


Initial  Data 
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After  the  P0LIcy  subroutine  is  written,  the  input  data  for  setting 
up  the  initial  conditions  of  the  problea  aust  be  read  in.  The  initial 
position  and  velocity  of  the  vehicles  aust  be  specified.  In  this  ex- 
aaple,  vehicle  1  la  situated  at  the  origin  at  an  altitude  of  15,000  ft. 
Vehicle  3  is  8000  ft  down  the  y-axls  at  an  altitude  of  20,000  ft.  The 
velocity  of  vehi.cle  1  is  1340  ft/sec;  it  is  headed  down  the  -x-axls 
(V(l,5)  »  180.0  deg).  Vehicle  2  has  a  velocity  of  1000  ft/sec,  and  it 
is  headed  down  the  y-axls  (V(3,S)  »  90  deg).  See  Fig.  16.  Other  data 
used  in  this  exai^le  are  the  following: 

Starting  value  for  integration  step  size  (DT0)  «  0.01. 

Structural  lateral  acceleration  liad.t  of  vehicles  >  7.3,  7.5. 

Area  of  vehicles  (ASFA)  «  196.0,  385.0. 

Initial  weight  of  vehicles  (UO)  >  16,669.0,  33,283.0. 

Time  Interval  for  printing  output  (DTP0)  *  0.5. 

Tlae  value  at  which  program  is  to  stop  (T0IAL)  «  25.0. 

Figure  17  shows  how  the  data  are  esitered  on  the  input  font. 

EXAMPLE  2:  LAUHCHIXG  OHE  MISSILE 

The  first  vehicle  (fighter)  is  to  pursue  the  third  vehicle  (tar¬ 
get)  until  a  launching  position  is  reached.  At  this  tine,  the  missile 
is  launched  frou  the  fighter  and  the  program  is  to  find  the  closest 
miss  distance  of  the  missile  and  then  termlnste.  Aircraft  tables  are 
to  be  used  to  simulate  the  aerodynamics  of  vehicles  1  and  3. 

P0LICT  Subroutine 

The  following  P0LICY  statements  define  the  problem: 

Vehicle  1  (F-104  Intercepting  Aircraft) 

o  Fly  lead  colllaion  navigation  course;  military  chrust. 

If  the  range  to  target  is  less  than  6500  ft  and  the  angle  off 
the  target's  tail  is  less  than  30  deg,  launch  missile. 


o 


Fig.  16  — Three-<«meniior«l  diogrom  of  imtial  conditions  for  Exomple  1 
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Fig.  17  —  Input  data  for  Example  1 


o  Pull  constant-Mach-numbert  constant-altitude  right  turn; 
thrust,  afterburner. 

Vehicle  2  (Proportional  Navigation  Missile) 

o  Fly  captive  flight  until  launch  criterion  is  satisfied. 

o  Launch,  boost,  fly  unguided  and  then  guided  in  accordance 
with  guidance  and  aerodynamic  characteristics  specified  in 
special  missile  subroutine. 

o  When  range  rate  (missile-target)  becomes  greater  than  zero, 
initiate  process  for  finding  miss  distance  and  end  program. 

Vehicle  3  (F-IOS  Target  Aircraft) 

o  Fly  straight  and  level;  thrust,  80-percent  throttle  setting, 
military  power. 

o  lifhcn  missile  is  launched,  perform  5.5-g  constant-Mach-number 
left  turn;  thrust,  afterburner. 

See  Fig.  18  for  the  actual  P0LICY  subroutine.  The  output  speci¬ 
fied  is  standard  output,  aerodynamics,  and  attitude-angle  sections. 

LEADCL  (lead  collision)  is  called  for  the  fighter,  with  IAER0  =“  2, 
indicating  tables  for  aerodynamic  functions,  and  ITHR  “  2,  indicating 
tables  for  military  thrust.  ILAUN  =  3  is  a  flag  indicating  that  the 
missile  has  been  launched;  it  is  used  in  this  case  as  the  criterion 
for  the  fighter  switching  to  RTRN2  (constant-Mach-number,  constant- 
altitude  right  turn);  ITHR  »  1  indicates  tables  for  afterburner  thrust. 

The  missile  is  being  held  in  CAPFLT(2,1),  captive  flight  on  the 
fighter,  until  the  range  between  fighter  and  target  (RREL(2,4))  is 
less  than  6500  ft  and  the  bearing  angle  between  fighter  and  target  is 
less  than  30  deg.  The  missile  Is  launched  with  boost  (DELV)  equal  to 
zero,  no  aerodynamics  (IAER0  *  3),  and  zero  thrust  (ITHR  =  5).  The 
MISILX  routine  is  called  to  simulate  the  aerodynamics  and  flight  of  a 
proportional  navigation  missile. 

The  program  automatically  computes  the  closest  miss  distance  if 
the  missile  comes  within  MINMR  (DATA  124)  of  the  target.  This  data 
number  should  be  set  at  a  large  value  (e.g.,  1000  ft).  The  program 
will  terminate  after  finding  the  miss  distance  unless  otherwise 
specified. 


ooo  non  oooo 
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SAMPLr-  POLICY  FOR  LAUNCHING  ONE  MISSILE 


♦♦♦♦  SPECIFY  OUTPUT 
NPR!MT*3 
IPRINTIIMI 
IPRINTI2)*2 
IPRINTn)*3 


**««••««**  FIGHTER  COMMANDS  ♦♦♦♦♦♦♦♦♦♦ 

60  TO  U10«120tl30l«JPaL 
110  CONTINUE 

IF  IILAUN  .EQ.  3)  GO  TO  120  «ILAUN«3  AT  LAUNCH 

CALL  LEAOCLf lf0«0.6500.0«2,21 

GO  TO  190  ^FIGHTER  FLYING  LEAD  COLLISION  NAVIGATIONt 

120  CONTINUE  MILITARY  THRUST*  UNTIL  LAUNCH  -  THEN 

CALL  RTRN2U*2tl)  SWITCHES  TO  CONSTANT  MACH  RIGHT  TURN* 

JPOL*2  AFTER-BURNER  THRUST 

130  CONTINUE 
190  CONTINUE 


«#«*«*****  MISSILE 
GO  TO  (210*220 
210  CONTINUE 

tF(A6S(BEARNG( 
1  GO 

CALL  CAPFLT(2, 
GO  TO  290 
220  CONTINUE 

CALL  LAUNCH! 2* 
230  CONTINUE 

CALL  MISILX(2) 
KPOL«3 

290  CONTINUE 


COMMANDS  ********** 

•230)«KPOL 

2)}  .LT.  30.0«RA0  .AND.  RRELt2*4)  oLT.  6500.01 
TO  220 

1)  ♦MISSILE  LAUNCHED  IF  RELATIVE  RANGE  BETWEEN 

FIGHTER-TARGET  IS  LESS  THAN  6500  FT  AND 
FIGHTER  IS  WITHIN  30  DEG  OFF  TARGETS  TAIL 

0.0«3,5I 


C 

c 

c 


c 

c 


**********  target 

60  TO  1310*320 
310  CONTINUE 

(F  (ILAUN  .EQ. 
THR0TLI3)«.8 
CALL  STRFLT(3f 
60  TO  390 
320  CONTINUE 

CALL  LTRN3I3*5. 5*2*1) 
LP0L*2 

330  CONTINUE 
390  CONTINUE 


RETURN 

END 


COMMANDS  ♦♦♦♦♦♦♦♦♦♦ 

*330)*LPQL 

31  GO  TO  320 

♦THROTL  PROPORTIONS  THRUST 

•  2*2) 


♦TARGET  ON  STRAIGHT  FLIGHT  UNTIL  LAUNCH, 
THEN  PULLS  5.5  6  CONSTANT  MACH  LEFT 
TURN*  AFTER-BURNER  THRUST 


Fig.  18 


Sample  P0LiCY  subroutine  for  launching  one  missile 


Target  connands  are  the  following:  Call  STRFLT(1AER0  *2,  ITHR  * 

2)  until  the  missile  is  launched  (ILAUN  3).  THR0TL(3)  is  set  equal 

to  0<8  to  give  an  80-percent  throttle  setting.  The  target  then  pulls 
LTRN3  (constant-Kach-number  left  turn),  where  the  g's  are  specified  as 
5.5  in  the  argument.  ITHR  »  1  indicates  that  tables  are  to  be  used 
for  obtaining  a  value  for  afterburner  thrust.  See  Appendix  D  for  fur¬ 
ther  instructions  on  the  calling  of  individual  maneuvers. 

Initial  Data 

Since  table  values  are  used  for  the  aerodynamic  and  propulsion 
computations,  data  decks  for  Che  F-104  and  F-105  aircraft  are  entered 
following  Che  JVEH  flag  card,  which  indicates  whici\  vehicles  have  tables. 

Figure  19  shows  the  initial  position  and  velocity  of  the  vehicles 
as  entered  in  the  data.  The  fighter  is  at  the  origin  at  an  altitude  of 
20,000  ft.  The  target  is  5000  ft  down  the  -x-axls,  9000  ft  along  the 
y-axls,  and  at  25,000  ft  altitude.  In  this  case,  velocity  is  read  in 
as  Mach  number  instead  of  as  fc/sec.  To  do  this,  JATM0S  (DATA  20)  is 
set  equal  to  1.  Only  the  magnitude  of  the  fighter’s  velocity  (i.e., 

Mach  1.2)  is  read  in,  and  the  program  automatically  aims  the  fighter 
at  the  target  for  a  lead  collision  course.  The  aiming  routine  is 
triggered  by  setting  IVF  (DATA  16)  equal  to  2.  (See  Appendix  E  for 
details  on  thf  flags  for  reading  position  and  velocity.)  The  target 
is  flying  Mach  0.9,  heading  down  the  y-axis  (V(3,5)  ■  90  deg).  Other 
data  used  in  this  example  are  Che  following: 

Starting  value  for  Integration  step  size  (DT0)  *  0.01. 

Structural  lateral  acceleration  limit  of  aircraft 

(ASMAX)  »  7.3,  8.0. 

Navigation  constant  for  closed-loop  guidance  routines 

(LAMDAO)  -  4.0,  4.0. 

Area  of  aircraft  (AREA)  »  196.0,  385.0. 

Initial  weight  of  aircraft  (WO)  ••  16,699.0,  33,287.0. 

Time  interval  for  printing  output  (DTP0)  ■  0.5. 

Time  value  at  which  program  is  to  stop  (T0T<iL)  ••  25.0. 

Flag  specifying  that  Initial-condition  value  of  velocity 

of  aircraft  is  expressed  as  Mach  number  (JATM0S)  >■ 


ig.  19  —  ThreeKliniensional  diagram  of  initial  conditions  for  Example  2 
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Missile  range  to  target  within  which  program  will 
automatically  Initiate  process  for  miss-distance 
computation  (HINMR)  »  1000.0  ft. 

Since  missile  parameters  such  as  weight  and  area  are  defined  within 
the  missile  subroutine,  their  values  do  net  have  to  be  read  in.  The 
inltjal  position  and  velocity  of  the  missile  are  set  equal  to  that  of 
the  fighter. 

Figure  20  shows  the  initial  data  for  this  example  entered  on  the 
input  form.  Refer  to  Section  X  for  instructions  on  preparing  such  a 
data  package. 

EXAMPLE  3;  GROUMD-IAUNCHED  MISSILE 

A  SAM  is  launched  from  the  ground  at  a  constant-velocity  target. 
Analytic  functions  are  used  to  compute  the  missl'le’s  aerodynamics.  The 
closest  miss  distance  between  missile  and  target  is  to  be  found  and  the 
program  terminated. 

P0LT.CY  Subroutine 

The  following  POLICY  statements  define  the  problem: 

Vehicle  1 

o  Not  used. 

Vehicle  2  (SAM) 

o  Laimch  missile  at  time  zero  with  a  constant  thrust  of 
5397  lb  and  an  initial  velocity  of  100  ft/sec. 

o  Fly  in  0-g  dive  until  burnout  at  1.0  sec. 

o  Continue  0-g  dive  until  time  to  guide;  thrust:  0.0. 

0  Begin  guidance  1.0  sec  after  launch,  fly  proportional 
navlgatlcn  with  a  time  lag  of  0.2  see;  thrust;  0.0. 

Vehicle  3  (Target  Aircraft) 

o  Fly  straight  and  level  at  constant  speed,  no  aerodynamic 
or  propulsion  computations. 


TACTICS  PROGRAM  -  INPUT  FORM' 


f 
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Figure  21  shows  the  P0LICY  subroutine.  The  printout  is  to  ocn- 
sist  of  standard  output,  aerodynamics,  and  attitude  angles.  Becaut-e 
the  fighter  is  not  used  in  this  run,  its  quantities  are  set  to  zero  by 
calling  CAPFLT(1,2). 

The  missile  is  launched  at  time  zero;  DIVEl  with  0  g's  is  called 
(l.e.,  a  ballistic  trajectory).  Analytic  functions  are  used  to  compute 
aerodynamics,  so  IAER0  »  1;  constant  thrust  is  indicated,  so  ITHR  4. 
TLAIJN(2)  is  the  time  of  launch  and  can  be  used  as  a  criterion  in  P0LtCY. 
Here  it  is  used  both  to  determine  burnout  (TBURK(l))  and  the  time  the 
missile  is  to  begin  guiding  (TGIIIDE(2)) .  Missile  guidance  begins  by 
calling  PR0NAV  (proportional  navigation).  Thrust  now  equals  zero,  so 
ITHR  «  5.  A  time  lag  is  added  to  the  missile  at  this  point  by  setting 
ITAU(2)  “  1,  which  indicates  that  only  one  time  lag  is  being  used; 
TAU(2,1)  ■  0,2,  indicating  a  0.2-sec  lag.  STRFLT  is  called  for  the 
target  with  no  aerodynamic  or  thrust  computations,  so  IAER0  «  3  and 
ITHR  »  5.  See  Appendix  D  for  further  Instructions  on  the  calling  of 
individual  maneuvers. 

Initial  Data 

Figure  22  is  a  three-dimensional  diagram  of  the  initial  positions 
and  velocities  of  the  missile  and  target  specified  in  the  input  data. 
The  program  is  capable  of  reading  in  values  for  the  position  and  ve¬ 
locity  of  the  missile,  but  in  this  case  the  initial  conditions  are  set 
to  those  of  the  fighter  instead,  i.e.,  the  missile's  values  are  read 
in  as  the  fighter's  position  and  velocity.  The  fighter  values  are 
later  zeroed  out,  since  they  are  not  used  in  this  problem.  This 
seemingly  needless  complexity  may  be  used  to  advantage  for  automatic 
aiming  of  vehicle  1  against  3  or  vice  versa  using  the  flags  I’/F  *  2  or 
IVT  *  2.  This  will  automatically  call  the  AIM  routine  and,  in  this 
case,  provide  for  initially  aiming  2  against  3  (see  the  description  of 
subroutine  AIM  in  Appendix  C) . 

The  missile  is  initially  placed  at  the  origin,  and  the  altitude 
component  is  accordingly  zero.  The  target  begins  its  flight  3000  ft 
down  the  -x-axis,  10,000  ft  along  the  y-axis,  and  at  a  500-ft  altitude. 


nor*  noo  nr»o  ono 


S«^FL£  PCIECY  FCR  GRCLKL  LAtWEC  FISSILE 

•••  SPECIFY  CUTPLT 
fcFi»IfcT*2 
IPRIliTf  llsl 
IPRIfcTI2J»2 
IPPIX1«3I=3 


FIGFIEP  CCPPAKCS  4444444444 
GC  TC  I11C,12C,13C),JPCL 
lie  CCKIIKLE 

CALL  CAPFLTtl,2}  4FiGhIER  ACT  8EIAG  USED  SC  CUA^TITIES 

GC  TC  ISC  ZEBCEC  CLt 

i2C  CCMIALE 
13C  CCATiAkE 
ISC  CCMIALE 


4444444444  FISSILE  CCPPAKCS  4444444444 


21C 

22C 

23C 

2AC 

25C 

29C 


GC  TC  f.21Ct22C«23C,24C*25C}«KPGL 
CCKTIAUE 

CALL  LALACH(2,C.Ct3«5) 

CCAllALE 

IF  (ITIPE^TLAUKiail  .GT. 

CALL  CIVE1(2«C.C«1,4} 

KPCL«2 
GC  TC  2SC 
CCMIftUE 

if  CITIPE-ll5tM2))  .GT, 

CALL  CIVEl<2«C.C«lt5} 

KPCl*3 
GC  TC  290 
CCKTIACE 
lTAtl2>=l 


♦FISSILE  LALNCKEO  AT  TIFE  2EltC 
TELRMl  GO  TC  23C 

4FISS3LE  FLIES  IK  ZERO  G  DIVE  klfH 
TFRLST»CCKSTANT  CKTIL  BURkCUTITBURKI) 
,AT  liUCh  TIF.E  ThROST  SET  TC  ZERC 
AKC  FISSILE  COKIINLES  UKGUICEC 
TGLICEIZn  GC  TC  2SC 


♦A  TIFE  LAG  CF  .2  SEC  IS  ADDED  TO  FISSILE 
AFTER  IT  BEGINS  GUIDING 


TAL(2tl}=.2 

CCKTIKUE 

CALL  FRCKAVl2tlt5l 

KPCL*5 

CCKTIKUE 


♦kFEK  TIFE  JINCE  LAUNCH  IS  GREATER  THAN 
TIFE  THE  FISSILE  IS  TC  FLY  UNGUiCEC 
(TGLICE(2»}«  PRCPORTICNAL  NAVIGATICK  IS 
CALLED 


44«44444«4  TARGET  CCFFANCS  ♦4444#444# 

C-C  TC  {2IC,32C,23C),IFCL 
31C  CCNTIKUt 

CALL  STRFLT(3,3,5)  ♦lARGfT  FLIES  STRAIGHT  FLIGHT 

GC  TC  39C 
32C  CCMIKLE 
:3C  CCNTINCE 
39C  CCMIKtE 
C 
C 

RETLRK 

ENC 


Fig.  21  —  Sample  P0LICY  subroutine  for  ground-launched  missile 


■ 
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'Hbe  obissilc's  Solciftl  Telccicjr  3s  100  £c/scc.  asod  ic  is  pcisted 
alractlj  <5oiic  uSm  ICS  ^csswea  aiissile  sad  tartsc,  S^viaf  ic  a  fttrsdlpg 
a£  ^B.9  d«s  is  ^ae  Easrizcoral  piaae  (^Cl>5]^)  aaad  2.55  cck  la  the  ver- 
cicsl  plaas  C^Cl,6)).  l&e  cargec  cats  a  Telocity  of  13«0  fc/sec  beaded 
dews  Cae  x-axls  (VCS,£}  >  90  def). 

A  slgaificaac  quaaciCf  of  aexedjfaaBic  data  omsc  be  eatered  for 
tbis  case,  sloce  aaaljcic  ftsaedetes  are  csed  co  coapate  the  arlssile 
aerodjraaaics.  aad  eqaatloa  conscaats  aese  cserefore  be  specified  for 
use  with  Eq.  (26).  the  follovisi  data  used  la  this  ezaaple: 

Starciag  value  for  latesratioa  st^  size  isfit)  *  0.91. 

Tioe  iotezvsl  missile  is  co  fljr  tmguided  afcer  launch 

(TCijID£(2))  »  l.a. 

Firsc-scate  ba^miag  else  of  adsslle  (IKTESil)  «  1.0  sec. 

Stmccural  lateral  acceleraclac  llmic  of  missile  acd 
target  (ASHA£)  «  15.0,  7.3  g*s. 

Ccastaats  to  be  used  for  thrust  (7BC9!8(2)}  *  5397. 

SEavlgatloa  constant  for  closed-loop  guidance  routines 
of  missile  (LAMIttO(2))  «  4.0. 

Maximum  aerodynamic  lift  coefficient  of  missile  (CI2ttX(2))  *  1.4. 

d(CQ)/d(£^),  constant  of  missile  (BC^(2))  ^  0.G041666667. 

dC  /do  (assuced  to  be  a  constant)  SL#FE(2)  ~  0.1. 

^  * 

Specific  impulse  of  rocket  motor  (1MPLS£(2))  »  250.0. 

Initial  vei^t  of  missile  (UO)  50. C. 

Area  of  missile  (AREA)  »  0.264. 

Integration  node  (in  this  case,  fixed-step  Adams-Moulton) 

(JINTEG)  -  2. 

Tine  interval  for  printing  output  (DTPO)  ■9.5. 

Time  value  at  which  program  is  to  step  (T0TAL)  ■  25.0. 

Missile  range  to  target  within  which  program  will  automatically 
initiate  process  for  miss-distance  computation  (MINMR)  ■  1000.0. 

Figure  23  shows  the  initial  data  for  this  example  entered  on  the 
Input  form.  Refer  to  Section  X  for  instructions  on  preparing  such  a 
data  package. 


These  quantities  must  be  set  for  computing  missile  weight  during 
the  propulsion  Interval. 


cuMXf  I)  vfMacus» 
sf7l 


I  t 


Fig.  23  —  Inpuf  data  for  Example  3 


EXAMPLE  4:  lAUilCHINC  TWO  MISSILES,  RECALL 


The  recall  features  of  Che  program  provide  the  capability  of 
launching  more  chan  one  missile.  After  Che  missile  has  been  launched 
the  first  time  and  the  miss  distance  determined,  the  program  continues 
instead  of  terminating,  and  the  missile  can  be  recalled  either  to 
vehicle  1  or  to  vehicle  3.  The  aerial  combat  is  then  continued  and 
the  missile  is  launched  again.  Therefore,  if  vehicle  1  misses  on  the 
first  launching,  it  can  repeat  the  process  by  launching  another  missile, 
or  vehicle  3  can  retaliate  by  launching  a  missile  at  vehicle  1. 

In  this  example,  vehicle  1  is  pursuing  vehicle  3  and  launches  a 
missile  when  the  launch  criteria  have  been  satisfied.  After  the  missile 
has  missed  the  target  it  is  recalled  to  vehicle  3,  and  now  vehicle  1 
is  the  target.  The  missile  is  now  launched  against  vehicle  1. 

PtolCY  Subroutine 

The  following  P0L1CY  statements  define  the  problem: 

Vehicle  1  (P-104  Intercepting  Aircraft) 

o  Fly  pursuit-course  navigation;  thrust,  military  power. 

o  If  the  range  to  target  is  less  than  7000  ft,  launch  missile. 

o  After  launch  perform  constant-Mach-number,  constant-altitude 
left  turn;  thrust,  afterburner. 

o  When  heading  angle  in  horizontal  plane  is  greater  than 
270  deg,  fly  straight  and  level;  thrust,  military. 

Vehicle  2  (First  Missile) 

o  Hold  in  captive  flight  by  fighter  until  launch  criterion 
is  satisfied. 

o  Launch,  boost,  fly  unguided  and  then  guided  in  accordance 
with  guidance  and  aerodynamic  characteristics  specified 
in  special  missile  routine. 

o  When  range  rate  (missile-target)  becomes  greater  than  zero, 
initiate  process  for  finding  miss  distance. 


Vehicle  3  (F-IOS  Target  Aircraft) 

o  Fly  straight  and  level;  thrust,  military. 

o  When  missile  is  launched,  perform  5-g  left  turn;  thrust, 
military. 

o  After  fighter-launched  missile  has  missed,  switch  to  pursuit- 
course  navigation  tracking  vehicle  1;  thrust,  afterburner. 

o  If  range  to  fighter  is  less  than  10,000  ft  and  greater  than 
7000  ft  and  the  target  is  within  30  deg  off  fighter's  tail, 
launch  second  missile,  which  has  been  recalled  to  vehicle  3. 

Vehicle  4  (Second  Missile) 

o  Missile  recalled  to  vehicle  3  after  first  missile  has  missed. 

o  Hold  in  captive  flight  by  target  until  launch  criterion  is 
satisfied. 

o  Launch,  boost,  fly  unguided  and  then  guided  in  accordance  with 
guidance  and  aerodynamic  characteristics  specified  in  special 
missile  subroutine. 

o  When  range  rate  (misslla-flghter)  becomes  greater  than  zero, 
initiate  process  for  finding  miss  distance  and  terminate 
program. 

Figure  24  shows  the  P0LICY  subroutine.  The  printout  is  to  con¬ 
sist  of  star  lard  output,  aerodynamics,  and  attitude  angles. 

Fighter  flies  PRSUIT  (pursuit-course  navigation)  with  table  values 
for  aerodynamics  (IAER0  ■  2)  and  table  value-  for  military  thrust 
(ITHR  *  2).  At  launch  time  (ILAUN  ■  3  Indicates  that  missile  has  been 
launched)  LTRN2  (constant-Mach-number,  constant-altitude  left  turn)  is 
called  with  afterburner  thrust  (ITHR  =  1).  Fighter  flies  LTRN2  until 
its  horizontal  heading  (V(l,5))  is  greater  than  270  deg.  This  is  used 
as  the  criterion  for  switching  to  STRFLT. 

The  first  missile  is  held  in  CAPFLT(2,1),  captive  flight  by 
fighter,  until  the  range  between  fighter  and  target  (RREL(2,4))  is 
less  than  7000  ft.  The  missile  is  launched,  and  the  MISILl  routine 
is  called  to  simulate  the  aerodynamics  and  flight  of  a  proportional 
navigation  missile. 

To  initiate  the  recall  feature  of  the  program,  it  is  necessary  to 
set  the  flag  IMISS  =  1,  which  causes  the  program  to  continue  after 
finding  the  first  missile  miss.  It  is  reset  to  zero  for  the  second 
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C  SAMPLE  PCLICY  FCR  L4LhCHK6  TkC  PISSILES 
C 

c 


C*'  <  SPECIFY  CLTPLT 
KPRINT«3 
IPRIMIl)*! 

IPRIKTI2»»2 

IPRIM<3)»3 

C 

C 

C  ♦♦♦♦♦♦♦#♦♦  FIGETER  CCPHANCS  ♦♦♦♦♦♦♦♦♦♦ 


12C,t3C)*JPG: 


GC  TC  (110, 
lie  CCMHLt 

IF  (ILALN  .EC.  3)  GC 
CALL  PRSL1T(1,2,2) 

GC  TC  ISC 
12C  CCKTIKL'E 

IF  (V(l,5)  .GT.  2T0.0«RAC) 
CALL  LTRN2(1,2,1) 

JPCL*2 
GC  TC  ISC 
13C  CCNTIhUE 

CALL  STRFLT(l,2,2l 
JPCL*3 

ISC  CCNTIKUE 


TC  12C 
♦FIGFTER 
THBLST, 


♦lLAtN*3  AT  LALKCH 
FLIES  PRSUIT  COURSE  NAVIGATION, 
LNTIL  MISSILE  LALNChEC 


MILITARY 


GC  TC  13C 
♦  AT  LALNCF  FIGHER 
CCNSTANT  ALTITUCE 


SRITCF.ES  TC  CCNSTANT  MACH, 

LEFT  TURN,  AFTER-BURNER  THRUST 


C 
C 

C  ♦♦*♦♦♦♦♦♦♦  MISSILE  CCMMANCS  ♦♦♦♦♦♦♦♦♦♦ 


♦NFEN  HEACING  ANGLE  IN  HCRIZCNTAL  PLANE  (THETA) 
IS  GREATER  THAN  270  CEG,  FIGHTER  FLIES  STRAIGHT 
FLIGHT,  MILITARY  THRUST 


GC  TC  (210,22C,23C),KPCL 
210  CCNTINUE 

IF  (KMISS  .EC.  1)  GC  TC  2SC 

IMISS»l  ♦FLAG  FCR  PROGRAM  TO  CCNTINUE  AFTER  MISSILE  MISS 

IF  (RREL(2,A)  .LE.  7CCC.C)  GC  TC  220 

CALL  CAFFLT(2,n  AIF  RELATIVE  RANGE  BETWEEN  FIGHTER-TARGET  IS  LESS 

GC  TC  2SC  THAN  7CCC  FT.,  LAUNCH.  MISSILE.  CTHERNISE  MISSILE 

220  CCNTINUE  HELC  IN  CAPTIVE  FLIGHT  BY  FIGHTER 

CALL  LALNCH(2,C.C,3,5) 

23C  CCNTINUE 

IF  IIMISS  .EC.  2)  GC  TO  2AC 

CALL  MISILi{2)  AIF  IMISS*2  PROGRAM  IS  READY  TC  CCNTINUE  AFTER 

KPCL=3  FINCING  MINIMUM  MISS  DISTANCE  CF  FIRST  MISSILE 

GC  TC  290 
2AC  CONTINUE 

KMISS=1  ♦RMISS-l  INDICATES  SECCND  MISSILE  IS  TC  BE  USED 

KPCL=l 

29C  CCNTINUE 


C 

C 

Q  target  ccmmancs  ♦♦♦♦♦♦♦♦♦♦ 


Continued 


Fig.  24  —  Sample  P0LICY  subroutine  for  launching  two  missiles 
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39C 

C 

C 

C  *** 


GC  TC  (J10,32C, 

CCKTINL'E 

IF  (ILAL^  ,EC. 

CALL  STRFLT(3,2 

GC  TC  39C 

CC^TI^UE 

IF  (Kf'ISS  .EC. 

CALL  LTRM(3,5. 

LPCL=2 

GC  TC  3SC 

CCNTINUE 

CALL  PRSUIT(3,2 

LPCL=3 

CENTIME 


33C),LPCL 
3)  GC  TO  32C 

,2)  ♦TARGET  FLIES  STRAIGHT  FLIGHT,  MILITARY  THRUST 


1)  GC  TC  33C 

C,2,2)  ♦AFTER  LAUNCH,  TARGET  PULLS  5  G  LEFT  TURK* 
MILITARY  THRLST 


♦AFTER  FIGHTER  LAUNChEC  FISSILE  HAS  FISSEC,  TARGET 
TARGET  ShITCHES  TC  PLRSIIT  CCURSE  NAVIGATION, 
AFTER-ELRNER  THRLST 


*#*#«i»*  SECCNC  F 
GC  TC  (AI0,42C, 
CCNTINC'E 
IF  (KFISS  .NE. 
IFISS»C 
CONTINUE 
IF  (PRELt2,4{  . 
lABS(eEARNG(2}  . 
CALL  CAPFLT(2,3 
FPCL*2 
GC  TC  A9C 
CCNTINUE 
CALL  LALNCH(2,C 
CCNTINUE 
CALL  FISIL1(2) 
FPCL=A 
CCNTINUE 


ISSUE  CCFFANCS  ♦♦♦♦♦♦♦♦♦♦ 

43C,A4CJ ,FPCL 

U  GC  TO  ASC 

♦FLAG  FCR  PRCGRAF  TO  STOP  AFTER  FINDING  FINIFUN 
FISS  DISTANCE  FCR  SECCNC  FISSILE 
LT.  ICCCC.C  .ANC.  RREL(2tA)  .GT.  7CCC.0  .AND. 

GT.  ISC.C^RACl  GC  TC  430 

)  ♦FISSILE  LAUNCHED  IF  RELATIVE  RANGE  BETWEEN  TARGET 

-FIGHTER  IS  GREATER  THAN  70CC  FT.  AND  LESS  THAN 
1C,CCC  FI.,  AND  THE  TARGET  IS  WITHIN  30  DEG  CFF 
FIGHTERS  TAIL.  OTHERWISE  FISSILE  IS  HELD  IN 
.0,3,5)  CAPTIVE  FLIGHT  BY  TARGET 


RETURN 

END 


Fig.  24  (continued) 
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misslle,  and  the  program  will  terminate  after  the  miss-distance  calcu¬ 
lation.  IMISS  ■  2  is  a  flag  set  in  the  program  to  indicate  that  the 
miss  distance  for  the  first  missile  launch  has  been  computed  and  that 
the  program  is  ready  to  continue.  This  flag  is  used  as  the  criterion 
in  P0LICy  for  triggering  the  flag  KMISS  ^  1,  which  indicates  that  the 
first  missile  is  no  longer  in  flight  and  that  the  second  missile’s 
cennands  should  now  be  followed.  (See  the  first  and  second  missile 
command  sections  in  F0LICY,  Fig.  24.) 

The  second  missile  is  held  in  captive  flight  by  target  CAPFLT(2,3) 
until  the  target-fighter  range  (RREL(2,4))  is  greater  than  7000  ft  but 
less  than  10,000  ft  and  the  fighter-target  bearing  angle  is  greater 
than  150  deg.  (See  Section  XI  for  che  definition  of  bearing  angle.) 

The  missile  is  launched  against  the  fighter  and  MISILl  routine  called. 
Flag  IMISS  =  0  Indicates  that  the  program  will  stop  after  determining 
miss  distance. 

Target  flies  STRFLT  (IAER0  =•  2,  ITHR  =  2)  until  the  first  missile 
is  launched  (ILAUN  *3).  It  then  pulls  LTRNl  (simple  left  turn)  with 
5  g's  specified  until  the  first  missile  has  missed  and  is  recalled  to 
the  target  (this  is  indicated  by  KMISS  ■  1) .  Target  now  flies  PRSUIT 
against  the  fighter  attempting  to  arrive  within  firing  range  at  after¬ 
burner  thrust,  ITHR  «  1.  See  Appendix  D  for  instructions  on  calling 
the  optional  subroutines  used  in  F0LICY. 

Initial  Data 

Since  table  values  are  used  for  computing  aerodynamics,  data 
d^cks  for  the  F-104  and  F-105  are  entered  following  the  JVEH  flag 
card,  as  shown  in  Fig.  11. 

Figure  25  is  a  three-dimensional  diagram  of  the  initial  positions 
and  velocities  used  in  this  example.  As  can  be  seen,  the  fighter  is 
6000  ft  along  the  y-axis  at  15,000  ft  altitude.  It  has  a  velocity  of 
Mach  0.92  (velocity  is  given  in  Mach  number  instead  of  ft/sec  if  JATM0S 
(DATA  20)  is  set  eqvial  to  1)  and  a  horizontal  heading  of  210  deg 
(V(l,5)).  The  target's  initial  position  is  8000  ft  down  the  -x-axis 
and  at  15,000  ft  altitude.  Initial  velocity  is  Mach  0.87  with  180-deg 
heading  angle  (V(3,5)).  Other  data  needed  for  this  example  are  the 
following: 
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Startlng  value  for  Integration  step  size  (DT0)  =»  0.01. 

Structural  lateral  acceleration  limit  of  aircraft 
(ASMAX)  «■  6.8,  7.0. 

Navigation  constants  for  closed-loop  guidance  routines 
(LAMDAO)  -  40.0,  40.0. 

Area  of  aircraft  (AREA)  =  196.0,  385.0. 

Initial  weight  of  aircraft  (WO)  »  16,699.0,  33,283.0. 

Time  interval  for  printing  output  (DTP0)  *  0.5. 

Time  value  at  which  program  is  to  stop  (T0TAL)  =  60.0. 

Flag  specifying  that  initial-condition  value  of  velocity 
of  aircraft  is  expressed  as  Mach  number  (JATM0S)  =*  1. 

Missile  range  to  target  within  which  program  will  auto¬ 
matically  initiate  process  for  miss-distance  computation 
(MINMR)  »  1000.0. 

Because  missile  constants  and  variables  are  defined  within  the 
missile  subroutines,  their  values  do  not  have  to  be  read  in.  The 
initial  position  and  velocity  of  the  missile  are  set  equal  to  those 
of  the  fighter. 

Figure  26  shows  the  initial  data  for  this  example  as  entered  on 
the  input  form.  Refer  to  Section  X  for  instructions  for  preparing  such 
a  data  package. 

EXAMPLE  5;  USING  RESTORE  FEATURE  TO  FIRE  20-MM  PROJECTILES 

After  vehicle  2  is  launched  and  a  hit  or  miss  has  occurred,  this 
option  enables  the  program  to  restore  all  numerical  values  existing  at 
launch  time.  Now  events  may  proceed,  a  new  launch  may  take  place,  and 
characteristics  and  parameters  may  change.  In  this  case  the  restore 
feature  is  used  to  fire  consecutive  projectiles  at  O.Ol-sec  intervals 
at  a  ground  target.  A  projectile  is  launched  first  at  a  given  range 
and  the  miss  distance  computed.  The  numerical  values  are  then  restored 
to  those  existing  at  launch  time.  Integration  then  occurs,  and  the 
next  projectile  may  be  launched  at  a  subsequent  time. 

P0LICY  Subroutine 


The  following  P0LICY  statements  define  the  problem: 


Fig.  26  —  Input  data  for  Example  4 
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Vikicle  i  {F-iJ-4  Purauinsj  Aivanift) 

o  Fly  pursuit-course  navigation  against  a  stationary 
ground  target;  thrust,  ssilitary. 

o  Uhen  a  ground  range  of  3000  ft  has  been  covered,  begin 
launching  20-ca  projectiles  at  0.01-sec  Intervals. 

o  If  altitude  falls  below  3000  ft,  perfons  a  4-g  climb; 
thn»t,  afterburner. 

Vehicle  2  (20-m  Projectile) 

o  Fly  captive  flight  until  launch  criterion  is  satisfied. 

o  Launch  with  an  incremental  velocity  of  &V  and  fly  in  accor¬ 
dance  with  characteristics  specified  in  B20MH  subroutine, 
i.e.,  a  ballistic  trajectory. 

o  Find  miss  distance  at  point  of  ground  impact,  restore  program 
to  values  cf  launch  time,  launch  second  projectile  0.01  sec 
later,  and  continue  until  five  projectiles  have  been  fired- 

Vehicle  3  (Stationary  Ground  Target) 

o  Re^in  stationary  on  ground. 

Figure  27  shows  the  PjiLICY  subroutine.  The  printout  is  to  con- 
cist  of  standard  output,  aerodynamics,  and  attitude  angles.  To  use 
the  restore  feature  of  the  program,  flag  IST0RE  =  1  must  be  set.  Tlag 
IMISS  =  1  must  also  be  specified  so  that  the  picgram  will  not  stop  after 
computing  the  missile  miss  distance.  Both  these  flags  are  set.  the  first 
time  thrcjgh  P0LICY  and  must  be  reset  for  each  separate  launch. 

The  fighter  is  flying  PRSl'IT  (pursuit-course  navigation)  with  table 
values  for  aerodynamics  (IAER0  =  2)  and  table  values  for  military  thrust 
(ITHR  =2).  If  its  altitude  falls  below  5000  ft,  the  fighter  begins  a 
4-g  climb  (CLlMBl)  with  afterburner  thrust  (ITHR  =  1)  in  order  to  pull 
out.  The  projectiles  are  held  in  captive  flight  by  the  fighter 
(CAPFLT(2,1)) ,  until  ;he  fighter  hac  ^lown  3000  ft  along  the  y-axis; 
i.e.,  the  y-component  (R(l,2))  is  greater  than  3000.  Each  projectile 
is  launched  with  a  boost  velocity  (DELV)  of  2800  ft/sec,  and  the  B20MM 
subroutine  is  called  to  determine  the  aerodynamics  and  flight  path. 
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C  SAMPLE  POLICY  USING  STORE  TG  FIRE  20MM  CANNON  PROJECTILES 

C 

C****  SPECIFY  OUTPUT 
NPRINT*3 
IPRINT(i)«l 
1PRINT421«2 
IPRINT(3)c3 

C 

C  •♦♦♦♦♦♦♦♦♦  FIGHTER  COMMANDS  ♦♦♦♦♦♦♦♦♦♦ 

GO  TO  (110*120«130)»JPQL 
110  CONTINUE 

IMISS*1  ♦FLAG  FOR  PROGRAM  TO  CONTINUE  AFTER  MISSILE  MISS 

IST0RE>1  ♦FLAG  FOR  STORING  VALUES  AT  LAUNCH 

120  CONTINUE 

IF  (Rfl»3l  .LE.  5000.0)  GO  TO  130 
CALL  PRSUITI 1,2,2)  ♦FIGHTER  ON 
JP0L«2 
GO  TO  190 
130  CONTINUE 

CALL  CLIMBK  1,A.0,2,  n 
140  CONTINUE 
190  CONTINUE 

C  #«****«#*«  MISSILE  COMMANDS  ♦•♦♦♦♦♦♦♦♦ 

GO  TO  (210, 220*230, 240, 250),KPOL 
210  CONTINUE 

IF  (R(l,2)  .GE.  3000.0)  GO  TO  220 
CALL  CAPFLT(2,1) 

GO  TO  290 
220  CONTINUE 

CALL  LAUNCH(2»2SOO.O,3,5) 

230  CONTINUE 

IF  (ISTORE  .EG.  0)  GO  TO  240 
CALL  820MM<2) 

KP0L*3 
60  TO  290 
240  CONTINUE 
JINTEG-2 
IMISS-1 

IF  (II  .EQ.  3)  IMISS^cQ 
IST0R£»1 
ICP0I.*5 
KK*0 

250  CONTINUE 

IF  (KK  .NE.  0)  GO  TO  220 
CALL  CAPFLT(2,1) 

KK='i 

ii»n*i 

GO  TO  290 
260  CONTINUE 
290  CONTINUE 

Qt*************************J^RQiT  COMMANDS  ♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

GO  TO  (310, 320, 330, 340), LPOL 
310  CONTINUE 

GO  TO  390  ♦USED  AS  STATIONARY  GROUND  TARGET 

390  CONTINUE 
RETURN 
END 


♦FIXED  STEP  INTEGRATION  FOR  PROJECTILES 
♦IMISS  RESET  FOR  EACH  LAUNCHING 
♦IMISS*0,  STOP  AFTER  5TH  PROJECTILE 
♦ISTORE  RESET  FOR  EACH  LAUNCHING 


♦AFTER  RESTORING  PROGRAM  TO  VALUES  AT 
LAUNCH  TIME,  INTEGRATE  AHEAD  .01  SEC 
BEFORE  LAUNCHING  NEXT  PROJECTILE 


PRSUir  COURSE  NAVIGATION,  TABLE 
TABLE  VALUES  FOR  AERODYNAMICS, 
MILITARY  THRUST 

♦IF  ALTITUDE  IS  LESS  THAN  5000 
FT,  SWITCH  TO  4  G  CLIMB  TO 
PULL  OUT,  AFTER-BURNER  THRUST 


♦IF  FIGHTER  HAS  FLOWN  3000  FT 
ALONG  Y-AXIS,  LAUNCH  PROJECTILE 
WITH  0ELV*2800.  OTHERWISE  HELD 
IN  CAPTIVE  FLIGHT  BY  FIGHTER 


♦ISTORE  SET  BACK  TO  ZERO  WHEN 
PROGRAM  HAS  RESTORED  ITSELF  TO 
LAUNCH  VALUES  AFTER  COMPUTING 
THE  PREVIOUS  PROJECTILE  MISS 


I  /: 

I 


4  I 


Fig.  27  —  Sample  P0LICY  subroutine  for  firing  20-mm  piojectiles 


i 


At  ground  impact,  the  integration  routine  backs  up  to  obtain  the 

it 

distance  between  the  impact  point  and  the  target.  Since  IMISS  >1, 
the  program  then  continues  instead  of  ending,  and  the  values  are  re¬ 
stored  to  those  of  launch  time.  At  this  point  the  flag  IST0RE  >  0 
Indicates  the  program  has  been  restored.  This  is  used  as  the  criterion 
for  switching  back  to  CAPFLT  so  that  the  program  can  integrate  ahead 
0.01  sec  before  launching  again.  This  process  continues  until  five 
projectiles  have  been  launched.  After  the  fifth  launching,  IMISS  =  0, 
and  the  program  ends  after  computing  closest  miss.  Vehicle  3  is  a 
stationary  ground  target  and  therefore  receives  no  velocity  or  accel¬ 
eration  commands.  See  Appendix  D  for  further  instructions  on  the 
calling  of  individual  maneuvers. 

Initial  Data 

Since  table  values  are  used  for  computing  the  aerodynamics  of  the 
fighter,  the  F-104  data  deck  is  entered  immediately  following  the  JVEH 
flag  card,  as  shown  in  Fig.  11. 

Figure  28  is  a  three-dimensional  diagram  describing  the  initial 
positions  and  velocities  of  the  vehicles  in  this  example.  The  fighter 
is  located  at  the  origin  at  an  altitude  of  7000  ft.  It  is  headed 
directly  down  the  y-axis  (V(l,5)  *  90  deg)  with  a  velocity  of  1000 
ft/sec.  The  target  is  located  on  the  ground  6000  ft  along  the  y-axls. 
Since  it  is  stationary,  it  has  no  initial  velocity.  Other  data  needed 
for  this  example  are  the  following: 

Starting  value  for  integration  step  size  (DT0)  =  0.01. 

Structural  lateral  acceleration  limit  of  fighter  (ASMAX)  ■  7.3. 

Navigation  constant  for  closed-loop  guidance  routine 

(LAMDAO)  -  40.0. 

Area  of  fighter  (ARFA)  *  196.1. 

Initial  weight  of  fighter  (WO)  *  19,470.0. 

A 

This  is  not  the  same  as  the  point  of  closest  approach  or  miss 
distance  as  usually  defined.  (See  Section  XIII.)  TACTICS  automat¬ 
ically  makes  the  distinction  if  the  altitude  of  the  target  R(3,3)  is 
exactly  zero. 


Aiming  error  in  launching  projectiles  (DVPHI)  »  3.4. 

Time  interval  for  printing  output  (DTF0)  ■  0.5. 

Time  va.  ue  at  which  program  is  to  stop  (T0TAL)  *  20.0. 

Missile  range  to  target  within  which  program  will  automatically 

initiate  process  for  ground  impact  miss  computation  (MINMR)  == 

1000.0. 

Because  ballistic  constants  and  variables  are  defined  within  the 
B20MM  subroutine,  values  do  not  have  to  be  read  in.  The  initial  posi¬ 
tion  and  velocity  of  the  projectiles  are  set  equal  to  the  fighter 
values . 

Figure  29  shows  the  initial-condition  data  for  this  example  as 
entered  on  the  input  form.  Refer  to  Section  X  for  instructions  on 
preparing  such  a  data  package. 


XIII.  NUMERICAL  INTEGRATION 


MODES  OF  OPERATION 

TACTICS  uses  an  integration  subroutine  that  is  an  updated  and 

modified  version  of  Adams-Moulton,  Runge-Kutta  SHARE  subroutine 

(2) 

RU  INT.  In  addition  to  certain  other  advantages,  this  updated 
version  is  written  in  FORTRAN  IV  source  language,  whereas  RW  INT  is 
written  in  a  machine-oriented  language.  An  excellent  simplified  de¬ 
scription  of  the  Adams-Moulton  predictor-corrector  method  of  integra¬ 
tion  is  given  in  Ref.  1,  and  more  detailed  information  is  given  in 
Refs.  2  and  3. 

The  integration  subroutine  has  four  optional  modes  of  operation: 

o  Adams-Moulton  predictor -corrector  method  using  a  variable  step 
size. 

o  Runge-Kutta  classical  fourth-order  method  using  a  fixed  step 
size. 

o  Adams-Moulton  predictor-corrector  method  using  a  fixed  step  size, 
o  Adams-Moulton  predictor-corrector  method  using  a  variable  step 
size,  controlled  to  allow  printout  exactly  at  specified  intervals. 

The  optional  modes  may  be  selected  by  setting  the  flag  JINTEG  equal  to 

"ft 

0,  1,  2,  or  3,  respectively,  either  in  initial-condition  data  (DATA  122) 
or  in  P0LICY.  If  no  value  is  set,  JINTEG  will  automatically  be  zero 
and  the  variable-step-size  mode  will  be  used  by  the  program. 

Each  process  has  relative  advantages.  The  Adams-Moulton  method 
(JINTEG  =  0)  generally  requires  the  least  execution  time,  since  the 
step  size  is  automatically  controlled  (doubled  or  halved)  to  maintain 
a  value  consistent  with  a  preset  truncation  error  test  limit.  A  pos¬ 
sible  disadvantage  is  that  substantial  overshoots  may  occur,  so  that 
the  printout  time  only  approximates  the  specified  time  (e.g.,  1.22, 

2.34,  3.56  sec,  instead  of  1.0,  2.0,  3.0  sec,  respectively).  This 
disadvantage  has  been  removed  with  only  a  slight  cost  in  execution 

■k 

See  Appendix  E. 
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time  by  modifying  the  variable-step  mode  of  operation  so  that  printout 
will  occur  as  specified  by  the  user.  The  fixed-step  Adams-Moulton  and 
Runge-Kutta  methods  generally  require  the  most  execution  time,  but 
there  are  certain  computational  advantages  in  having  a  predictable  step 
size.  Experience  with  trial  problems  indicates  that  the  advantages  of 
both  variable-  and  fixed-step  methods  may  be  exploited  by  changing  from 
one  to  the  other  during  different  phases  of  the  flight  simulation,  par¬ 
ticularly  just  before  and  after  missile  launch.  This  subject  will  be 
discussed  in  further  detail  in  the  following  subsections. 

INTEGRATION  ACCURACY 

Double  precision  procedures  are  used  to  reduce  round-off  errors 
in  accumulating  variables.  Truncation  errors  are  a  function  of  the 
step  olze  used.  A  significant  feature  of  the  variable-step  Adams- 
Moulton  mode  of  operation  is  that  the  step  size  is  controlled  by  the 
integration  subroutine  so  that  it  will  be  less  than  a  preset  trunca¬ 
tion  error  test  limit.  This  limit,  designated  ERTEST,  is  the  maximum 
allowable  relative  error  in  any  one  of  the  dependent  variables  in  a 
local  region,  as  distinguished  from  errors  accumulating  with  each  step. 
Practical  values  for  ERTEST  range  from  10  to  10  ;  this  value  is 

automatically  set  to  10  ^  unless  otherwise  specified  by  input  data  or 
by  P0LICY.  With  reference  to  the  input  form  of  Fig.  9,  the  value  for 
ERTEST  is  set  by  specifying  the  number  of  significant  digits  of  accu¬ 
racy  required  in  data  location  123  (for  example,  if  7  is  entered,  ERTEST 
will  be  set  at  10  .  It  should  be  emphasized  that  the  truncation  error 

test  limit  is  only  applicable  when  either  of  tne  two  variable-step  inodes 
of  operation  is  used  (JINTEG  ■  0  or  JINTEG  »  3) . 

FIXED-STEP  APPROACH  TO  LAUNCH 

As  previously  indicated,  the  use  of  the  variable-step  options  will 
generally  reduce  execution  time,  which  will  result  in  an  appropriate 
step  size  for  a  specified  truncation  error.  Accordingly,  these  modes 
should  be  used  whenever  the  problem  permits — e.g.,  in  time-consuming 
aerial  acrobatics.  However,  a  disadvantage  exists  due  to  the  inability 
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to  predict  which  step  will  be  used  at  a  particular  phase  of  the  problem, 
particularly  when  missile-launch  criteria  are  about  to  be  satisfied. 

A  problem  arises  because  the  step  may  be  large,  allowing  an  undesirable 
overshoot  of  the  criteria  (similar  to  the  inexactitude  of  printout  in 
the  JINTEG  »  0  mode  of  operation) .  A  remedy  has  been  provided  which 
automatically  causes  launch  criteria  (time,  range  to  target,  etc.)  to 
be  approached  with  a  fixed  minimum  step  size.  Actually,  the  procedure 
is  to  approach  launch  on  a  variable  (or  fixed)  step  and,  when  overshoot 
occurs,  back  up  to  the  previous  time  step,  restore  all  conditions,  and 
approach  again  using  small  fixed  steps.  Moreover,  these  fixed  integra¬ 
tion  steps  will  be  maintained  during  the  missile  flight  time  until  the 
point  of  closest  approach  to  the  target  or  miss  distance  has  been  deter¬ 
mined.  The  previous  integration  mode  of  operation  (whether  fixed  or 
variable)  will  then  be  restored  unless  the  option  is  used  to  terminate 
the  problem  run.  If  the  mode  of  cperr’’ion  was  fixed-step  in  the  first 
place,  back-up  and  fixed-step  approaches  to  launch  would  not  occur. 

MISS-DISTANCE  COMPUTATION 

The  point  of  closest  approach  is  taken  to  be  the  value  of  the 
relative  missile  range  to  target  when  its  time  derivative  equals  zero. 
That  is,  when  approaching  the  target  the  range  rate  will  be  negative, 
when  closest  to  the  target  it  will  be  zero,  and  after  passing  the  tar¬ 
get  it  will  be  positive.  This  criterion  is  used  to  determine  miss  dis¬ 
tance  by  "backing  up"  the  integration  process  in  a  manner  similar  to 
that  described  for  fixed-step  approach  to  launch.  The  procedure  is  to 
return  to  the  previous  time  step  whenever  the  missile-to-target  range 
rate  becomes  positive,  restore  all  conditions,  shrink  the  time  step  by 
a  factor  of  ten,  and  repeat  this  process  until  the  time  step  becomes 
smaller  than  DTMIN  (DATA  134).  If  this  value  is  not  read  in  by  the 
user,  It  is  automatically  set  to  10  .  In  order  to  allow  positive  range 

rates  to  occur  without  initiating  the  process  for  miss-distance  compu¬ 
tation,  an  adoLtloual  condition  is  imposed:  The  missile  range  to  tar¬ 
get  must  be  less  than  MINMR  (DATA  L24) .  Otherwise,  miss  distance  will 
not  be  calculated.  For  ASM  or  SSH  applications  it  is  usually  necessary 
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to  determine  the  ground  range  to  target  at  ground  Impact  (i.e.,  zero 
missile  altitude),  which  (as  mentioned  earlier)  is  generally  not  the 
same  as  the  point  of  closest  approach.  TACTICS  will  automatically 
compute  this  ground  range  by  using  the  altitude  of  vehicle  2  rather 
than  range  rate  as  the  criterion  for  backing  up  to  determine  the  time 
and  range  at  ground  impact.  This  feature  is  triggered  if  the  altitude 
of  vehicle  3  is  exactly  0.0. 

INITIALIZATION 

The  Adams-Moulton  predictor-corrector  modes  of  integration,  both 
fixed-  and  variable-step,  use  previous  values  of  the  variables  to  pre¬ 
dict  and  to  correct  values  for  each  new  step.  In  other  words,  past 
history  is  used  to  predict  the  future.  This  procedure  is  not  appro¬ 
priate  in  many  cases  because  of  discontinuities  or  large  step  changes 
(e.g.,  a  large  instantaneous  reversal  in  acceleration  from  right  to 
left  turn,  or  from  climb  to  dive).  Large  step  changes  can  not  only 
cause  a  significant  waste  of  execution  time  (particularly  when  the 
variable-step  modes  are  used)  but  can  also  cause  important  errors. 

This  possibility  may  be  eliminated  by  reinitializing  the  integration 
process  whenever  interfaces  occur,  e.g.,  changes  in  P0LICY  or  propul¬ 
sion.  Initializing  may  be  thought  of  as  simply  stopping  and  starting 
over  again,  not  using  previous  values  for  prediction  of  new  values. 

An  automatic  reinitializing  feature  has  been  built  into  the  framework 
of  the  P0LICY  subroutine  using  the  integer  variobles  JP0L,  KP0L,  LP0L, 
I'1P0L,  and  NP0L.  A  change  in  any  one  of  these  will  trigger  this  process. 
Accordingly,  in  formulating  a  P0LIcy  subroutine,  the  user  should  be 
careful  to  separate  drastic  changes  in  acceleration,  e.g.,  thrust, 
boost,  or  maneuver  changes,  by  using  these  FORTRAN  variables.  (See 
the  examples  in  Section  XII.) 


< 


-99- 


Appendlx  A 

AIRCRAFT  ATTITUDE  AND  ORIENTATION  ANGLES 


The  relationships  given  in  this  appendix,  contained  in  subroutine 
ATITUD,  are  applicable  to  conventional  airframes  where  coordinated 
turns  are  assumed,  as  defined  in  Section  IV  (see  Fig.  7).  Thus  the 
use  of  this  version  of  ATITUD  is  optional,  depending  on  a  particular 
vehicle's  airframe  characteristics.  For  example,  if  a  cruciform  mis¬ 
sile  airframe  is  to  be  represented,  the  attitude-angle  relationships 
corresponding  to  those  in  this  appendix  may  be  included  in  the  missile 
control-law  routine  without  calling  upon  ATITUD,  or  a  different  version 
of  the  routine  may  be  formulated. 

A  right-hand  set  of  orthogonal  roll,  pitch,  and  yaw  orientation 
axes  is  assumed  as  shown  below. 


where  Ip  is  oriented  in  the  direction  of  the  right  wing,  1^  is  oriented 
in  the  direction  of  the  thrust  vector,  and  forms  a  right-hand  orthog¬ 
onal  system  with  Ip  and  1„  (approximately  in  the  opposite  direction  of 
the  lift  vector  L,  i.e.,  neglecting  the  angle  of  attack  a). 

YAW  AXIS 

With  the  assumption  that  the  vehicle  velocity  vector  is  in  the 
plane  formed  by  1^,  vectors  (i.e.,  no  sideslip  angle),  the  unit 
vector  1^,  may  be  determined  by  examination  of  the  sketch  on  the  fol¬ 
lowing  page; 


PRECEDING  PAGE  BLANK 


i-.  *  1.^  sin  5-1^  CCS  o 


«her«,  frcra  Fig.  8, 


-  B  *«Jl  +  <*oV  +  S  “•  It)  Id 

2_  »  L/L  *  —  - = - 

T-  g 


ROLL  AXIS 


SiailarlTp  the  roll  axis  1_  is  defined  by 


T  *  "S? 


cos  o  +  sin  o 


PITCH  AXIS 


Kith  and  deterained.  the  unit  vector  1^  aiay  b?  readily 
found  froB  the  right-t  id  orthogonal  relationship 


AZl»J7d  A2£?  ELEVATION  ANGLES 

The  unit  vectors  1^,  1^,  and  1^  define  the  attitude  or  orients  ..ion 
if  te  vehicle.  Kov  ws  aay  conveaien*:ly  perform  a  coordinate  transfonca- 
tion  to  resolve  the  relative  range  vector  r  (i  -  1,2,3;  J  ■  1,2,3) 
icto  this  fraoe  of  reference.  The  prinary  objective  is  to  detersiine 
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che  orlent£cion  o£  r(l,j),  the  LOS,  in  terae  of  an  aziauth  angle  n  and 
an  elevation  angle  St  defined  as  shown  in  Fig,  30. 


Fig,  30 — Azimuth  end  elevetien  coordination  tronsformotion 

The  negative  sign  for  ^  is  an  unfortunate  result  of  the  conven¬ 
tional  aarodynaaic  orientation  of  the  pitch,  roll,  and  yaw  vectors. 

(We  wish  4^  to  correspond  to  the  up  direction.) 

The  relative  range  vector  r^^  aay  be  transformed  into  tliis  systea 
by  the  follovlcg  equation: 

^ij  *  ^^ij  * 

where 

*  Ij.  -  r^  ^  cos  S  cos  n 

Cy  •  Ip  -  r^^  cos  S  sin  n  (49) 

'ij  ■  \  ■  'ij  ^ 

Froa  the  above,  we  cay  conclude  that 
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where 


u  * 


-1 

tan 


1 

r 


1 

~P 


1 


r 


(50) 


BAKK  AKCLE 

The  bank  angle  ^  is  defined  with  respect  to  a  wind  axis  or 

s 

velocity-vector  reference  systea  (see  Figs.  7  and  8  and  Eq.  (19)). 
TACTICS  calculates  the  angle  froa  the  expression 


-  tan 


-1  ^  ^ 


^  • 


(51) 


where  Ij^  is  the.  unit  vector  along  the  iift  vector  L  (norcal  to  V). 


ROLL  Ah'GLE 

.loll  angle  ^  (ROLL(I))  is  arbitrarily  defined  as  the  angle  existing 
between  the  pitch  axis,  l.e.,  the  winge  of  the  vehicle,  and  a  line 
through  the  c.g.  of  the  vehicle  both  noraal  to  the  longitudinal  (i.e., 
roll)  axis  and  parallel  to  the  horizontal  plane.  The  detemination  is 
accoaplished  as  follows.  First,  the  spherical  coordinates  6^  and  <Pj, 
of  the  unit  vector  1^  are  computed: 


-  ^  ^ 


(52) 
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Kexc,  the  unit  vectors  and  analogous  to  the  vectors  1^  and 
are  deterained.  That  Is,  1^  is  noraal  to  1^.  and  In  the  horizontal 
plane,  and  1^^  Is  nocaal  to  1^  and  In  a  vertical  plane  so  as  to  fora 
a  right-handed  orthogonal  systea.  The  coaponente  of  these  vectors  are 


1.,  -  -  sin  8_ 

Alx  T 

^1.  -  -  9r  “•  ®i 

1.-1  “  CO*  S.T. 

Aly  T 

^ly  *  -  sin  %  ‘i”  «T 

(53) 

1.1  *  0 

Alz 

The  roll  angle  p  aay  now 

Fig.  31  as  viewed  frcs  the  tip 

be  computed  frc«  the  geosetry  shown  in 

of  unit  vector  1_  (directed  out  of  the 

page) . 


Fig.  31  —  RoU-ongle  geometry 


t 


-104- 


tan 


-1  h  ' 

h  * 


(54) 


As  defined  above,  a  turn  to  the  right  results  in  a  positive  roll  angle 
varying  fros  0  to  ISO  deg,  and  a  turn  to  the  left  results  in  a  negative 
roll  angle  varying  from  0  to  180  deg. 
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Appendix  B 

INTEGRATING  TEE  EQUATIONS  OF  MOTION 

In  representing  the  flight  of  three  vehicles  in  motion  siaulta- 
neously,  there  ere  18  differential  equations  to  be  numerically  integrated, 
nine  involving  the  accelerations  and  the  other  nine  involving  the 

velocities  These  equations  sight  be  expressed  in  a  variety  of 

different  ferns,  each  having  certain  advantages  in  certain  sittiations. 

The  TACTICS  prograa  expresses  the  differential  equations  in  two  different 
forms:  One  Involves  the  flat-earth  representation  (less  coiiq>lex  and 
faster),  and  the  other  Involves  a  round  rotating  or  nonrotating  earth 
(essential  for  space  applications,  long  ranges,  or  high  speeds).  Unless 
specified  by  the  KINTEG  option,  GATA(llO),  the  flat-earth  equations,  are 
used.  In  order  to  decrease  execution  time,  only  12  equations  are  inte¬ 
grated  when  vehicle  2  is  in  captive  flight  (subroutine  CAPFLT). 


FLAT-EARTH  CARTESIAN  FORM 

The  accelerations  are  deterstined  by  resolving  the  applied  forces 
(e.g.,  aerodynamic,  gravitational,  or  propulsive)  into  three  components. 
Two  of  these  coaqionente,  and  normal  to  the  velocity  vector 

V,  and  their  vector  sum  is  the  net  lateral  acceleration  of  the  vehicle. 
The  direction  of  the  third  component,  V,  is  parallel  to  or  along  the 
velocity  vector  V.  When  integration  is  tc  occur,  the  three  components 
of  acceleration  well  as  the  unit  vectors  1^,  1^, 

and  ly,  have  been  determined,  so  that  the  expression  for  the  net  ac¬ 
celeration  vector  a  is 
o 

•o  ■  ‘oh  h  +  »oV  ^  «5) 

Converting  to  Cartesian  coordinates  yields  the  expression 


x*a.l,  +a„l^  +  VL 
oh  Ax  oV  Dx  va 


(56) 


and  similarly  for  y,  z.  (See  Eqs.  (9).) 
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Integratlon  the  equations  for  a^,  and  a^  corresponding  to 
each  of  the  three  vehicles  will  result  in  V^*  and  V^,  respectively. 
At  this  point,  the  use  of  relative  velocities  and  positions  is  intro¬ 
duced  to  mlninlze  errors  in  relative  position  due  to  differencing  large 
numbers.  Accordingly,  the  velocities  of  vehicles  1  and  2  relative  to 
vehicle  3  are  formed  and  Integrated  to  obtain  the  relative  positions. 
The  equations  are  as  follows: 


(57) 


Integration  of  the  above  equations  yields  the  relative  position 
vectors  r^^^  and  T^2y  velocity  of  vehicle  3,  (i.e.,  an  absolute 

velocity  relative  to  the  fixed  origin  of  the  x,  y,  z  frame),  is  inte¬ 
grated  to  yield  the  absolute  position  vector  Xy  The  absolute  posi¬ 
tions  r^^  and  are  subsequently  determined  from 


This  procedure  emphasizes  the  accuracy  of  the  relative  rather  than 
absolute  position  geometry. 

ROUND-EARTH  SPHERICAL  FORM 

The  dynamic  equations  for  a  particle  whose  motion  is  observed 
from  a  rotating  reference  frame  involve  centrifugal  and  Coriolis-force 
terms.  Moreover,  since  the  earth  is  represented  as  spherical,  a  number 
of  coordinate  transformations  are  necessary  to  relate  the  various 
frames  of  reference.  For  convenience,  the  basic  equations  of  motion 
and  certain  unit  vector  notations  have  been  extracted  from  the  ROCKET 
program^^^  in  order  to  provide  compatibility.  However,  since  T.4CTICS 
is  primarily  concerned  with  the  relative  motion  of  possibly  three 


f 
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different  vehicles,  there  are  significant  differences  between  coordi¬ 
nate  systems  and  computational  methods  employed. 

The  basic  differential  equations  of  motion  are  derived  using  three 
unit  vectors  defined  as  follows: 

1^,  directed  radially  from  the  earth's  center  to  the  vehicle,  a 
point  mass. 

1  ,  normal  to  1  and  directed  eastward  along  a  parallel  of 

la  K 

latitude. 

Ip,  directed  northward  along  the  local  meridian  so  as  to  form  a 
right-handed  orthogonal  system  with  1  and  1  . 

R  la 

The  acceleration  of  the  vehicle  with  respect  to  a  nonrotating  inertial 
X,  Y,  Z  coordinate  system  is 

a-r  =  Xi  +  Yj+Zk  (59) 


This  acceleration  a  corresponds  to  a  net  result  of  all  forces  applied, 
F/m.  The  following  derivation  develops  expressions  for  the  net  accel¬ 
eration  resulting  from  F  in  terms  of  a  rotating-earth  fixed-coordinate 
system,  l.e.,  longitude  A,  latitude  (p,  radial  distance  r,  and  time 
derivatives  thereof. 

First,  the  following  simple  device  may  be  used  to  transform  coordi- 

I  •  =  (Ij^  •  I)  X  +  ♦  J)  W  .  k)  z  (fo> 


and  similarly  for  a  *  1  and  a  •  1  . 

Li  It 


Referring  to  Fig. 

L. 

32 ,  the 

required  dot  products 

are  by  inspection 

•  i  =•  cos  <p  cos  0 

i  *  -  sin  0 

i  »  -  sin  q)  cos  0 

♦  j  »  cos  (p  sin  0 

j  *  cos  0 

j  *  -  sin  9  sin  0  (61) 

\ 

•  k  *  sin  rp 

k  -  0 

* 

k  *  cos  (p 

,  Ip  sys  tem: 


nates  to  the  1^, 


« 
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(Note  the  analogy  with  the  definition  of  unit  vectors  1^,  1^,  and  Ij^j 
see  Eq.  (6).) 

The  desired  transforoatlon  may  now  be  accomplished  using 

I  «  (I  •  y  \  +  G  •  \)  \  +  (a  •  ip)  ip  (62) 

Expanding  and  substituting  the  dot  products  results  in 

a  -  Ij^  (X  cos  <p  cos  0  +  Y  cos  <p  sin  0  +  Z  sin  (p) 

+1^  (-  X  sin  0  +  Y  cos  0+0)  (63) 

+  Ip  (-  X  sin  (p  cos  0  -  Y  sin  ip  sin  0  +  Z  cos  ^) 

The  inertial  coordinates  X,  Y,  Z  may  be  expressed  in  terms  of  spherical 
coordinates  as 

X  ■  r  cos  (p  cos  0 

Y  *  r  cos  <p  sin  0  (64) 

Z  ■  r  sin  <p 

Differentiating  these  expressions  twice  with  respect  to  time  and  sub¬ 
stituting  in  Eq.  (63)  will  yield 

_  _  ••  ^2  *2  2 

a  "  Ip  (r  -  r  9  -  r  0  cos  9) 

—  **  ♦  •  • 

+  Ij^  (0  r  cos  <p  +  2  f  0  cos  <p  -  2  r  9  0  sin  9)  (65) 

—  /•  .? 

+  Ip  (9  r  +  t  0  sin  9  cos  9  +  2  r  9) 

Referring  to  Fig.  32,  we  note  that  0  is  a  time  variable: 

0  “  X  +  0)  t 
e 


where 


X  *  an  arbitrarily  defined  angle  of  longitude  oi:  the  earth's 
surface 

0)^  »  the  earth's  angular  rate  of  rotation  (7.29211585  •  10  ^  rad/sec) 
t  =  the  time  interval  of  rotation 

Accord *ngly. 


0  =  X  +  0) 

e 


and 


0  »  X 


(66) 


(67) 


The  resultant  force  F,  causing  acceleration  a  with  respect  to  the 
inertial  Y,  Z  frame,  may  be  resolved  into  components  Fp  in 
the  Ip,  Ij^,  and  Ip  directions,  respectively.  Substituting  a  »  F/m  in 
the  preceding  equations  and  equating  components  results  in  the  desired 
set  of  equations  of  motion: 


Fp/m  +  r 


•2  ^ 

9  +  r 


(X  +  u)  ■  cos  <p 


X 

9 


jp^^/m  -  2  r  (X 

[- 


X  +  lo)  cos  9  +  2  r  9  (X  +  uj)  sin  9  /r  cos  9 


Fp/m  -  r  (X  +  lo)  sin 


.  .1 

9  cos  9  -  2  r  9  /i 


'] 


(68) 


These  equations  are  a  set  of  second-order  differential  equations 
in  terms  of  geocentric  radial  distance  r,  longitude  X,  latitude  9,  and 
components  of  resultant  force  F.  (The  methodology  of  force  definition 

is  discussed  in  Section  IV.)  Numerical  integration  of  these  second- 

$  • 

order  equations  will  yield  the  first-order  velocity  components  r,  X, 
and  9  at  a  time  t  +  At.  Theoretically,  a  second  integration  will  yield 
the  position  quantities  r,  X,  and  9,  but  we  wish  to  avoid  this  proce¬ 
dure  to  minimize  numerical  errors  arising  from  differencing  large 
numbers  associated  with  the  geocentric  distance  r.  Accordingly,  posi¬ 
tion  components  are  computed  by  the  following  method: 


4 
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O 

O 


o 


Integrating  the  second-order  equations  to  determine  r.  A,  and  qj. 
Forming  the  topocentrlc  velocities  i»e.,  the  velocities 

with  respect  to  an  arbitrarily  chosen  reference  point  (or 
coordinate  system  origin)  on  the  earth's  surface. 

Integrating  the  topocentrlc  velocities  to  determine  the 

topocentrlc  range  position  vectors  . . 


The  topocentrlc  reference  point  or  local  coordinate  system  origin 
is  defined  by  input  data  in  terms  of  latitude  <jp^  and  longitude 
so  that 


R  =  R  cos  <p  cos  (A  +  to  t) 
ox  o  o  o  e 


R  =  R  cos  G>  sin  (A  +  to  t) 
oy  o  ^o  o  e 


(69) 


oz 


R  sin  <p 
0  ^o 


where  R  is  taken  to  be  2.0925861  •-<10  ft. 

0  __ 

The  velocity  of  R^  with  respect  to  an  inertial  frame  is 


R“-R  toi  +  R  loj 
o  oy  e  ox  e 


(70) 


Note  that  in  the  second-order  differontlal  equations  the  centrifugal 
and  Coriolis  accelerations  are  a  function  of  latitude  angle  <p;  hence, 
numerical  solutions  will  be  affected  by  assumed  initial-condition  values 
for  vehicle  latitude.  On  the  other  hand,  Initial  assumptions  for  vehicle 
longitude  position  may  be  purely  arbitrary.  Figure  33  shows  the  topo- 
centric  nonrotating  coordinate  system  with  the  X^,  Y^,  axis  always 
remaining  parallel  to  the  X,  Y,  Z  inertial  axis.  With  the  R^  origin 
established,  the  topocentrlc  range  vectors  (from  the  topocentrlc  origin 
to  the  vehicles)  are 


Pi  “  ^  -  ^o 


1,  2.  3 


(71) 


and  the  velocities  are 
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acd 


im 

+  <ij)  *  r 

e 

local 


dc 


•  r  +  r  1  CCS  5»  +  r  9  Ij, 

local 


(73) 


is  *  r 

e 


r  u> 
e 


coa  9  Ij^ 


(74) 


so  that 


+  r(A  +  M^)  cos  9  1^  +  r  9  Ip  (75) 

* 

The  Integratios  cf  the  topoecntric  velocities  will  yield  the 
positions  in  terns  of  tha  Inertial  coordinates  Z„,  Although 

this  coordinate  systea  is  convenient  for  astrononical  or  space  applica¬ 
tions,  it  is  not  desirable  for  aany  other  problens  in  which  the  hori¬ 
zontal  plane  and  the  local  vertical  are  nore  significant.  Referring 
back  to  the  flat-earth  x,  y,  z  system,  ve  would  like  the  z-con^onent 
to  correspond  to  altitude  and  t'r,  z-y  plane  to  correspond  approximately 
to  the  horizontal  plane.  Accordingly,  all  coordinates  are  trans¬ 
formed  Into  an  azimuth-elevation  system  with  the  origin  at  the 
position,  the  z-axis  in  the  direction  of  and  the  y-axis  directed 
eastward,  as  shown  in  Fig.  34.  Subroutine  T0P0CN  is  used  for  trans¬ 
forming  from  either  eqxiatorlal  plane  to  azimuth-elevation  coordinates 
or  vice  versa. 

The  coordinate  transformations  are  as  follows.  The  transformation 
from  topocentrlc  X„,  Y„,  Z  coordinates  to  azimuth-elevation  x,  y,  z 
coordinates  is  given  by  the  following  matrix: 


n 
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Ibe  r-evers*  trassfsnaiftca,  £r«ra  ajIaotSs-eleyasioa  x,  y,  z  coordftates 
to  to;»5C£acric  3^,  2^  cocrdiczCes,  is  g;lvi£s  by 


1 

sin  ^  cos  5^ 

-  sin  e_ 

T 

cos  8^  cos  1 

s 

sin  sin  3^ 

cos  8^ 

sin  6  cos 

L 

-  cos  eu. 

0 

*5ln  9^ 

(77) 


where  anii  3^  are  the  latitude  and  longitude,  respecti’/ely,  of  the 
rotating  reference  origin  at  point  T,  defined  by 


-  <Pj.(V 

®T  -  +  “e  -  "o> 


(78) 


In  addition  to  the  "standard"  printout  in  terms  of  x,  y,  z  azimuth- 
elevation  (horizontal-plane)  system  coordinates,  the  quantities  of 
latitude  <p,  longitude  A,  and  altitude  h  for  each  vehicle  are  available 
as  an  optional  printout  section.  Note  that  because  of  the  earth  cur¬ 
vature  and  the  use  of  a  fixed-earth  reference  point  (R  ),  the  magnitude 

o 

of  the  z-position  component  will  not  generally  correspond  exactly  to 
altitude  (because  of  the  x-y  displacement  from  the  origin) . 


I 
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Apptsdix  C 

CgmSMtCE-AiGMamiaSL-IAK  DEFIXTTIOXS 


Ifels  appcsd  'x  coQCaitts  descrlpcioss  asd  sxachescscical  defiaitioss 
of  aosc  of  Che  slsniflcanc  guidance  asd  conCcol  laws  developed  for  use 
with  X&CZICS.  The  llsc  is  open-ended  and  izicoaplece,  since  new  ideas 
are  continuously  being  iaplcmsted  as  new  subroutines.  See  Section  V 
for  a  suanary  of  definitions  and  the  applicable  acceleration  equations. 

OPBi-LOOP  OailRDL  LAWS 

Straight  Flight 

The  coaaanded  lateral  acceleration  z„  is  zero.  The  vehicle  will 
fly  a  straight-line  path  (but  not  necessarily  ''straight  and  level"). 
However,  an  acceleration  or  deceleration  dlon^  this  path  eay  occur 
due  to  Che  thrust-drag  zelaclonship. 

Straight  and  Level  Flight 

This  routine  is  used  to  return  a  vehicle  to  level  flight  in  Che 
horizontal  plane  from  a  diving  or  climbing  condition.  An  arbitrary 
control-law  assumption  is  made  that  the  vehicle  will  comence  the 
return  to  lavel  flight  at  maximum  g  capability  subject  to  structural 
or  constraints.  The  magnitude  of  this  acceleration  tapers  off 

linearly  as  the  velocity  vector  reaches  an  angle  of  10  deg  from  the 
horizontal  plane  and  becomes  zero  within  0.1  deg  of  level  flight.  The 
direction  of  the  acceleration  is  accordingly 

T,  -  ±  \  (79) 

(positive  for  climbing  and  negative  for  diving) .  The  magnitude  of  the 
acceleration  follows  the  arbitrary  law 


(30) 


»c  *  Us^:1(t/io)»5_„I) 

a  *  0  y  <  0.1 


Captive  FUjrbt 

This  routine  Is  U2e;d  to  zero  out  coaputatlooa  and  printed  values 
for  vehicles  which  art.  In  captive  fllj^it.  There  are  three  nodes : 

o  Vehicle  2  locked  to  vehicle  1. 

o  Vehicle  1,  2,  or  3  locked  to  the  zero  origin  (l.e. « 
all  values  for  position,  ^^ilocity,  and  acceleration 
are  zero). 

o  Vehicle  2  locked  to  vehicle  3. 

Captive  flight  is  not  a  guidance  lav  in  the  sense  of  the  preceding 
discussion  but  rather  a  device  to  elisinate  unnecessary  ccaputations 
and  Isprove  the  appearance  of  the  printed  values. 

Launch 

This  control  .aw,  which  simulates  the  launch-~boost  phase  of  a 
missile  flight,  may  be  applied  to  any  of  the  three  vehicles.  The  call 
for  "launch"  is  usually  based  on  some  criteria  stated  in  P0LICY  (e.g., 
range,  range  rate,  geometry,  accelerations,  time,  and — most  Importantly — 
combinations  thereof).  When  this  routine  is  called,  the  boost  velocity 
AV  must  be  specified  as  a  constant.  The  commanded  lateral  acceleration 
Is  gravitational  only: 


-  -g  cos  Y  Ip 


<81) 


Aim  (1) 

This  routine  was  provided  for  applications  in  which  it  is  necessary 
to  simulate  the  aiming  process,  e.g.,  SAM  launch.  It  involves  the 
following  problem:  Given  a  target  velocity  vector  and  a  vehicle 
speed  (both  assumed  to  remain  constant)  separated  by  a  range  vector 
r^j,  find  the  orientation  of  i.e.,  aiming  angles,  so  that  an 


intercept  will  occur.  Idke  ceptive  It  is  a  device  for  defining 

angles  rather  than  a  lav  for  defining  a^.  ?ke  solution  is  outlined 
below: 


r  +  V  T  *  V  T 
ij  T  M 


(82) 


where  T  is  the  "time-co-go"  until  impact.  The  magnitude  of  is 

fl 

known,  but  the  components  and  are  to  be  determined. 

From  Che  above  equation. 


?  7\^ 


,I|  -  .  2ry  ■  V,I  +  vyj 

(  -  b  ±Vb^  -  4ac) 


(83) 


T  >  0 


2a 


where' a,  b,  and  c  are  the  corresponding  coefficients  in  the  quadratic 
equation.  Hence, 


V  «  — i.1  4.  V 
M  T  T 


\x  “  T  '*■  ^Tx 


(84) 


and  similarly  for  and 
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The  alaing  angles  tp.  and  6  nay  be  detcrnined  from 

A  n 


This  routine  is  similar  in  purpose  to  the  previous  item.  In  cer¬ 
tain  problems  it  is  necessary  to  simulate  the  launching  of  a  missile 
(or  firing  of  a  projectile)  along  the  LOS  to  the  target.  The  aiming 
process  is  more  compiler,  ted  if  the  launching  platform  motion  is  taken 
into  account.  In  order  for  the  projectile  to  travel  along  the  LOS 
(ignoring  ballistic  drop) ,  the  orientation  angles  of  the  vesuttant 
projectile  velocity  and  platfcrm  .-elocity  must  be  the  same  as  those  of 
the  LOS.  Hence,  if  is  tite  projectile  velocity  and  is  the  launch¬ 
ing  platform  velocity,  where  is  the  resultant  veloc¬ 

ity  from  which 


(86) 


The  magnitude  of  determined  from 

"L  -  =°’(\0S  -  *Vl)  “'(’los  '  m)  -  *1  -  i  ■  0  (87) 

Knov;lng  the  magnitude  and  orientation  of  'v^2' 


V 

V 

V 


2x 

■  V 

12x 

-"lx 

■  V  - 

-  V 

2y 

I2y 

ly 

"  V  ^ 

-  V 

2z 

12z 

Iz 

(88) 
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The  orientation  angles  of  are 

{\2  ’los  -  '’lx) 

Y  m  ..  1,1-  1,1,  fm: 

^V2  V2 

(89) 

g  ('■'12  Ijos  %3  -  ''2v) 

(''12  “•  'fias  ®L0S  -  ''2.) 

and  the  sagnitude  is  the  boost  velocity 


IV^I  -  AV 

The  incremental  angles  which  should  be  added  to  the  la  inching  platform 
angles  for  aiming  are  then 


Ay, 


V2 


'V2 


-  Y, 


VI 


AG. 


V2 


V2 


-  e 


VI 
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Left  or  Right  Turn  (1) 

The  commanded  lateral  acceleration  vector  a  is  in  the  horizontal 
plane  and  has  a  constant  value  ac  specified  when  calling  the  routlne(s). 
The  two  routines  (left  and  right)  arc.  identical  except  for  an  algebraic 
sign  corresponding  to  the  direction  of  the  turn  (±  1^).  It  was  decided 
to  specify  the  magnitude  of  the  turning  acceleration  in  terms  of  the 
resultant  normal  acceleration,  expressed  in  number  of  g’s  or  (see 

Eq.  (21)).  Accordingly, 

“c  ■  (w)y 

Left  or  Right  Turn  (2) 

The  commanded  latssral  acceleration  vector  a^  is  in  the  horizontal 
plane.  The  magnitude  a^  depends  on  the  excess  thrust  available  to 
perform  a  turn  with  constant  Mach  number.  In  other  words,  these 
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< 


maneuvers  arc  conscant-Mach-number  horizontal  turns;  a  necessary  aero¬ 
dynamic  condition  for  initiating  such  a  turn  is  that  thrust  must  exceed 
drag,  because  additional  drag  will  be  induced  by  the  maneuver  itself. 
These  routines  require  Iterative  processes  if  aerodynamic  tables  are 
involved,  since  the  tables  must  be  read  in  reverse  (see  item  2  belov). 

A  brief  outline  of  the  problem  is  as  follows: 

a  a 

1.  Since  Hach  number  is  constant  with  altitude,  H  V  •*  0,  and 
the  necessary  aerodynamic  relationship  if 


D  +  W  sin  Y  +  T  cos  a  »  0 


(92) 


from  which  D  may  be  calculated.  Knowing  D, 


"D  Aq 


(93) 


2.  Find  and  a  correaponding  to  and  Mach  number.  For  example, 
using  the  familiar  analytic  expressions 


2 

S "  S  — r-  \ 

°  %  diet)  ^ 


da 


s. - 


'(S  -  S 

o. 


d(cj) 


(dCj^/da)  ^  % 


(94) 


T,  *  aAq 

ia 


F  •  L  +  T  sin  a 
n 


The  equivalent  of  the  preceding  operations  must  be  performed  if  tabular 

2 

values  are  used.  In  effect,  values  for  dC^^/da  and  dCp/d(Cj^)  must  be 
numerically  determined. 
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3.  Fled  toe  oagnitude  sad  dlreccion  of  corresponding  to  the 

airailable  foxce  F  . 

a 

‘C  *  (v)/^n  ■  V 

idiere  (4-)  is  i;sed  for  a  left  turn  and  (-)  for  a  right  turn.  Note  that 

the  Bagnitude  of  a.  is  not  explicitly  specified  but  is  dependent  on  F 
c#  n 

and  Y»  i«e«>  it  is  not  known  a  priori. 

Left  or  Right  Turn  (3) 

These  maneuvers  are  similar  to  the  Left  or  Right  Turn  (2)  maneuvers 
described  above  except  for  the  constant-altitude  requirement.  By  means 
of  a  dive,  gravity  may  be  utilized  to  hold  Mach  nianber  constant  during 
the  Cum.  If  there  is  insufficient  thrust  to  meet  the  conditions  for 
constant  Hach  nianber  and  a  specified  number  of  g's  In  Che  turn,  the 
vehicle  will  dive  in  order  to  satisfy  these  conditions.  On  Che  other 
hand,  if  excess  thrust  is  available  (or  after  a  sustained  dive),  the 
vehicle  will  climb.  If  a  maximum  lift  coefficient  (C_  )  limit  is 

LlDHX 

reached,  the  vehicle  will  turn  at  the  corresponding  allowable  accel¬ 
eration  for  C-  ,  and  Mach  number  will  equal  a  constant.  Since  i-Iach 
Lmax 

number  varies  with  altitude,  the  following  relationship  must  be 
satisfied: 


•  e 


Vs 

M  ■  V/s - 2  “  ^ 

s 


(96) 


Referring  to  the  relationships  given  in  Appendix  I  and  performing  the 
necessary  arithmetical  operations, 

p  -  -  6.9697419  x  10"®(1  -  6.8865741)^*^^®  z 


p  -  -  32.207899  pz 


a  «  0 


(97) 
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From  the  preceding  expression  for  M  =  0, 

V  -  Vs/s  0 


(98) 


The  necessary  aerodynamic  relationship  is 


D  +  W  sin  Y  +  T  cos 


(99) 


from  which  the  drag  D  may  be  determined  by  means  of  the  applicable 
expressions  as  given  by  Eq.  (26).  The  components  of  commanded  accel¬ 
eration  are 


acv  -./(g/W)  ^ 


(100) 


where  F’  (or  rather  F’/W  in  terms  of  number  of  g’s)  is  to  be  specified 
n  n 

when  the  maneuver  is  call'i.J  and  F  corresponds  to  the  L,  C_,  C,  ,  and 

n  u  L 

a  relationships  of  Eq.  (27).  The  component  is  positive  for  left 
turns  and  negative  for  right  turns.  Special  cases  arise  if  the  speci¬ 
fied  acceleration  F^/W  exceeds  structural  or  aerodynamic  constraints. 

If  F'  should  exceed  the  structural  constraints,  the  value  is  altered 
n 

to  correspond  to  the  maximum  value  a„  .  If  F*  exceeds  the  aerody- 

Smax  n 

namic  value 


F 

nmax 


Lmax 


Aq  +  T  sin  a 

max 


(101) 


then  F'  is  altered  to  correspond  to  F  .  Since  C,  varies  with 
n  nmax  Lmax 

Mach  number  and  by  definition  the  Mach  number  should  not  change,  F' 

n 

will  remain  equal  to  F  under  this  constraint. 

^  nmax 


Left  or  Right  Turn  (4) 

These  maneuvers  are  identical  to  the  Left  or  Right  Turn  (3)  maneu¬ 
vers  except  under  the  F  constraint  discussed  above.  Since  C, 

nmax  Lmax 
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varles  with  Mach  number,  a  correspondence  must  be  established  between 
the  specified  Mach  number  and  the  specified  number  of  g's  in  the  turn. 
The  Left  or  Right  Turn  (3)  maneuvers  l^jnit  the  g's  to  the  corresponding 
Mach  number  and  condition.  The  Left  or  Right  Turn  (4)  maneuvers 

allow  the  vehicle  to  dive  and  increase  Mach  number  (and  hence  C,  ) 
in  order  to  obtain  the  specified  ntniber  of  g's  in  the  turn.  That  is, 
if  the  specified  F'  exceeds  F  ,  the  constant-Mach-number  maneuver 
requirement  is  initially  abandoned  to  increase  speed  and  F^^^^  so  that 

v  ^  v  * 
nmax  “  n 

Under  the  constraint  conditions,  the  initial  dive  is  defined  arbitrarily 
by 


Left  or  Right  Turn 


These  are  climbing  or  diving  turns  as  specified  by  arguments  of 

the  number  of  g's  required,  i.e.,  F  /W  and  the  bank  angle  As  shown 

n  D 

in  Figs.  7  and  8,  the  commanded  acceleration  components  are 


A  right  turn  results  for  positive  values  from  0  to  180  deg  and  a  left 

turn  results  for  negative  i()„  values  from  0  to  -180  deg.  Diving  or 

o 

climbing  occurs  for  absolute-value  magnitudes  greater  or  less  than  90 
deg,  respectively. 


f 
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Cllmb  or  Dive  (1) 

These  tnaneuvers  are  identical  to  the  Left  or  Right  Turn  (1)  de¬ 
scribed  above  except  that  the  acceleration  vector  a  Is  in  the  vertical 

0 

plane .  That  is , 


The  plus  sign  corresponds  to  the  climb  maneuver  and  the  minus  sign  to 
the  dive  maneuver. 

Climb  (2) 

This  maneuver  Is  a  constant-Mach-nuciber  climb  depending  on  excess 
thrust  available  at  the  time  it  is  initiated.  The  procedure  is  the 
same  as  that  for  Left  or  Right  Turn  (2)  described  above  except  that 
the  vector  a^  is  in  the  vertical  plane.  That  is, 

F  \  _ 

~  -  cos  Y  J  Ijj  (105) 

where  is  not  specified  but  rather  calculated  from  the  relationships 
described  in  Left  or  Right  Turn  (2)  and  (3)  (see  Eq.  (94)). 


Barrel  Roll  (1) 

In  attempting  to  describe  the  trajectory  of  a  barrel  roll  within 
the  framework  of  coordinated  turn  definition  (as  described  in  Section 
IV) ,  there  does  not  seem  to  be  universal  agreement  among  pilots  and 
engineers  as  to  precise  mathematical  definition.  The  following  equa¬ 
tions  were  extracted  from  a  Target  Generator  Program  received  from 

,  * 

Eglln  Air  Force  Base: 


A 

Private  communication. 


where  the  bank  angle  i}.-  Is  obtained  from  a  specified  banking  rate 

D  15 

so  that 

’I'b  *  J  5’b  (107) 

Referring  to  Fig.  8,  we  note  that  the  force  rotates  about  the  veloc¬ 
ity  veci.or  V.  The  accele .itlon  component  a^j^  varies  sinusoidally,  where¬ 
as  the  component  a^^  has  a  varying  gravitational  term  g  cos  y-  In 
calling  for  Barrel  Roll  (1) ,  it  is  necessary  to  specify  (1)  the  number 
of  360-deg  rolls  required,  (2)  the  number  of  g's  for  and  (3)  the 

banking  rate  ifi  in  degrees  per  second. 

O 

Barrel  Roll  (2) 

In  order  to  decrease  the  altitude  loss  due  to  gravitational  effects, 
this  routine  assumes  that  the  net  acceleration  vector  a^  rotates  about 
the  velocity  vector  V  at  a  constant  specified  rate  ij)'.  The  net  accel¬ 
erations  in  the  horizontal  and  vertical  planes  are 


*Ch  '  " 


®cv  “ 


(108) 


where  is  given  by 


S  dt 


(109) 


Referring  to  Fig.  8,  we  see  that  the  acceleration  vector  a^  will  rotate 
uniformly  about  the  velocity  vector  V,  but  tne  force  F^  will  vary  with 
the  rotation.  F^  is  given  by 


F 

n 


© 


4h-^ 


(a^^  +  g  cos  y)‘ 


(110) 


In  calling  for  Barrel  Roll  (2),  it  is  necessary  to  specify  (1)  the 

number  of  360-deg  rolls  required,  (2)  the  number  of  g's  for  a„  (not 

« 

F  /W) ,  and  (3)  the  angular  rate  il>„  in  degrees  per  second, 
n  D 


f 
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Constant-Roll-Angle  Climbing  or  Diving  Turn  (GRP(?LL) 

ITie  term  roll  an^le,  as  used  here,  is  an  angle  referenced  to  the 

vehicle's  pitch,  roll,  and  yaw  coordinate  system,  as  distinguished  from 

* 

bank  angle  i|)„,  which  is  referenced  to  a  wind  or  velocity  system.  In 

D 

TACTICS,  the  roll  angle  ili  is  taken  to  be  the  angle  between  the  vehicle 
pitch  axis,  i.e.,  the  wings  of  the  vehicle,  and  a  line  through  the 
c.g.  of  the  vehicle  both  normal  to  the  longitudinal  (i.e.,  roll)  axis 
and  parallel  to  the  horizontal  plane.  The  constant-roll-angle  maneuver 
requires  this  angle  to  be  held  constant.  The  normal  force  in  g's,  Fj^/W, 
is  also  specified.  For  values  of  ip  between  -90  deg  to  -180  deg  and  be¬ 
tween  +90  deg  Co  +180  deg,  a  diving  turn  to  the  left  or  right,  respec¬ 
tively,  will  result. 

Since  and  t  (see  Appendix  A  for  definition)  are  specified,  the 
subroutine  solves  the  problem  of  finding  the  components  a^j^  and  a^^  to 
force  the  following  two  conditions: 


and 


=  constant 


(111) 


ip  ■  *oV^  "*  (112) 

The  functional  relationship  1*  not  amenable  to  solving  these 

two  equations  for  a^j^  and  a^y  directly,  so  numerical  methods  are  used. 

The  technique  employed  is  to  rotate  the  lift  vector  L  about  V  by  an  amount 
(the  absolute  value  |l(  Is  known  but  the  direction  is  not).  With 
each  Incremental  rotation  the  angle  ip  is  computed  until  the  required 
conditions  are  fulfilled  within  an  error  AiJ^  less  than  1  mrad.  Note  that 
by  holding  a  roll  angle  of  absolute  magnitude  greater  than  zero,  the 
vehicle  will  continuously  be  climbing  or  diving  until  the  longitudinal 
roll  axis  becomes  vertical.  However,  the  velocity  vector  will  be 


Apparently,  there  is  no  universal  agreement  in  texts  or  among 
engineers  on  terminology  or  exact  definition. 
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lagging  by  approximately  the  angle  of  attack  a.  When  this  condition 
is  reached,  roll  angle  becomes  singular  and  a  warning  message  is  printed. 
Thereafter,  the  vehicle  will  continue  its  climb  (or  dive)  until  the 
velocity  vector  becomes  vertical;  then  straight  flight  up  (or  down)  is 
assumed . 

Ballistic  20-mm  Projectile 

This  routine  simulates  the  trajectory  of  a  20-mm  type  M56A1  pro¬ 
jectile.,  The  drag  coefficient  versus  Mach  number  characteristics 
are  represented  by  several  linear  functions  over  various  Mach-number 
regimes.  The  commanded  accelerations  correspond  to  a  0-g  ballistic 
trajectory! 


a.  =  0.0 
h 

=•  -  g  cos  Y  (113) 

F  »  0.0 
n 

The  projectile  weight  is  0.22  lb,  and  the  reference  area  is 
0.0033842  ft^.  The  muzzle  velocity  ranges  from  3350  to  3450  ft/sec 
and  should  be  set  as  4V  in  the  launch  subroutine  (see  Launch  above) . 

CLOSED-LOOP  CONTROL  LAWS 

Proportional  Navigation 

The  commanded  lateral  acceleration  is  proportional  to  the  space 
rate  of  rotation  of  the  LOS  between  missile  and  target.  Expressed  in 
vector  notation, 


a^  -  XVo)  X  1  (114) 

C  r  V 

where 

>\  >  the  "navigation  constant"  (it  may  be  treated  as  either  a 
constant  or  a  variable) 
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V  »  missile  speed 

oj^  »  relative  angular-rate  vector  as  defined  by  Eq.  (5)  in 
Section  III. 

1^  ■  unit  vector  along  the  missile  velocity  V,  i.e.,  ly  *  V/V 
The  direction  of  the  acceleration  may  be  defined  by 


u  ^  1,,/w 
r  V  r 


(115) 


The  commanded  acceleration  a  may  be  resolved  into  horizontal  and 
vertical  components  by 


®Ch  “ 


*CV  “  *C 


(‘l  •  \)  ■ 


(116) 


Biased  Proportional  Navigation 

The  commanded  lateral  acceleration  a^  is  proportional  to  a  biased 
(or  "modified")  space  rate  of  rotation  of  the  LOS  between  interceptor 
and  target.  The  bias  term  may  be  either  a  constant  or  a  time  variable; 
and  its  general  purpose  is  (1)  to  provide  for  accelerations  involved  in 
the  problem  and/or  (2)  to  shape  (l.e.,  straighten  or  curve)  the  resul¬ 
tant  trajectory.  The  calculations  to  determine  this  quantity  may  range 
in  complexity  from  the  estimation  of  a  constant  to  the  solution  of  com¬ 
plicated  second-order  prediction  equations.  The  present  version  of 
this  routine  uses  a  bias  term  to  account  for  an  average  boost  velocity 
increment,  such  as  would  occur  for  a  high-specJ.flc-lmpulse  missile 
launch.  It  may  be  used,  for  example,  in  simulating  a  lead  collision 
course  where  the  interceptor  aircraft  is  to  steer  so  as  to  account  for 
the  AV  occurring  at  missile  launch.  The  equations  used  are 


'("t  -  “b)  " 


(117) 
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where 


/aV  r  X  1  \ 

X  IJ  V/ 
2 


(118) 


r 


AV  »  terminal  boost  velocity 

=  relative  Interceptor-target  range  vector 

ly  =  unit  vector  along  the  Interceptor  velocity 


Note  that  the  bias-term  evaluation  described  above  requires  an  a  priori 
knowledge  of  only  AV  and  LOS  orientation  and  hence  requires  minimal 
hardware  on  board  the  vehicle. 


Lead  Collision 

This  guidance  law  may  be  considered  as  a  more  sophisticated  ver¬ 
sion  of  biased  proportional  navigation.  Since  range  and  range-rate 
Information  are  assumed  available,  it  is  feasible  to  make  more  accurate 
predictions  in  providing  for  the  missile  launch.  Moreover,  the  effec¬ 
tive  gain,  1,  may  be  varied  as  a  function  of  a  computed  time-to-go,  T. 
The  miss  vector  (lead  collision)  form  of  guidance  is  briefly  described 
in  the  following  paragraphs.  Consider  the  following  miss  vector  dia¬ 
gram  (the  vectors  are  not  necessarily  co-planar) ,  in  which 

T  *  time-to-go,  e.g.,  from  "now"  to  the  missile  launching  point 

M  »  a  miss  vector — to  be  driven  to  zero. 

F  =  vector  representing  the  distance  and  direction  a  missile  will 
travel  after  launch  in  a  flight  time  t^ 

V  “  launching  aircraft  velocity  vector 

r 

V^  “  target  velocity  vector 

r  »  r^j ,  the  relative  range  vector  (LOS) 


F  =  1 


a 
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The  vector  M  nay  now  be  resolved  into  components  parallel  a^td 
perpendicular  to  the  LOS  vector  r.  First,  consider  the  parallel  com¬ 
ponent,  designated  : 


M||«K*1  -  r-  F‘l  +V_*1  +  rT 

II  r  r  F  r 


(123) 


since 


'  -(vt  -  Vp)  •  1, 


(124) 


Next  arbitrarily  set  Mii  -0  and  solve  for  T,  the  time-to-go. 


-  "f)  •  ^  - 


r  •  1  -  r 

f  r 


(125) 


since 


F  -  (Vp  +  Avjtp  (126) 

For  notational  convenience  we  will  define 

F*  -  AVtp  Ty  (127) 

Because  of  the  arbitrary  definition  of  T  in  Eq.  (125),  the  parallel  com¬ 
ponent  of  miss  M|j  is  zero  at  T.  The  problem  is  reduced  co  finding  an 
expression  for  Interceptor  lateral  acceleration  so  that  at  time  T 
is  driven  to  zero  as  shown  below: 
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Note  that  the  lateral  acceleration  la  by  definition  nonaal  to  The 

lateral  acceleration  coaaand  nay  (arbitrarily)  be  made  proportional 
to 


~  5 


(128) 


where 

“  V  T  +  F* 

F 

a^  -  AV6 

where  X  represents  the  guidance-loop  gain.  The  magnitude  and  orienta¬ 
tion  of  is  determined  as  follows: 

-  1m  X  T^l  -  l-F*  X  +  t(v^  -  Vp)  X  1^1  »  j-F*  X  Tj.  +  wrlj  (129) 

The  orientation  of  the  acceleration  a^  is  given  by 

V/ 

(m  X  I  )  X  I 

^  (130) 

The  complete  expression  for  may  now  be  written  as 


“C  ■  'yT  F*  [-(f*  >-  1,)  +  »  rl]  X  1,  (131) 

This  relationship  may  be  expressed  in  a  different  form  to  show  the 
analogy  with  biased  proportional  navigation: 


(132) 


where 


I 
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A’ 


Ar" 


VpT  +  F* 


“b  - 


pa  X  1 


ri 


(133) 
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Koce  chat  A’  and  sa  are  tiE«.-varying  functions- 
Pure  Pursuit  course  (1) 

The  coaoanded  lateral  acceleration  is  proportional  ro  the  angu¬ 
lar  differeuce  iy  betveen  the  interceptor’s  velocity  vector  V  and  the 
interceptor-target  range  vector  This  saay  oe  expressed  as  follows; 


Ay  =  sin  ^jl^.  X  1^1 

^1  “  [h  ^r) 

=  AVAyl^  (135) 

®Ch  “  ^c(^i  '  ^a) 

®cv  *  “c(^i  '  -d) 

It  is  significant  to  mention  that  a  pursuit-course  navigation  law 
is  not  usually  applied  for  missile  terminal  homing  guidance,  since  ja^j 
may  become  infinitely  large  as  ;he  relative  "ange  approaches  zero. 
However,  this  control  law  is  very  useful  for  describing  fighter  air¬ 
craft  "gunsight  aiming"  flight  paths.  The  gain  constant  A  is  arbi¬ 
trarily  selected,  since  it  represents  the  pilot's  reaction  and  skill 
in  keeping  the  target  in  the  crosshairs  (practical  values  range  from 
about  4  to  40) . 


Pursuit  (2) 


Pure  pursuit-course  navigation  is  defined  in  terms  of  maintaining 
the  velocity  vector  of  the  pursuing  vehicle  along  the  LOS  to  the  target. 
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In  practice,  a  pilot  would  only  be  able  to  approximate  a  pure  pursuit 
course  without  angle  of  attack  or  velocity  information.  Rather  than 
the  velocity  vector,  the  most  convenient  reference  is  the  airframe  it¬ 
self,  l.e.,  the  longitudinal  axis.  With  modern  jet  aircraft  in  high-g 
maneuvers,  angle-of-attack  differences  may  range  from  15  to  20  deg. 

The  Pursuit  (2)  routine  may  be  termed  a  deviated  pursuit  course  where 
the  bearing  angle  B  between  the  longitudinal  axis  of  the  pursuing  ve¬ 
hicle  and  the  LOS  Is  maintained  at  a  near-zero  value.  A  necessary  con¬ 
dition  for  beginning  the  deviated  pursuit  guidance  is  that  the  magni¬ 
tude  of  the  bearing  angle  be  less  than  which  corresponds  to  the 

*^Lmax  routine  calls  for  pure  pursuit  guidance  until 


The  trajectory  we  wish  to  simulate  is  described  mathematically  as 
follows:  Find  an  a^^  and  an  a^  resulting  in  a  net  acceleration  vector 
a  and  a  normal  force  vector  F^.  This  normal  force  will  require  a  lift 
vector  L  with  a  corresponding  angle  of  attack  a.  The  problem  requires 
solving  for  the  three  unknown  components  of  F^,  all  functions  of 
and  a.  There  are  two  functions  which  should  be  minimized: 

^1  "  Y  (136) 

where  y  Is  the  angle  between  V  and  the  LOS  and  a  is  the  angle  between 
V  and  the  longitudinal  axis  1^,  and 

•2  (*h’“v)" 

In  other  words,  we  would  like  a  =  y  and  |b(  “0  for  an  Idealized 
trajectory.  Since  nonlinear  relationships  (i.e.,  table  values)  and 
transcendental  relationships  are  Involved,  numerical  iterative  tech¬ 
niques  are  used  to  determine  the  acceleration  a  and  its  components  a^^ 
and  a^.  As  a  first  approximation,  the  proportional  navigation  law  is 
used  to  obtain  a  first  guess  for  1^.  Next,  the  ~  “  correspondence 
is  established  for  a  ■  y,  thereby  determining  L  and  F^.  Computations 
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are  then  performed  to  obtain  a^,  a^,  and  ultimately  B.  The  vector 
is  then  rotated  by  a  small  angle  ^  to  minimise  |Bj. 


On-Off  (1) 

This  guidance  law  (sometimes  called  hang-bang)  is  analogous  to 
proportional  navigation.  The  direction  of  the  lateral  acceleration 
unit  vector  1^  is  determined  identically.  However,  the  magnitude  of 
the  acceleration  is  i.ot  proportional  to  the  angular  rotation  rate  |u^| 
of  the  LOS,  but  dual-valued.  That  is. 


a  “  0.0  when  0  <  ^ 


(138) 


a  »  a 


Smax 


when  ju^l  >S 


where  ^  is  some  small  specified  threshold  value  of  angular  rate  (e.g., 
1  mrad/sec)  required  for  stability. 


On-Off  (2) 

This  guidance  law  is  analogous  to  pure  pursuit-course  navigation, 
and  the  direction  of  is  determined  identically.  As  in  Eq.  (138), 

a  *  0.0  when  0  <  |y1 

®  "  *Smax  "hen  IyI  (139) 

where  ^  is  some  small  specified  threshold  value  of  angle  (e.g.,  1  mrad) 
required  for  stability. 


Missile  (X) 

The  commanded  lateral  acceleration  is  proportional  to  the  space 
rate  of  rotation  of  the  LOS  between  Interceptor  and  target,  l.e.,  pro¬ 
portional  navigation.  Many  of  the  significant  guidance  and  aerodynamic 
parameters  vary  with  M.ich  number  and  are  unique  to  this  hypothetical 
missile  design,  as  ar-i  the  boost/burn/guide-time  intervals.  Accordingly, 


these  detelled  characteristics  are  packaged  into  s  specialised  sub¬ 
routine.  The  pertinent  characteristics  i  igcciated  vith  this  ad.ssile 
performance  are  listed  below. 


Constants 

A  -  0.13635  ft^ 

V  (initial  weight)  «•  187  lb 
o 

W  (burnout  weight)  *  125  lb 

D 

t_  (burn  time)  *  4.75  sec 

W  -  (125  -  187) /4. 75  -  -13.06  ib/sec 
Thritst  -  3053  lb 

Time  Constants 

» 

*  0.1 

T2  *  0.15 

•  0.15  (Pq/p) 

Parameters 

3/2 

A  -  0.35  +  12.22IH 

15. OM  (p  /p) 

^Smax  *  - 0~95 -  *  (16.52  -  1.52  M)  (M  »  Hach  nmsber) 

o 

2 

The  values  of  dC. /do  and  dC  /d(C^)  with  respect  to  Mach  nuniber  are 

id  1/  l« 

given  in  the  tables  on  the  following  page: 
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dCj^/dc  (per  rad) 

M^ch  Ko. 

23.52  +  J.09.9/M^ 

i  1.85 

55.5 

i  1.5,  <  1,85 

12  +  45M 

i  1.1,  <  1.5 

37.5 

<  1.1 

0.51  +  2.52/M^ 

i  1.5 

2.19 

i  1.0,  <  1.5 

0.515  +  1.675M 

i  0.6.  <  1.0 

1.5 

<  0,6 

dCjj/d(C^)  (per  rad) 

Mach  Ko. 

10  +  93/M 

i  1.5 

72 

<  1.5 

24.33  +  43.33M 

0.5,  <  1.1 

46 

<  0.5 

f 
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Appendix  D 

INSTRUCTIONS  FOR  CALLING  OPTIONAL  SUBROUTINES 


Many  optional  subroutines  are  available  with  TACTICS,  Including 
maneuver  and  guidance  routines.  (Launch  and  captive-flight  routines 
are  considered  special  cases  in  this  category,  as  shown  in  Table  2.) 

When  a  maneuver  or  guidance  subroutine  is  called,  the  following 
three  arguments  must  always  be  specified  (unless  the  routines  are 
formulated  to  represent  specific  vehicles  such  as  MISILX) : 


I; 

LA.ER0: 


ITHR: 


In  addition, 
DELV: 
M0DE: 


EPSL0N: 
LEVEL (I): 


GF0RC : 


Vehicle  to  be  used  (1,  2,  or  3). 

Type  of  aerodynamics  to  be  used. 

IAER0  =  1,  analytic  functions. 

IAER0  =  2,  aerodynamic  tables. 

IAER0  =  3,  VD0T  =  0.0. 

Thrust  to  be  used  (lb). 

ITHR  =  i,  afterburner  thrust  obtained  from  tables. 

ITHR  =  2,  military  thrust  obtained  from  tables. 

ITHR  =  3,  afterburner  thrust  =  constant  (DATA  96,  101,  106). 
ITHR  =  4,  military  thrust  =  constant  (DATA  95,  100,  105). 
ITHR  =  5,  thrust  =  0.0. 
the  following  arguments  may  be  required; 

Boost  velocity  of  missile  (ft/sec). 

Flag  indicating  captive-flight  option. 

M0DE  =  1,  holds  missile  in  captive  flight  on  fighter. 

M0DE  =  2,  sets  all  quantities  related  to  designated 
vehicle  equal  to  zero. 

M0DE  =  3,  holds  missile  in  captive  flight  on  target. 
Threshold  value  used  in  on-off  control  laws  for  stability. 
Flag  set  in  subroutine  STRLVL  used  to  communicate  to  P0LICY 
that  the  vehicle  velocity  vector  is  in  a  horizontal  plane 
within  a  tolerance  of  0.002  rad. 

Total  lateral  acceleration,  including  gravitational  effects, 
for  a  maneuver  (g's). 


-140- 


Tab  ie  2 

MANEUVER  AND  GUIDANCE  ROUTINES 


Routine  and  Argument  Listing 

Function 

BIASPN(I,  DELV,  IAER0,  ITHR) 

Biased  proportional  navigation 
guidance  law 

BRLRLld,  R0LLS,  IR«)LL,  GF0RC, 

R0LLRT,  IAER<#,  IIHR) 

Constant-g  barrel  roll  where  the 
force  Fjj  rotates  about  veloc¬ 
ity  vector 

BRLRL2(I,  R(?LLS,  IR0LL,  GP0RC, 

R0LLRT.  IAER'<,  ITHR) 

Constant-g  barrel  roll  where  the 
lateral  acceleration  rotates 

about  the  velocity  vector 

B20MM(I) 

Simulating  ballistics  for  20-inm 
cannon  projectiles 

CAPFLTd,  MlJDE) 

Captive  flight 

CLIMBld,  GPtfRC,  IAER0,  ITHR) 

Simple  cons tea c~o  climb 

1 

CLIMB2®d,  IAER0,  ITHR) 

Cons tan t-Mach  climb  depending  on 
available  excess  thrust 

DIVEld,  GPSiRC,  IAER0,  ITHR) 

Simple  constant-g  dive 

GFR^LLd,  GFCJRC,  R0LLL,  lAERjJ,  ITHR) 

Desired  GF0RCE  AND  RgiLL  (GF0RC  and 

R0LLL  in  argument  listing)  are 
specified  and  corresponding  ac¬ 
celerations  are  computed 

LAUNQId,  DELV,  IAER(?,  ITHR) 

Launches  missile 

LEADCLd,  DELV,  IAER(i,  ITHR) 

Lead  collision  guidance  law 

LTRNKI,  GP0RC,  IAER0,  ITHR) 

Simple  constant-g  left  turn 

LTRN2®(I,  IAER0,  ITHR) 

Constant-altitude,  cons tan t-Mach 
left  turn  depending  on  avail¬ 
able  excess  thrust 

LTRN3®d,  GF0RC,  IAER0,  ITHR) 

Cons tan t-Mach  left  turn  where 

GF0RC  is  specified 

LTRNA^d,  GFCJRC,  IAER0,  ITHR) 

Cons  tan t-Mach  left  turn  unless 
aerodynamic  flight  conditions 
(CLMAX)  prohibit  obtaining  spec¬ 
ified  value,  in  which  case  air¬ 
craft  will  maneuver  at  maximum 
acceleration  in  conformity  with 

CIMAX  limitations,  keeping  Mach 
number  constant 

These  subroutines  require  the  routine  MACHR8. 
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Table  2  (continued) 


Routine  and  Argument  Listing 

LTRN5®(I,  GP0RC,  R0LLL,  IAER0, 
ITHR) 

MISILX(I) 

MISIL2(I) 

0N0FF(I,  EPSL0N,  IAER0,  ITHR) 
0N0FF2(I,  EPSL0N,  IAER0,  ITHR) 
PR0NAV(I,  IAER0,  ITHR) 


Function 


Cons tan t-g  climb ing  or  diving 
left  turn  as  defined  by  bank 
angle  R0LLL  (see  Appendix  C) 

Proportional  navigation  missile 

Reaivard-launched  MISILX(I) 

Bang-bang  (on-off)  control  law, 
analogous  to  PR0NAV 

Bang-bang  (on-off)  control  law, 
analogous  to  PRSUIT 

Proportional-navigation  guidance 
law 


PRSUIT (I,  IAER0,  ITHR) 
PRSUIT2(I,  IAER0,  ITHR) 


Pursuit-course  navigation,  veloc¬ 
ity  vector  pointed  down  LOS 

Pursuit-course  navigation,  thrust 
vector  pointed  down  LOS;  must 
be  used  with  subroutine  FUNCTN(I) 
(See  Appendix  C) 


RTRNld,  GF0RC,  IAER0,  ITHR) 
RTRN2^(I,  IAER0,  ITHR) 

RTRN3®(I,  GF5IRC,  IAER0,  ITHR) 
RTRN4®(I,  GF0RC,  lAERSf,  ITHR) 


RTRNS^d,  GF0RC,  R0LLL,  IAER0, 
ITHR) 

STRFLTd,  IAER0,  ITHR) 


Simple  constant-g  right  turn 

Constant-altitude,  constant-Mach 
right  turn  depending  on  avail¬ 
able  excess  thrust 

Constant-Mach  right  turn  where 
GF0RC  is  specified 

Constant-Mach  right  turn  unless 
aerodynamic  flight  conditions 
(CLMAX)  prohibit  obtaining 
specified  value,  in  which  case 
aircraft  will  maneuver  at  maxi¬ 
mum  acceleration  in  conformity 
with  CLMAX  limitations,  keeping 
Mach  number  constant 

Constant-g  climbing  or  diving 
right  turn  as  defined  by  bank 
angle  R0LLL  (see  Appendix  C) 

Straight  flight 


STRLVLd,  IAER0. 


TTHR,  LEVEL) 


Level  off  to  horizontal  position 


^These  subroutines  require  the  routine  MACHR8. 


Appendix  E 


DESCRIPXIOS  OF  INPUT  DATA* 


Data 

No. 

Program 

Variable 

Symbol 

Description 

1 

IRF 

Position  flag,  vehicle  1 

(Cartesian  coordinates) ;  IRF  =  0 

2 

R(l.l) 

*1 

x-coordinate,  vehicle  1 

J 

R(l,2) 

y-coordioate,  vehicle  1 

4 

R(l,3) 

2-coordinate,  vehicle  1 

5 

W0(1) 

«G1 

Initial  weight,  vehicle  1 

6 

IRF 

Position  flag,  vehicle  1 

(spherical  coordinates);  IRF  »  1, 

7 

R(1.4) 

^1 

Range  vector  magnitude,  vehicle  1 

8 

R(l,5) 

Range  angle  measured  in  horizontal 
plane,  vehicle  1  (see  Fig.  3) 

9 

R(l,6) 

‘Pi 

Range  angle  measured  in  vertical 
plane,  vehicle  1  (see  Fig.  3) 

10 

AREA(l) 

^1 

Reference  area,  vehicle  1 

11 

IVF 

Velocity  flag,  vehicle  1 

(Cartesian  coordinates) ;  IVF  *  0 

12 

V(l,l) 

x-component  of  velocity,  vehicle  1 

13 

V(l,2) 

^1 

y-component  of  velocity,  vehicle  1 

14 

V(l,3) 

z-component  of  velocity,  vehicle  1 

15 

ASMAX(l) 

Smaxl 

Maximum  lateral  acceleration  limit 
(structural),  vehicle  1 

16 

IVF 

Velocity  flag,  vehicle  1  (spher¬ 
ical  coordinates) ;  IVF  *  1 
Magnitude  only  (calls  AIM) ;  IVF  * 

17 

V(l,4) 

"l 

Velocity  vector  magnitude,  vehicle  1 

18 

V(l,5) 

\l 

Velocity  angle  measured  in  hori¬ 
zontal  plane,  vehicle  1  (see 

Fig.  4) 

19 

V(l,6) 

Flight-path  angle  measured  in 
vertical  plane,  vehicle  1  (see 

Fig.  4) 

no.) , 


Units  are  distance  (ft),  time  (sec),  velocity  (ft/sec  or  Mach 
acceleration  (g's),  angles  (deg),  weight  (lb),  area  (ft^). 
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Data 

Ho. 

Program 

Variable 

Symbol 

Description 

20 

JATM0S 

Flag  for  reading  in  velocity  mag¬ 
nitudes  (DATA  17,  33,  and  56) 
JATM^S  -  0  for  ft/sec 

JATMiJS  *  1  for  Hach  number 
(applies  to  all  vehicles) 

21 

TGUIDE(l) 

‘'gl 

Time  interval  that  vehicle  1  (if 
a  missile)  is  to  fly  ungulded 
after  launch 

22 

23 

24 

25 

WO(2) 

”02 

Initial  weight,  vehicle  2 

26 

AREA(2) 

"2 

Reference  area,  vehicle  2 

27 

ASMAX(2) 

^SBiax2 

Maximum  lateral  acceleration 
limit  (structural),  vehicle  2 

28 

ICAP 

Flag  to  indicate  initial-condition 
flight  status  of  vehicle  2 
(see  Appendix  H) 

29 

RMTMAX 

Maximum  flight  range,  vehicle  2 

30 

R(2,l) 

^2 

x-coordinate,  vehicle  2 

31 

R(2,2) 

^2 

y-coordinate,  vehicle  2 

32 

R(2,3) 

^2 

2-coordinate,  vehicle  2 

33 

V(2,4) 

Velocity  vector  magnitude,  ve¬ 
hicle  2 

34 

V(2.5) 

®V2 

Angle  measured  in  horizontal 
plane,  vehicle  2  (see  Fig.  4) 

35 

V(2,6) 

^2 

Flight  path  angle  measured  in 
vertical  plane,  vehicle  2 
(see  Fig.  4) 

36 

TBURNl 

ft 

First  stage  rocket  motor  burning 
time  (may  be  used  for  vehicle  1, 
2,  or  3) 

37 

TBURN2 

‘'b2 

Second  stage  rocket  motor  burnins 
time  (may  be  used  for  vehicle  1, 
2,  or  3) 

38 

KLAUN 

Decimal  fraction  of  vehicle's  maxi 

mum  range  at  which  it  is  to  be 
launchei  (may  be  used  for  vehicle 
1,  2,  or  3) 


-144- 


Data 

No. 

Program 

Variable 

Symbol 

Description 

39 

TGUIDE(2) 

’'g2 

Time  interval  that  vehicle  2 
(if  a  missile)  is  to  fly  un¬ 
guided  after  launch 

40 

IRT 

Position  flag,  vehicle  3 

(Cartesian  coordinates) ;  IRT  =  0 

41 

R(3,l) 

^3 

x-coordinate,  vehicle  3 

42 

P-(3.2) 

^3 

y-coordinate,  vehicle  3 

43 

R(3,3) 

^3 

2-coordinate,  vehicle  3 

44 

W0(3) 

M 

03 

Initial  weight,  vehicle  3 

45 

IRT 

Position  flag,  vehicle  3 

(spherical  coordinates) ;  IRT  =  1 

46 

R(3.4) 

^3 

Range  vector  magnitude,  vehicle  3 

W1 

1.(3. 5) 

®3 

R  .nge  angle  measured  in  horizontal 
plane,  vehicle  3  (see  Fig.  3) 

48 

R(3,6) 

Range  angle  measured  in  vertical 
plane,  vehicle  3  (see  Fig.  3) 

49 

AREA(3) 

^3 

Reference  area,  vehicle  3 

50 

IVT 

Velocity  flag,  vehicle  3 

(Cartesian  coordinates) ;  I'T  =  0 

51 

V(3,l) 

^3 

x-component  of  velocity,  vehicle  3 

52 

V(3.2) 

^3 

y-component  of  velocity,  vehicle  3 

53 

V(3,3) 

^3 

z-component  of  velocity,  vehicle  3 

54 

ASMAX(3) 

Smax3 

Maximum  lateral  acceleration  limit 
(structural) ,  vehicle  3 

55 

IVT 

Velocity  flag,  vehicle  3 

(spherical  coordinates) ;  IVT  »  1 
Magnitude  only  (calls  AIM) ;  IRF  * 

56 

V(3,4) 

Velocity-vector  magnitude,  vehicle  3 

57 

V(3,5) 

CD 

<  ' 

Velocity  angle  measured  in  hori¬ 
zontal  plane,  vehicle  3  (see 

Fig.  4) 

58 

V(3,6) 

^3 

Flight-path  angle  measured  in 
vertical  plane,  vehicle  3 
(see  Fig.  4) 

59 

TGUIDE(3) 

*^83 

Time  interval  that  vehicle  3  (if 
a  missile)  is  to  fly  ungulded 
after  launch 

4 
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Data 

No. 

Program 

Variable 

Syzsbol 

Description 

60 

(Not  used) 

61 

(Mot  used) 

62 

TIME 

t 

o 

Running  time  (when  entered  as  in¬ 
put  data,  it  is  equivalent  to 
initial  time,  i.e.,  time  at 
which  problem  is  to  start) 

63 

DTP0 

Print  interval,  i.e..,  time  incre¬ 
ment  for  printout 

64 

T0TAL 

Time  limit  placed  on  internal 
program  running  time  if  it  is 
not  to  be  teminated  after 
miss  calculation 

65 

ITAU(l) 

Flag  indicating  number  of  first- 
order  time  lags,,  vehicle  1 
(sea  Section  V) 

66 

TAU(1,1) 

"ll 

First  time  lag,  vehicle  1 

67 

TAU(1,2) 

Second  time  lag,  vehicle  1 

68 

TAUa,3) 

"13 

Third  time  lag,  vehicle  1 

69 

LAMDAO(l) 

^01 

Navigation  constant  for  guidance, 
vehicle  1 

70 

ITAU(2) 

Flag  Indicating  number  of  first- 
order  time  lags,  vehicle  2 
(see  ScctJon  V) 

71 

TAU(2,1) 

^21 

First  time  lag,  vehicle  2 

72 

TAU(2.2) 

■^22 

Second  time  lag,  vehicle  2 

73 

TAU(2,3) 

^23 

Third  time  lag,  vehicle  2 

74 

LAMDAO(Z) 

^02 

Navigation  constant  for  guidance, 
vehicle  2 

75 

ITAU(3) 

Flag  indicating  number  of  first- 
order  time  lags,  vehicle  3 
(see  Section  V) 

76 

TAU(3,1) 

"^31 

First  time  lag,  vehicle  3 

77 

TAU(3,2) 

^32 

Second  time  lag,  vehicle  3 

78 

TAU(3,3) 

’^33 

Third  time  lag,  vehicle  3 

79 

LAMDA0(3) 

^03 

Navigation  constant  for  guidance, 
vehicle  3 

— |4ib— 


Saca 

So* 

Probleni 

Variabla 

SjB&Ol 

Cescripcion 

SO 

CLSSAXCD 

Sjia»l 

SzxhBssa  aerodynassic  lift  co¬ 
efficient 

81 

coosmii} 

Zero  or  profile  drag  coefficient: 
Cjj  ,  vehicle  1 

82 

0 

Coefficient  used  vith  parabolic 
approxiiitation  for  drag  coef¬ 
ficient  as  a  function  of  lift 
coeff iciest,  vehicle  1 

83 

Slope  of  vs  a  curve,  assused 

to  be  a  ccnstast  (analytic 
functions),  vehicle  1 

Si 

ALPHAOd) 

“oi 

Total  zero-lift  angle  cf  attack, 
vehicle  1 

85 

CLS!AX<2) 

^Lasax2 

Kaxisssm  aerodynamic  lift  coef¬ 
ficient,  vehicle  2 

86 

CC0qlK(2> 

^302 

Zero  or  profile  drag  coefficient: 
Cjj  ,  vehicle  2 

37 

BC0!f(2) 

0 

Coefficient  used  with  parabolic 
approxixation  for  drag  coeffi¬ 
cient  as  a  function  of  lift 
coefficient,  vehicle  2 

88 

SL0?E(2) 

Slope  of  Cj^  vs  a  curve,  assused 
to  be  a  constant,  vehicle  2 

S9 

.AIPHA0(2) 

^02 

Total  zero-lift  angle  of  attack, 
vehicle  2 

90 

CLHftX{3) 

^Ls3x3 

Haxiausi  aerodynaoic  lift  coeffi¬ 
cient,  vehicle  3 

91 

CD0CftJ(3) 

^D03 

Zero  or  profile  drag  coefficient, 
vehicle  3 

92 

BCfHOy 

Coefficient  used  with  parabolic 
approximation  for  drag  coef- 
flcienc  as  a  function  of  lift 
coefficient,  vehicle  3 

93 

SL0PE(3) 

dC^/do3 

Slope  of  Cj^  vs  c  curve,  assumed 
to  be  a  constant,  vehicle  3 

94 

AL?aA0(3) 

®03 

Zcro-llft  angle  of  attack, 
vehicle  3 

* 

Note 
ore  used 

that  DATA  80-109 
instead  of  table 

are  reqoir 
values. 

ed  only  when  analytic  function': 

< 
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Data 

No.* 

Program 

Variable 

Symbo! 

95 

THC0N(1) 

T 

Ml 

96 

TABC0N(1) 

^abl 

97 

IMPLSE(l) 

98 

WBURN(l) 

W 

31 

99 

AB00ST(1) 

100 

THC0N(2) 

^M2 

101 

TABC0N(2) 

^ab2 

102 

IMPLSE(2) 

103 

WBURN(2) 

104 

AB00ST(2) 

^2 

105 

THC0N(3) 

T 

M3 

106 

TABC0N(3) 

^ab3 

107 

IMPLSE(3) 

108 

WBURN(3) 

109 

AB00ST(3) 

110 

KINTEG 

111 

ALT(i) 

h 

*Note  that  DATA  80-109  are 


are  used  Instead  of  Cable  values. 


Description 


Constant  for  military  thrust, 
vehicle  1 

Constant  for  afterburner  thrust, 
vehicle  1 

Specific  impulse  of  rocket  motor, 
vehicle  1 

Weight  at  rocket  motor  burnout, 
vehicle  1 

Boost  acceleration  (assumed  to  be 
a  constant),  vehicle  1 

Constant  for  military  thrust, 
vehicle  2 

Constant  for  afterburner  thrust, 
vehicle  2 

Specific  Impulse  of  rocket  motor, 
vehicle  2 

Weight  at  rocket  motor  burnout, 
vehicle  2 

Boost  acceleration  (assumed  to 
be  a  constant),  vehicle  2 

Constant  for  military  thrust, 
vehicle  3 

Constant  for  afterburner  thrust, 
vehicle  3 

Specific  impulse  of  rocket  motor, 
vehicle  3 

Weight  at  rocket  motor  burnout, 
vehicle  3 

Boost  acceleration  (assumed  to 
be  a  constant) ,  vehicle  3 

Flat-  or  round-earth  coordinate 
system  flag 

Flat-earth;  KINTEG  »  0 
Round-earth;  KINTEG  »  1 

Altitude,  vehicle  1 

(geocentric  coordinates) 


lived  only  when  analytic  functions 


1 
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Data 

No. 

Program 

Variable 

Symbol 

Description 

112 

L0NG(1) 

'‘1 

Longitude,  vehicle  1 

(geocentric  coordinates) 

113 

LAT(l) 

1 

Latitude,  vehicle  1 

(geocentric  coordinates) 

114 

INERF 

Flag  for  reading  vehicle  1 

velocity  in  local  or  inertial  system 
Local:  INERF  =  0 

Inertial:  INERF  =  1 

115 

IR0T8 

Flag  for  earth’s  rotation 

Nonrotating:  IR0T8  *  0 

Rotating:  IR0T8  *  1 

116 

ALT (3) 

hi 

Altitude,  vehicle  3 

(geocentric  coordinates) 

117 

L0NG(3) 

^3 

Longitude,  vehicle  3 

(geocentric  coordinates) 

118 

LAT(3) 

^3 

Latitude,  vehicle  3 

(geocentric  coordinates) 

119 

INERT 

Flag  for  reading  vehicle  3  veloc¬ 
ity  in  local  or  Inertial  system 
Local:  INERT  =  0 

Inertial;  INERT  =  1 

120 

L0NGO 

% 

Longitude  of  the  local  coordinate 
system  origin 

121 

LATO 

^0 

Latitude  of  the  local  coordinate 
system  origin 

122 

JINTEG 

Integration  flag 

Variable-step  Adams-Moulton: 

JINTEG  =  0 

Fixed-step  Runge-Kutta:  JINTEG  “  1 
Fixed-step  Adaras-Moulton:  JINTEG  =  2 
Variable  step  with  exact  print¬ 
out:  JINTEG  =  3  (see  Section  XIII) 

123 

ERTEST 

Number  of  significant  digits  of 
accuracy  required  (see  Section 

XIII)  of  numerical  integration 

124 

MINMR 

^IIN 

Vehicle  2  range  to  target  within 
which  program  will  automatically 
initiate  process  for  miss-dis¬ 
tance  computation 

125 

DVTH(l) 

Mi 

Assumed  error  in  6^  for  aiming 
vehicle  1 
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Kb. 

Prograa 

Variable 

Sysicol 

Oescripcion 

126 

D»153(2) 

Asscsed  error  irs  3^.,  for  aiaci^g 
vehicle  2 

127 

£K7H(3) 

3 

Assisaed  error  is  for  aianing 

vehicle  3 

128 

DVPHl<i) 

Asscned  error  is  ’f-g  aianigig 

vehicle  1  * 

129 

DVPHI(2) 

Assisted  error  is  for  aiarfsg 

vehicle  2 

130 

D¥Pai(3) 

4Y3 

Assisned  error  in  ^cr  adstiag 

vehicle  3 

131 

3ETA(i) 

“1 

Balliscic  coefficient,  vehicle  1 

132 

BETA(2) 

Saliistic  coefficient,  vehicle  2 

133 

Ei:TA(3) 

=3 

Ballistic  coefficient,  v^iicle  3 

134 

Hicinua  integration  step  size 
for  backup  to  ccstpute  miss 
distance 

135 

KHIS 

Mini  BUB  allowable  size  of  inte¬ 
gration  step  used  in  variable 
Adass-Houlton  node  integration 
(excluding  aiss  distance  com¬ 
putation) 

136 

DT0 

Initial  (starting)  value  of 
integration  step  size 

137 

H>IX 

Haximiia  specified  value  for  inte¬ 
gration  step  size  used  in  vari¬ 
able  Adaas-Houlton  code  inte¬ 
gration 

138 

ELEVMX(l) 

^maxl 

Maxicics  elevation  angle  limit, 
vehicle  1 

139 

ELEVMX(2) 

^nax2 

Maximun  elevation  angle  limit, 
vehicle  2 

.’40 

£LEVMX(3) 

^max3 

Maximum  elevation  angle  limit, 
vehicle  3 

141 

AZMAX(l) 

*^raaxl 

Maximum  azimuth  angle  limit, 
vehicle  1 

142 

AZMAX(2) 

Maximum  azimuth  angle  limit, 
vehicle  2 

143 

AZMAX(3) 

^raai:3 

Maximum  azimuth  angle  limit, 
vehicle  3 

Agyeadla  F 

Asaasi^Rjagc  /jgp  gao?BLSigs  •atsiEs 


Tables  fcraisbfiag  'vari^les  for  cotsgiuciag  die  aorodyaanlc  fusc- 
ciocs  aad  ibe  propolsioe  £±aracceristics  of  specific  aircraft,  iaclud- 
icg  rha  F-'10&  asd  F-1G5,  are  available.  These  are  to  be  used  «hes 
UlEi^  »  2  (tihea  che  cable  s<.ccioa  of  is  being  used). 

An  exasple  sec  of  tables  for  cn  aircraft  is  given  below;  ssets 
of  Cables  for  ocher  aircraft  staj  be  cosscrv£ced  following  a  siallar 
foncac. 

o  MavliwoB  life  coefficient  as  a  function  of  Mach  nmdier. 

o  Drag  coefficient  as  a  function  of  Mach  nusber  and  lift 
coefficient. 

o  ^3.  change  in  when  external  stores,  e.g. ,  fuel  tanks,  are 
jettisoned  (the  prograa  does  not  new  have  the  feature 
for  using  £C^,  but  cards  must  be  included  in  data  deck). 

o  Afterbumar  thrust  as  3  function  of  Mach  cuaber  and 
altitude. 

o  Fuel  flow  (afterburner  thrust)  as  a  function  of  Mach 
nusd^ar  and  altitude. 

o  Military  thrust  as  a  function  of  Mach  uuaber  and  altitude. 

Q  Fuel  flow  (ailltary  thrust)  as  a  function  of  Mach  number 
and  altitude. 

o  Placard  limit  (maxicinn  permissible  Mach  number)  as  z 
function  of  altltuie. 

o  Angle  of  attack  as  a  function  of  Mach  number  and  lift 
coefficient. 

See  Table  3  for  the  format  of  an  aerodynamic  table  using  the 
above-mentioned  functions.  Figure  35  shows  the  actual  data  deck  for 
the  F-104  aircraft. 
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Tablc  3 

FORMATS  FOR  AERODYNAMIC  TABLES 


Table 

No.  of 
Cards 

Card 

No, 

Quantity 

Format 

No ,  of 
Values 

c 

3 

1-3 

Mach  riumbei 

18F4.3 

i  42 

Lmsx 

i 

4-6 

^Lmax 

18F4.3 

42 

s 

1 

7 

Mach  number 

18F4.3 

16 

2 

8-9 

18F4.3 

32 

32 

10-41 

s 

18F4.4 

32<16 

*  c 

»j 

9 

42-50 

i6F4 .4 

9-16 

A3T^  and 

fuel  flow 

1 

51 

Mach  number 

18F4.3 

16 

2 

52-53 

Altitude 

14F5.0 

16 

32 

54-85 

ABT 

14F5.0 

16\16 

1 

32 

1 

86-117 

Fuel  Flow 

14F5,0 

16vl6 

Military 
thrust  and 

1 

fuel  flow 

1 

118 

Mach  number 

18F4.3 

16 

2  1 

119-1201 

Altitude 

14F5.0 

16 

1 

32 

121-152 1 

Thrust (military) 

14F5.0 

16x16 

32 

153-184 

Fuel  Flow 

14F5.0 

16x16 

Placard 

1 

1 

i 

limit 

1 

185 

Altitude 

14F5 .4  1 

14 

1 

186 

Mach  number 

14F5.4 

14 

Angle  of 

i 

1 

attack 

2 

187-188 

Mach  number 

12F6.2 

22 

2 

189-190 

Cl 

12F6.2 

13 

44 

191-234 

Alpha 

12F6.2 

13x22 

Afterburner  thrust. 
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C  «•••  F-104  TAai.ES  FCX  AERODYNAMICS 

C 

C  MACH.  COLIFTMAX 

200  300  600  TOG  100  830  9C0  930  970100010301100113012001300140013001600  1040001 

17001800190020002200  .070002 

1040003 

733  737  741  748  766  794  882101010131000  938  911  902  886  840  786  738  697  1040004 

634  639  612  391  374  1040003 

1040006 

C 

C  MACH.  CCLIFT.  COORAG 

000  800  850  900  923  930  973100010301100113012001400160018002400  1040007 

CCO  030  IOC  130  200  230  300  330  400  300  600  723  830  93010301130  1040008 

1040009 

161  166  173  197  228  236  310  372  460  669  94113131726208224622866  1040010 

1040011 

161  166  173  197  228  236  310  372  460  669  94113131726208224622866  1040012 

1040013 

165  17C  179  199  229  237  314  382  473  684  93813371737212223132931  1040014 

1040015 

I/O  176  183  205  234  278  336  403  497  719  98913731814220426283087  1040016 

1040017 

191  193  199  218  231  297  363  429  317  751102214121839223826933164  1040018 

1040019 

230  231  243  268  306  354  417  490  376  800108114821938234327833238  1040020 

1040021 

309  322  338  357  382  428  489  363  636  881116913762033244028783330  1040022 

1040023 

379  391  403  420  452  308  576  651  735  981127216912171239930673573  1040024 

1040025 

469  472  485  308  534  603  672  757  8681117140318332343280933263893  1040026 

1040027 

300  304  318  330  398  657  731  821  9191172148119372469295034814061  1040028 

1040029 

500  508  322  552  599  662  737  830  9411206134220292591309436494244  1040030 

1040031 

491  500  318  352  602  669  732  848  95812361381208826793212380C4441  1040032 

1040033 

449  457  478  511  972  644  741  852  9821287164521782807337940114705  1040034 

1C40C35 

407  415  44C  48l  543  627  739  857  9981327169322732988365844135259  1040036 

1040037 

38C  39'  <  I'l  47C  537  63C  747  87410281388178724213204393947695695  1040038 

1040039 

369  370  194  438  392  678  800  9801 140160021C030634400574C73209140  1040040 

1040041 

C 

C  CELTA  CGCRAC 

27  27  29  32  35  39  45  51  91  46  43  30  30  30  30  30  1040042 

1040043 

1040044 

1040045 

1040046 

1040047 

1040048 

1040049 

1040050 

C 

C  FACH.ALTITUOE. THRUST,  (AFTER-BURNER) 

Coiifinued 


Fig.  35  ■—  F-104  aerodynamic  fables 
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2S0  600  too  90010001iCC22001300l4001S001«001TOOltu01«002C002200  10400S! 

0  500010000190002«00c25c00a0000390efi400004s00050000950006000065000  1C4C0S2 

:«40093 

12400 {1200  4690  6290  7300  6350  5300  4200  3200  2300  1790  1900  1290  790  1040094 

1040095 

19350139001180010200  8450  7290  6000  4890  3800  2890  2200  1700  1390  800  1040056 

1040C97 

16750149501320S11J90  8700  8290  6800  9600  4390  3390  2500  190A  1900  890  1040098 

1040098 

17450157501400012200104C0  88C0  7330  6100  4750  3700  2790  2100  1690  1090  1040060 

1C40061 

1820016950148001300011250  8600  8100  6690  9250  4000  3090  2300  1750  1200  1040062 

1040063 

18790174C81969013300132CC10490  8800  7350  9800  4400  3350  2900  1800  1400  1040064 

1040065 

187501830816450148001320011350  8700  8100  6350  4850  3750  2800  2200  1690  1040066 

1040067 

1839018300173901510014150.239010690  8800  7100  94C0  42C0  3200  2400  1800  1040061 

1040068 

17900178001780016800151001330011600  8800  7800  6100  4700  3600  2700  2100  1040070 

1040071 

175CC17500175C0\750CU2C0143CC12S0010750  8600  6700  5200  3800  2800  2300  1040072 

1040073 

16800169001680017900U6C0193C01340011600  8250  7250  5700  4300  3250  2500  1040074 

1040075 

16450164S0t64S016400164C0198CC1420012400  8800  7750  600C  4700  3550  2750  1040C76 

1040077 

1S8C015800198C015800158CC158CC14550131001C500  8200  640C  4900  3880  2800  1040078 

1040079 

1470014700147C01470C147CC147CC14S501329010700  8350  6550  500P  3850  2900  1040080 

1040081 

i350C13SOOI350013500133C0135COl3S001330010700  8350  6500  5050  3900  2850  1040082 

1040083 

13300133001330013300133C0133C013308t33001C700  8350  6580  5050  39C0  2850  1040084 

1040085 

C 

C  MCH«6LTlTU0itFUfL  FLOW  <6FTE«»ltRN6«l 

28100251002140018100159001330011500  8400  7800  6800  5700  4700  4000  3700  1040086 

1040087 


355 

307 

265 

225 

151 

16C 

134 

110 

50 

74 

60 

50 

43 

38 

1040088 

1040085 

396 

347 

301 

257 

218 

184 

153 

127 

103 

83 

67 

55 

46 

40 

1040090 

1040091 

418 

3F« 

322 

276 

235 

158 

166 

137 

Ill 

85 

73 

58 

50 

42 

1040092 

1040093 

437 

354 

334 

288 

255 

214 

181 

150 

121 

57 

75 

64 

53 

44 

1040084 

1040095 

448 

417 

366 

320 

277 

233 

15? 

169 

132 

106 

85 

70 

57 

47 

1040096 

1040097 

447 

435 

388 

332 

257 

254 

214 

178 

144 

116 

93 

76 

62 

51 

1040096 

1040085 

435 

4.35 

410 

365 

315 

275 

233 

155 

157 

126 

101 

82 

67 

55 

1040100 

1040101 

418 

415 

418 

388 

34C 

255 

253 

211 

171 

137 

110 

89 

72 

58 

1040102 

1040103 

403 

403 

403 

438 

361 

316 

272 

232 

186 

148 

119 

86 

77 

63 

1040104 

1040105 

387 

397 

387 

444 

415 

336 

251 

245 

260 

160 

128 

103 

83 

68 

1040106 

1040107 

Continued 


Fig,  35  (continued) 


369 

369 

369 

434 

434 

390 

309 

266 

214 

171 

137 

110 

09 

72 

1040100 

1040109 

393 

393 

393 

410 

410 

410 

360 

202 

227 

101 

149 

116 

93 

76 

1040110 

1040111 

337 

337 

337 

370 

370 

370 

3:c 

320 

247 

196 

196 

129 

101 

02 

1040112 

1040113 

310 

310 

310 

330 

330 

330 

330 

334 

271 

219 

171 

137 

110 

00 

1040114 

1040119 

276 

276 

274 

276 

216 

276 

276 

346 

276 

219 

179 

141 

113 

90 

1040116 

1040117 

C 

C  XACHt  £LTIVU0E*  TMUST  (NILITARYI 

200  100  400  SOO  000  100  100  406100011001200110010002000  10401 IS 

0  300010000190002CCC0250CQ3000035000400004S00050000950006000065000  1040119 

1040120 


0140 

7240 

6200 

9300 

4500 

3800 

3100 

2590 

2000 

1940 

1140 

800 

400 

100 

1040121 

1040122 

8140 

7240 

6220 

5330 

4530 

3020 

3140 

2980 

2030 

1960 

1160 

810 

41C 

110 

1040123 

1040124 

8140 

7240 

6240 

9360 

4560 

3840 

3180 

2620 

2060 

1990 

1190 

820 

420 

120 

1040125 

1040126 

8140 

7250 

6270 

9390 

4590 

3860 

3220 

2660 

2090 

1620 

1220 

830 

430 

130 

1840127 

1040128 

8140 

7250 

6330 

9460 

4670 

3930 

3308 

2730 

2190 

1660 

1260 

86C 

460 

160 

1040129 

1040130 

8140 

7260 

6420 

5960 

477C 

4030 

3400 

2020 

2220 

1720 

1320 

900 

900 

200 

1040131 

1040132 

8150 

7370 

6990 

9760 

4970 

4230 

3900 

2980 

235C 

1820 

1390 

990 

570 

240 

1040133 

1040134 

8160 

7490 

6760 

5970 

9180 

4440 

3770 

3140 

2480 

1920 

1470 

1080 

650 

280 

1040139 

1040136 

8260 

7590 

6890 

6170 

9430 

4690 

4010 

3360 

2690 

2060 

1960 

1200 

780 

310 

1040137 

1040138 

836C 

7700 

7030 

6390 

5700 

4970 

4230 

..600 

2040 

2210 

1660 

1330 

920 

340 

1040139 

1040140 

8460 

7830 

7200 

6630 

6010 

9280 

4970 

3070 

3060 

2380 

1780 

1480 

1080 

380 

1040141 

1040142 

76C0 

7600 

7600 

7490 

6490 

5840 

9290 

4590 

3660 

2060 

2160 

1860 

1460 

900 

1040143 

1040144 

6700 

6700 

6700 

6700 

6530 

6300 

9000 

4820 

3070 

3020 

2470 

2190 

1790 

800 

1040145 

1040146 

9200 

5200 

9200 

5200 

5190 

9090 

5000 

4090 

3960 

3080 

29C0 

2200 

1800 

900 

1040147 

1040148 

1040149 

1040150 

1040191 

1040152 


C 

C  l«ACH,  ALTITUOEt  FUEL  FLOM  (^IllTARYI 


83CO 

73CO 

6200 

5100 

4100 

3200 

2600 

2200 

1800 

1500 

1000 

500 

1040153 

1040154 

8400 

7400 

6300 

5200 

4400 

3400 

2800 

2300 

2000 

1600 

1100 

600 

1040155 

1040156 

8700 

7  700 

65C0 

5500 

4600 

3700 

3100 

2500 

2100 

1700 

1200 

700 

1040197 

1040198 

9000 

8000 

6800 

5800 

4900 

4CC0 

3300 

2800 

2200 

1800 

1300 

800 

1040159 

1040160 

9500 

8300 

7100 

6100 

5200 

43C0 

3600 

3000 

2400 

1900 

1400 

900 

1040161 

1040162 

Conlinued 


Fig.  35  (continued) 
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10000  1700  7SC0  *500  9600  47C0  3«00  3200  2600  2000  1600  1200 
10500  9100  7900  6900  6000  5000  4300  3300  2t0C  2200  1700  1300 
11100  9600  0400  7300  6400  5400  4600  3700  3000  2300  1600  1400 
1160010000  «IOO  7700  7000  StCO  4900  4000  3200  2900  2000  1600 
121001C900  9200  8100  7500  63C0  5200  4300  3400  2700  2200  1800 


C  4LT1TUCE.  KACHHAX  IfLACARC  LI)*IT) 

0  50001000015C0C2CCa02S0CC3C00034000650GO 
113  125  138  153  168  182  191  200  200 

C 

C  6ACH,  CCLlFTf  AL8HA 


1040163 

1040164 

104C165 

1040166 

1040167 

1040168 

1040169 

1040170 

1040171 

1040172 

1040173 

1040174 

1040175 

1040176 

1040177 

1040178 

1040179 

1040180 

1040181 

1040182 

1040183 

1040184 


1040185 

1040186 


.0 

1.5 

.5 

1.6 

.6 

1.7 

.7 

1.8 

.8 

1.9 

.9 

2.0 

.95 

2.1 

1.0 

2.2 

1.1 

2.3 

1.2 

2.4 

1.3 

1.4 

1040187 

1040188 

.4 

l.l 

.0 

.1 

.2 

.3 

•  4 

.5 

•  6 

.7 

.8 

.9 

1.0 

1040189 

1040190 

9.0 

• 

.2 

2.0 

3.8 

5.6 

7.4 

9.2 

ll.O 

13.0 

15.3 

17.9 

1040191 

1040192 

9.0 

• 

.2 

2.0 

3.8 

5.6 

7.4 

9.2 

11.0 

13.0 

15.3 

17.9 

1040193 

1040194 

9.0 

.2 

2.0 

3.7 

5.4 

7.2 

9.0 

10.8 

12.9 

15.3 

18.0 

1040195 

1040196 

8.2 

- 

.2 

1.8 

3.5 

5.C 

6.7 

8.5 

10.3 

12.6 

15.2 

17.1 

1040197 

1040198 

7.4 

— 

.2 

1.6 

3.0 

4.4 

5.8 

7.4 

9.3 

11.6 

15.0 

19.5 

1040199 

104020C 

6.6 

21.8 

•• 

.2 

1.4 

2.8 

4.0 

5.2 

6.4 

7.8 

9.6 

11.9 

14.7 

18.0 

1040201 

1040202 

6.5 

14.0 

- 

.1 

1.9 

2.7 

4.0 

5.1 

6.3 

7.8 

9.2 

10.7 

12.0 

13.1 

1040203 

1040204 

5.8 

16.0 

.2 

1.7 

2.8 

4.1 

5.4 

6.7 

8.0 

9.4 

10.9 

12.5 

14.2 

104C205 

1040206 

4.4 

.8 

2.1 

3.6 

4.9 

6.3 

7.7 

9.2 

10.7 

12.5 

14.3 

16.1 

1040207 

1040208 

4.9 

.3 

1.6 

3.1 

4.8 

6.3 

7.8 

9.3 

10.8 

12.7 

14.5 

16.2 

1040209 

1040210 

6.9 

•• 

.5 

1.1 

2.7 

4.2 

6.0 

7.6 

9.2 

10.9 

12.8 

14.7 

16.6 

1040211 

1040212 

7.0 

• 

.6 

1.0 

2.1 

4.9 

6.2 

8.0 

9.7 

11.6 

13.5 

15.5 

17.6 

1040213 

1040214 

7.8 

• 

,6 

1.2 

3.0 

4.8 

6.6 

8.5 

10.4 

12.4 

14.5 

16.6 

1040215 

1040216 

■  8.6 

- 

.6 

1.4 

3.3 

5.2 

7.2 

9.2 

11.2 

13.3 

15.6 

>7.8 

1040217 
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Fig.  35  (continued) 
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12.1 

14.4 

16.7 

1040218 

1040219 

13.0 

19.4 

17.8 

1040220 

1040221 

13.9 

16.6 

19.7 

1040222 

1040223 

14.9 

17.8 

21.2 

1040224 

1040229 

16.0 
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02S3  CAROS 


Fig.  35  (continued) 


Appendix  G 
PROGRAM  SUBROUTINES 


The  subroutines  in  Table  4,  which  constitute  the  aaln  body  of 
TACTICS,  must  always  be  used  when  running  the  program;  all  other  sub¬ 
routines  are  optional. 


Table  4 

TACTICS  SUBROUTINES 


Subroutine 

Function 

aercJdn 

Computes  aerodynamic  variables  and  outputs  change  in 
velocity 

AIM* 

Aims  fighter  to  obtain  initial-condition  values  for 

V(I,5)  and  V(I,6)  angles 

ATM<?S 

Computes  model  atmosphere 

ATTITUD 

Computes  attitude  angles  for  the  vehicles 

AUXC0M 

Calls  output  at  specified  times 

C00RD 

Makes  rectangular  and  spherical  coordinate  transfor¬ 
mations 

CR0SS 

Computes  vector  cross  product 

DAUX 

Computes  derivatives  for  position  and  velocity 

DECRD 

Reads  initial-condition  data 

D0T 

Computes  vector  dot  product 

GE0FRC* 

Computes  force  for  geocentric  integration 

GE0CEN* 

Determines  earth  rotation  rate  and  finds  unit  geometric 
vectors 

INC0ND 

Reads,  computes,  and  prints  initial  conditions 

INITS 

Initializes  conditions  for  integration 

INTGRT 

Calls  integration  and  computes  new  velocities  and  ranges 

LAG 

Imposes  time  lag  (see  Section  V) 

If  storage  is  a  problem,  dummy  decks  may  be  substituted  for  these 
subroutines,  since  they  are  only  used  when  special  options  are  selected. 
(Since  the  routines  in  the  main  body  refer  to  the  special  options,  how¬ 
ever,  some  type  of  deck  must  always  be  used  to  represent  them.) 
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7:£>le  & 


Ssbreatiae  | 

Fucctim 

umt  1 

1 

Inpeses  aetradjsaacic  aad  strrartcsral  Cl^Ceral)  accelera- 
tiraa  coosttaiats  «hea  applicable  (see  Sectioo 

905 

Coopdtes  nacsitode  of  a  vector 

SiADP 

Orders  opesatloas 

Prints  ostpnt 

njtOD 

lessees  throttle  settisc  if  there  is  a  danger  of  ex- 
ceedias  Mach  naxber  (placard  liMt) 

SUkUS 

i  Competes  Mgalar  rates 

SriKE®  1 

1  Stores  posidoo  erd  veloclt^r  «hes  i^ing  restore  optioo 

XABD^  I 

1  loterpolatrs  aerodynamic  tables 

Ti5L£g  I  leaSc  aets^fs^t^cic  uoles  if  rated 
TK2ST  I  Cbai^ces  e&rratc  of  velocities 

Ttff#C3i*  I  Makes  local  asd  inertial  coordinate  transforaatioss 
VEICE  I  Covpctes  «ei|^t  of  vi^icles 

*If  stora'^  is  a  probles,  dnaov  seeks  asaj  be  substituted  for  these 
£:d>rootl3es,  since  the^  are  only  used  vhes  special  options  are  selected. 
(Since  the  roatiaes  la  tiie  main  bodj  refer  to  the  specie!  options^  hev- 
ever,  soae  type  of  deck,  cost  elvays  be  used  to  represent  then.) 
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Appendlx  H 
FORTRAN  FLAGS 


As  is  usually  Che  case  with  large  computer  programs,  TACTICS  uses 
Esany  FORTRAN  flags  for  internal  control  and  operation.  Although  in 
cost  ^ases  the  user  need  not  be  concert!.!  with  their  definition  nr 
functional  purpose,  certain  key  flags  require  explanation.  Some  of 
chose  flags  aay  be  set  in  P0LICY  or  by  an  initial-condition  data  value 
to  select  a  program  option;  others  are  important  because  they  control 
or  are  indicators  of  various  phases  of  program  operation.  The  follow¬ 
ing  is  an  explanatory  list  of  those  FORTRAN  flags  considered  to  be  of 
primary  importance: 

IMISS:  Whether  program  is  to  continue  after  finding  closest 
missile  approach. 

IMISS  ^  0,  the  program  stops  after  finding  miss  distance. 
IMISS  =  1,  the  program  continues  after  finding  missile 
miss;  flag  is  set  in  P0LICY. 

IMISS  <=  2,  the  program  has  determined  the  miss  distance 
and  is  ready  Co  continue;  flag  is  set  in  AUXC0M, 
and  is  normally  uted  in  P0LICY  as  a  criterion. 

IST0RE:  Whether  position  and  velocity  values  are  to  be  restored 
to  those  existing  at  launch  time. 

IST0RE  G,  the  positions  and  velocities  are  not  to  be 

A 

stored. 

IST0RE  «  1,  Che  store  option  is  to  be  used.  The  values 
of  position  and  velocity  are  stored  at  launch, 
and  the  program  returns  to  these  values  after 
computing  the  closest  point  of  approach;  flag 
is  set  in  P0LICY  and  used  in  MAIN  to  trigger 
ST0RE  routine. 


IST0RE  is  set  at  zero  in  ST0RE  routine  to  Indicate  that  ’’alues 
have  been  restored,  and  is  normally  used  in  P0LICT  as  a  criterion. 
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JP0L,  KP0L, 
LP0L,  MP0L, 
NP0L 

JINTEG: 


ISTART: 


ILAUN: 


To  be  used  in  P0LICY;  initially  set  equal  to  1  in  INC0NO. 

The  type  of  integration  to  be  used;  set  by  DATA  122,  If 

integration  is  to  be  changed  during  a  run,  JINTEG  may  be 

reset  in  P0LICY. 

JINTEG  =  0,  variable-step  Adams-Moulton. 

JINTEG  =  1,  fixed-step  Runge-Kutta. 

JINTEG  =  2,  fixed-step  Adams-Moulton, 

JINTEG  =  3,  variable-step  Adams-Moulton  with  exact 
printout. 

Used  within  program  to  determine  specific  events. 

ISTART  =  0,  the  start  of  the  program  (time  =  0.0). 

ISTART  =  1,  value  is  set  after  the  first  integration 

step,  and  remains  at  this  value  until  closest 
missile  approach  or  until  the  end  of  the  problem 
run;  flag  is  set  in  MAIN. 

ISTART  =  2,  the  program  is  backing  up  to  compute  the 

closest  missile  approach;  flag  is  set  in  INTGRT. 

If  the  program  is  to  continue  after  computing 
miss  distance,  ISTART  is  automatically  reset  to 
1  in  AUXC0M. 

Used  within  program  to  indicate  status  of  missile  launching. 

ILAUN  =  1,  the  missile  has  not  been  launched;  flag  is 
set  initially  in  INC0ND  and  is  reset  in  CAPFLT 
if  recall  or  restore  options  are  used. 

ILAUN  =  2,  launch  criteria  have  been  overshot  because  of 
an  excessive  step  size,  and  program  is  to  back 
up  to  approach  launch  on  a  smaller  fixed  step; 
flag  is  set  In  LAUNCH  (see  Section  XI). 

ILAUN  =  3,  the  missile  has  been  launched;  flag  is  set 

in  LAUNCH  and  is  usually  used  as  a  criterion  in 
P0LICY. 


s 


JATM0S: 


KINTEG: 


ICAP: 
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ILAUN  =  4,  the  progran  has  backed  up  to  the  step  before 
launch  and  will  approach  on  fixed  step  sire 
integration;  flag  is  set  in  INITS.  This  is  to 
keep  CAPFLT  from  resetting  ILAU?<  to  1,  vhich 
would  prevent  the  aissile  derivatives  frca  being 
integrated. 

Whether  the  velocities  of  the  fighter  and  target  are  to 
be  read  in  ft/sec  or  Mach  number  for  initial  data. 

JATM0S  =  0  (ft/sec) 

JATM0S  =  1  (Mach  number) 

Wtiether  a  flat  or  spherical  earth  gravitational  force 
representation  is  to  be  used  (see  Section  IV  and  Appen¬ 
dix  B). 

KINTEG  *  0,  the  flat-earth  option  is  to  be  used. 

KINTEG  =  1,  the  round-earth  option  is  to  be  used. 

This  flag  may  serve  several  different  functions,  all 
related  to  the  flight  status  of  vehicle  2.  Its  primary 
purpose  is  to  avcid  needless  computations  and  Integration 
of  the  equations  of  motion  for  vehicle  2  when  it  is  not 
being  used,  i.e.,  captive  flight  or  undefined  motion. 

The  ICAP  value  should  be  read  in  as  0,  1,  2,  or  3 
(DATA  28). 

ICAP  *  0,  computations  and  integration  will  occur  for 

vehicle  2.  If  a  nonzero  value  for  V(2,4)  (DATA  33) 
has  been  entered  as  an  initial-condition  value, 
the  ILAUN  =  3  flight  status  is  assumed  to  ►;xist. 

If  DATA  33  is  an  exact  zero  value,  the  IIAUN  1 
status  Is  assumed  and  computation  will  occur,  pie- 
sumably  in  anticipation  of  the  ground  launch  of 
vehicle  2  at  some  subsequent  time.  If  an  intf r- 
cept  problem  is  involved,  the  target  vehicle  will 
be  vehicle  3  for  miss -distance  computation. 

ICAP  =  1,  vehicle  2  is  in  captive  flight  on  vehicle  1; 
if  launch  occurs,  it  will  intercept  against 


iCAP 


vehicle  3.  2io  coscjutacioos  for  vehicle  2  ere 
performed  prior  tr,  reaching  criteria  for  launch. 
*  2,  vehicle  2  is  not  being  used;  it  will  aurcsat- 
ically  be  located  at  the  zero  origin  with  zero 
velocity  and  acceleration.  Ko  conputations  or 
integration  will  occur  for  vehicle  2. 

ICAP  =  3,  vehicle  2  is  in  captive  f-llght  on  vehicle  3; 
if  launch  'iccursa  it  will  intercept  against 
vehicle  1.  No  ccsputations  for  vehicle  2  ar» 
performed  prior  to  reaching  criteria  for  launch. 
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Appendix  I 
MODEL  ATMOSPHERE 


A  nodel  atnosphere  may  be  postulated  In  either  tabular  or  analytic 
form  depending  on  the  degree  of  realism  required  by  the  problem.  For 
most  Intercept  problems,  an  exponential  analytic  form  is  sufficient,  and 
no  provision  is  made  in  TACTICS  for  incorporating  tabular  values.  The 
air  density  ;  ,  pressure  p,  and  absolute  temperature  T  may  be  approxi¬ 
mated  as  a  function  of  altitude  z  from  the  following  expressions,  which 
are  consistent  with  these  given  for  the  ICAO  standard  atmosphere  (al¬ 
though  constants  are  not  represented  to  the  same  degree  of  significance) 
For  the  troposphere  (0  <  z  <  z^  (ft)) 


t  =  59  -  0.00357  z  (‘‘F) 


P 


.  S.?5(. 
0.00357  z\ 

518. A  / 


(140) 


p 


0.002378 


(i  0.00357 
\  ~  518.4  } 


5lug/f t^ 


For  the  stratosphere  (z  <  z  <  86,000  ft) 

s 


t  =  -67'‘F 


p  =  489.456  e 


-(z-z  )/h 

®  Ib/ft 


T 

h  =  h  ft 

s  a  T 


(141) 


_E_ 


p  =  slug/  ft' 


(4) 
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where 

2 

p  “  pressure  at  z  ft  (lb /ft  ) 

T  =  absolute  temperature  (°R)  *  t  +  45/9.4  (°F) 
t  =  teraperature  at  altitude  z  ft  (z  ^  z^)  (°F) 

0  =  air  density  at  z  ft 

=  pressure  head  at  sea  level  (ft) 
p  =  pressure  at  sea  level  (Ib/ft^) 

d 

=  absolute  temperature  at  sea  level  (  R) 

h  =  pressure  head  at  z  (ft) 
s  s  2 

p  =  pressure  at  z  (calculated  from  Eqs.  (141))  (Ib/ft  ) 
s  s 

=  absolute  temperature  at  altitude  z^  (‘^R) 

z  =  an  altitude  defining  the  beginning  of  the  stratosphere  and 
isothermal  conditions  (ft) 

The  above  relationships  are  applicable  for  the  so-called  standard 
atmosphere,  with  a  temperature  gradient  of  -3.57°F/1000  ft  in  the  tropo¬ 
sphere.  Assuming  the  atmosphere  can  be  characterized  by  isothermal 
conditions  above  an  altitude  of  35,300  ft  at  a  temperature  of  -67°F 
yields  the  follw/ing: 

h  =  27,650  ft  h  =  20,916  ft 

a  s’ 

p  =  2,116  Ib/ft^  p  =  489.456  Ih/ft^ 

3  S 

T  =  518.4  °R  T  =  392. 4°R 

a  s 

t  =  59°F  t  =  -67°F 

a  s 

z  =  35,300  ft 
s 

The  speed  of  sound,  s,  is  equal  toV^l.4  p/p  ft/sec.  Mach  number  is 
related  to  vehicle  spead  V  by  M  =  V/s. 


« 
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